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1.0  INTRODUCTION  AND  SUMMARY 


1 . 1  Background 

Future  DoD  large  space  system  concepts  such  as  High  Altitude  Large 
Optics  (HALO),  High  Energy  Laser  Optics  (HELO)  and  Millimeter  (MM)  Wave 
applications  have  stringent  1 i ne-of-si ght  (LOS)  pointing  and  jitter  perfor¬ 
mance  reguirements  which  cannot  readily  be  met  with  existing  technology. 

These  systems  of  interest  have  common  features  in  terms  of  their  large  size, 
extreme  complexity,  low  mass,  low  stiffness  and  precision  structural  toler¬ 
ances.  Under  the  Active  Control  of  Space  Structures  (ACOSS)  Program,  the 
Defense  Advanced  Research  Projects  Agency  (DARPA)  has  sponsored  a  number  of 
studies  for  the  development  of  a  unified  structural  dynamics  and  control 
technology  base  to  support  the  future  development  of  these  large  space  system 
missions.  The  major  emphasis  in  these  studies  in  the  past  has  been  placed 
on  generic  control  law  development  for  active  vibration  suppression.  However, 
the  emphasis  is  now  shifting  towards  hardware  development  and  experimental 
verification  of  the  technology. 

This  document  reports  the  results  of  one  of  these  studies,  namely 
ACOSS-14,  onducted  by  TRW  Space  and  Technology  Group  for  DARPA  from 
June  1981  through  October  1982.  The  study  addressed  two  areas  of  interest: 

•  PREMOD  Contract  -  the  definition,  design  and  plan  for  an  on-orbit. 
Shuttle-based  Proof-cf -Concept  demonstration 

•  MOD  Contract  -  stability  ensuring  designs  which  extend  to  treat 
noncol ocated  actuator-sensor  control  systems,  while  remaining 
insensitive  to  modal  truncation/inaccurate  structural  models. 

1 . 2  Summary  of  PREMOD  Contract 

One  of  the  challenges  in  active  control  of  large  space  structures  (LSS) 
is  to  assure  the  stability  and  performance  of  structural  control  system 
designs.  The  problem  stems  from  the  relation  between  modeling  information 
and  performance.  Generally,  the  control  system  must  make  up  for  increasing 
modeling  uncertainty  by  commen sura tel y  giving  up  performance.  In  the  systems 
of  interest,  the  reguirement  for  performance  in  the  presence  of  uncertainty 
is  carried  to  an  extreme.  First,  the  finite  element  methods  currently  used 
to  predict  structural  behavior  truncate  and  provide  poor  knowledge  of  the 
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structural  modal  frequencies,  shapes  and  damping.  Moreover,  they  do  not 
account  for  nonlinear  structural  behavior  effects  which  can  often  dominate 
behavior.  Second,  the  desired  performance  levels  are  sub  urad  for  precision 
pointing  and  sub  um  for  structural  control.  Under  DARPA  sponsorship,  it  is 
currently  believed  that  the  control  theory  required  to  simultaneously  satisfy 
these  requirements  has  been  developed  but  it  has  not  been  experimentally 
demonstrated.  Therefore,  it  has  become  evident  that  hardware  tests  are 
needed  to  verify  the  existing  results. 

The  objectives  of  the  ACOSS-14  PREMOD  study  were  to  define,  design 
and  plan  an  on-orbit  experiment  to  further  develop  DARPA' s  understanding 
of  generic  structural  dynamics  and  control  of  LSS.  The  primary  reason  for 
wanting  to  conduct  the  tests  on-orbit  is  the  realization  that  currently 
neither  analysis  nor  the  combination  of  analysis  and  ground  test  data  pre¬ 
dicts  the  on-orbit  beha' ior  of  a  structure  well  enough  to  establish  that 
the  desired  on-orbit  performance  goals  can  be  met.  Also,  for  large  and 
flexible  structures,  ground  tests  are  not  always  feasible.  Therefore, 
since  verifying  performance  is  of  utmost  importance  to  DARPA,  on-orbit  tests 
have  two  immediate  goals.  One  is  to  learn  how  to  make  accurate  structural 
models.  The  other  is  to  verify  that  the  performance  objectives  can  be  met. 

The  major  results  of  the  PREMOD  study  were:  the  identification  of  the 
set  of  experiments  needed  to  meet  the  modeling  and  performance  verification 
goals;  the  selection  and  design  of  a  simple,  moderately  priced  structure 
which  can  be  used  as  the  test  bed  for  executing  these  experiments;  and  the 
analytical  verification  that  the  chosen  structure  and  experiment  concept 
are  feasible  and  compatible. 

Figure  1-1  and  Table  1-1  show  and  summarize  the  key  features  of  the 
resulting  experiment  system.  The  distinguishing  characteristic  of  the  con¬ 
cept  is  its  capability  to  provide  data  necessary  for  modeling  of  future 
spacecraft  of  interest,  while  demonstrating  a  technology 

1.3  PREMOD  Contract  Conclusions/Recommendations 

This  study  shows  that  a  moderately  priced,  on-orbit  experiment  system 
can  be  designed  to  provide  the  data  necessary  to  design  future  surveillance 
and  weapon  systems.  The  key  ingredient  to  keeping  the  cost  down  is  an 
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Figure  1-1.  Baseline  Offset  Cassegrain  Telescope  Test  Structu 


Table  1-1.  Salient  Features  of  Baseline  Offset 

Cassegrain  Telescope  Experiment  System 

Provides  the  data  necessary  to  design  future  surveillance  and  weapon 
systems 

•  Modeling  verification/extrapolation  objectives  are  met  together 
with  a  technology  demonstration 

•  Maximum  use  of  existing  hardware  is  made 

•  Fits  comfortably  on  one  Shuttle  pallet 

•  Program  lasts  for  4  years 

•  Estimated  cost  is  $43. 5M  (1982  $) 

Key  Features  and  Capabilities: 

•  Ground  and  on-orbit  experiments  are  compatible  at  the  outset 

•  State-of-the-art  performance  is  demonstrated 

-  1000:1  vibration  reduction 

-  structural  control  to  2  um,  3a 

-  open/closed  loop  LOS  control  to  .8/. 2  yrad 

•  Stellar  and  earth  imaging  are  the  final,  on-orbit  performance 
measures 

•  Passive  and  active  control  technologies  are  integrated,  validated 
and  evaluated 

•  Vibration  suppression  and  isolation  are  demonstrated  together 
with  high  speed  optical  autoalignment  and  figure  control 

•  Actuator  and  sensor  technologies  are  integrated  and  demonstrated 
Other: 

•  Concept  is  dynamically,  electronically,  thermally  and  power 
compatible  with  Shuttle  -  it  is  safe 

•  Experiments  are  autonomous,  with  manual  override  from  the 
astronauts  and  ground 
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integrated,  ground  and  on-orbit  test  plan,  with  focused  objectives  at  each 
stage.  The  hardware  required  in  the  test  are  not  off-the-shelf  but,  in 
many  cases,  has  been  built  and  has  been  demonstrated  to  work.  The  experi¬ 
ment  plan  exploits  this  situation  to  permit  execution  in  four  years. 

Execution  of  the  experiment  program  is  recommended.  If  this  is  not 
possible,  the  goals  and  objectives  set  forth  in  this  study  merit  considera¬ 
tion.  In  particular,  modeling  verification/extrapolation  experiments  may 
not  be  cost  effective  to  do  for  their  own  sake.  The  cost  of  a  Shuttle 
flight  is  just  too  expensive.  It  therefore  appears  to  be  more  cost  effective 
to  do  the  experiments  together  with  a  technology  demonstration.  Also,  the 
integrated  ground  and  on-orbit  test  approach  here  is  noteworthy.  The  ground 
experiment  will  be  required  to  provide  the  link  between  analysis  and  ground 
test  data  and  also  provide  the  basis  for  future  extrapolation  to  on-orbit 
behavior  (once  the  extrapolation  model  is  derived  in  the  post  experiment 
analysis  of  the  on-orbit  test).  The  ground  experiment  should  therefore  be 
designed  compatible  with  the  on-orbit  experiment  at  the  outset.  This 
simplifies  model  comparisons,  ensures  power,  thermal  and  computational 
system  compatibility  with  true,  on-orbit  capability,  and  reduces  integrated 
cost.  Finally,  it  is  recommended  that,  if  funding  limitations  are  severe, 
a  low  cost  ground  experiment  incorporating  the  essential  features  described 
in  this  document  is  still  useful.  Excluding  certain  elements  of  the  optical 
train,  for  example,  will  reduce  cost  significantly.  This  approach  will 
allow  for  growth  as  funding  becomes  available;  it  serves  to  compare  analytical 
and  ground  based  structural  models  in  the  near  term;  and  permits  selection 
of  the  most  promising  technologies  to  be  demonstrated  on-orbit  from  the 
several  that  are  available. 

1.4  Summary  of  MOD  Contract 

As  a  result  of  an  earlier  study  by  TRW,  ACOSS-8,  a  methodology  for 
designing  stable  control  systems  for  inaccurately  modeled  structures  was 
developed.  This  stability  ensuring  design  methodology  is  based  on  the 
positivity  of  operators  and  has  the  advantage  that  it  does  not  rely  on 
modal  truncation,  it  has  only  low  sensitivity  to  the  exact  knowledge  of 
the  structural  model,  and  can  accommodate  for  nonlinear  effects.  Moreover, 
it  leads  to  time  and  frequency  domain  solutions  which  can  be  implemented 
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using  current  microprocessor  technology.  Its  main  disadvantage  was  that 
it  required  all  the  actuators  and  sensors  in  the  control  system  to  be  co¬ 
located.  In  a  number  of  applications  of  interest,  however,  colocation  is 
not  possible  and  further  development  was  required. 

The  main  objective  of  the  ACOSS-14  MOD  study  was  to  extend  the  stability 
ensuring  design  methodology  of  ACOSS-8  to  include  systems  not  permitting 
colocated  actuators  and  sensors. 

A  major  result  of  this  study  is  the  successful  generalization  of  the 
ACOSS-8  design  technique  to  include  noncolocated  LSS  systems  with  possibly 
unequal  numbers  of  actuators  and  sensors.  In  order  to  do  this,  the  Sector 
Theory,  from  which  the  positivity  approach  was  derived,  was  reexamined  and 
the  assumptions  requiring  colocation  removed.  The  generalized  method 
degenerates  exactly  to  the  colocated  design  approach  as  a  special  case. 

In  addition,  the  revised  stability  ensuring  design  condition  provides  a 
quantifiable  indication  of  the  trade  between  design  insensitivity  and  per¬ 
formance.  A  flow  chart  is  shown  in  Figure  1-2. 
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Figure  1-2.  Flow  Chart  of  Generalized  Stability  Ensuring 
Design  Methodology 
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The  generalized  methodology  was  then  illustrated  on  two  design  examples 
and  shown  to  work  well.  In  the  "Flat  Plate"  example,  the  damping  augmenta¬ 
tion  of  the  critical  modes  improved  significantly.  Moreover,  when  the 
sensitivity  of  the  noncolocated  design  was  compared  with  an  existing,  very 
insensitive  colocated  design,  the  results  were  comparable.  Computational 
complexity  was  also  comparable.  The  results  were  also  similar  for  the 
Draper  Model  #1  (Draper  Tetrahedron)  example.  The  damping  of  the  tetra¬ 
hedron  apex  more  than  doubled  over  an  existing  colocated  design  which 
already  met  design  goals.  The  design  sensitivity  and  computational  re¬ 
quirements  here  were  also  found  comparable  to  the  colocated  design.  The 
price  paid  in  both  examples,  over  the  colocated  design,  was  the  addition 
of  a  single  sensor — no  new  actuators  were  added. 

1.5  MOP  Contract  Conclusions/Recommendations 

The  generalized  design  procedure  has  been  shown  to  work  well.  In 
particular,  it  is  concluded  that  noncolocated  control  systems  can  be  de¬ 
signed  to  be  insensitive  with  respect  to  stability  (stability  robustness). 
Also,  noncolocated  designs  may  actually  be  desirable  from  a  performance 
standpoint,  if  it  is  the  only  way  to  measure  the  variables  that  must  be 
controlled.  We  had  expected  the  colocated  designs  to  have  better  stability 
and  performance  properties  because  of  their  ability  to  operate  at  a  higher 
loop  gain  without  sacrificing  stability.  However,  this  was  not  true  in 
the  examples  considered.  We  found  that  the  colocated  system  must  sometimes 
work  harder  to  predict  and  control  poorly  observable  signals.  This  means 
that  it  may  need  a  higher  bandwidth  and  higher  loop  gain.  The  noncolocated 
design,  on  the  other  hand,  was  found  to  afford  to  be  less  robust  at  the 
outset  because  it  did  not  need  as  high  a  bandwidth  and  loop  gain.  It 
inherently  got  a  better  representation  of  the  critical  signals  because  of 
the  designer's  freedom  to  place  a  sensor  at  critical  locations.  Hence, 
when  the  overall  compensation  required  to  meet  performance  goals  were  con¬ 
sidered,  the  noncolocated  design  was  better.  The  price  one  pays  for  this 
is  the  additional  sensor  and  the  loss  of  decentralization  in  the  control 
law  computations. 
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A  corollary  which  includes  the  above  observation  and  the  ACOSS-8  result 
is  that  there  exists  a  quantifiable  trade  between  stability  robustness  of 
a  design  and  its  performance:  one  is  obtained  at  the  expense  of  the  other. 
Model  uncertainty  of  what  one  is  measuring,  controlling  and  their  relation 
may  be  large,  but  the  accuracy  must  be  known  within  specifiable  limits. 

If  the  model  uncertainty  is  small,  robustness  may  be  traded  for  performance. 
The  degree  to  which  one  succeeds  in  designing  performance  will  depend  on 
the  degree  to  which  one  is  measuring  and  controlling  what  one  really  wants 
to  control  and  measure.  Given  two  equivalently  accurate  models  (e.g., 
one  colocated  and  one  noncol ocated) ,  both  have  the  same  stability  robust¬ 
ness  and  output  control  performance  limitations  but  the  one  that  measures 
and  controls  the  desired  control  variable  more  directly  (noise-free, 
observable  and  controllable)  will  be  superior. 

In  as  far  as  lessons  learned,  it  is  concluded,  and  justified  from  the 
theory  developed  here,  that  a  good  rule  of  thumb  is  to  split  the  design 
problem  into  two  pieces.  First,  the  frequency  dependent  loop  gain  is 
allocated  (loop  by  loop)  proportionally  to  performance  and  the  knowledge 
of  the  model,  and  inversely  proportional  to  the  uncertainty  and  noise. 

The  loop  gain  matrix  is  then  scaled  to  the  highest  gain  (norm)  permitted 
by  the  positivity/Sector  condition  ensuring  stability  robustness. 

Finally,  we  conclude  that  the  ACOSS-8  proposed  three-step,  stability 
ensuring  design  approach  for  LSS  systems  appears  to  be  the  correct  approach 
in  the  sense  that  it  continues  to  minimize  risk  with  the  potential  for 
maximizing  performance. 

1)  It  assures  stability  robustness  at  the  outset,  when  it  is 
required  the  most. 

2)  It  identifies  model  parameters  through  on-orbit  tests  whenever 
the  starting  model  accuracy  is  too  high  to  a  priori  design 
the  required  performance. 

3)  It  tunes  the  controller  to  maximize  performance. 

It  is  therefore  recommended  that  the  developed  design  approach  be 
studied  further  in  order  to  streamline  it  and  automate  it.  This  will 
facilitate  its  use  in  more  complex,  realistic  problems.  It  is  also  recom¬ 
mended  that  a  representative  design  be  evaluated  on  actual  hardware. 
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1.6  Organization  of  Report 

Sections  2  through  9  address  the  technical  work  performed  under  the 
ACOSS-14  PREMOO  contract.  It  covers  the  ACOSS-14  baseline  experiment 
from  the  initial  conceptual  and  trade  stages  to  its  analysis  and  evaluation. 
The  flight  experiment  program  plan  and  cost  estimate  for  the  baseline  system 
has  been  inserted  in  Appendices  A  and  B  as  requested  in  the  contract's 
statement  of  work. 

Sections  8  and  9  address  the  technical  work  performed  under  the  ACOSS-14 
MOD  contract.  Section  8  covers  the  theoretical  development  of  noncolocated 
control.  Section  9  covers  the  design  examples  and  their  evaluation. 

Quoted  references  are  listed  following  Section  9.  Appendix  C  describes 
the  Orbiter's  digital  autopilot  model,  as  used  in  this  study.  Appendix  D 
contains  a  published  paper  detailing  the  digital,  colocated  design  of  the 
Draper  Tetrahedron  used  in  Section  9. 
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2.0  IDENTIFICATION  AND  DEFINITION  OF  EXPERIMENTS 

Task  one  identified  and  defined  a  generic  set  of  POC  experiments  that 
would  be  useful  to  both  DARPA  and  NASA.  The  implied  joint  venture  was  con¬ 
sidered  as  a  way  of  meeting  technical  objectives  at  a  reduced  cost.  This 
possibility  was  researched  by  examining  NASA's  and  the  AF 1 s  "generic"  LSS 
missions  listed  in  References  [1,  2].  The  interests  and  requirements 
described  in  these  documents  were  analyzed  and  the  set  of  common  technologies 
that  need  to  be  demonstrated  to  ensure  mission  success  were  identified.  The 
underlying  premise  in  this  research  was  that  a  joint  venture  would  only 
be  feasible  if  both  parties' interests  were  satisfied.  In  order  to  meet 
this  objective,  the  common  technical  unknowns  and  uncertainties  in  the 
systems  of  interest  were  identified  first  [the  identified  technologies 
directly  determine  the  knowledge  that  must  be  gained  through  experimenta¬ 
tion].  The  methods  needed  to  gather  the  data  to  reduce  or  eliminate  the 
uncertainties  then  make  up  the  experiments. 

2. 1  Determination  of  Technologies  to  be  Demonstrated 

The  results  of  the  NASA  document  search  leads  one  to  perceive  that 
its  interests  are  ideally  driven  by  the  desire  to  explore  and  exploit 
space  and  the  sciences  it  induces.  More  practically  however,  it  is  evident 
that  these  objectives  are  modulated  by  politics,  DoD  and  Congressional 
funding.  Some  of  the  "generic"  representative  systems  that  one  car  identify 
are:  solar  power  stations,  platforms,  communications  systems,  permanent 
man-in-space  systems,  terrestial /extra  terrestial  research  and  resource 
exploration,  and  astrophysics  research.  In  these  categories  one  can 
further  identify  systems  which  are  LSS.  Therefore,  any  true  desire  to 
build  these  systems  induces  a  need  for  demonstrating  LSS  technologies. 

The  most  important  technologies  identified  are:  structural  modeling  and 
verification,  dynamic  interaction  and  stability  during  assembly  and  con¬ 
struction,  man-machine  "loop"  in  space,  passive  and  active  damping  aug¬ 
mentation  in  space,  distributed  actuation  and  sensing,  and  ground  predicted 
performance  verification. 


2-1 


A  similar  search  through  AF  and  DARPA  documents  resulted  in  the 
perception  that  they  are  driven  by  "threats"  to  the  security  of  the  U.S. 

In  particular,  DARPA' s  role  appears  to  be  coordination  and  support  of  the 
technology  development  of  advanced  communication,  surveillance  and  weapon 
systems.  In  some  cases,  the  systems  of  interest  to  DARPA  and  the  AF  are 
LSS  but  the  requirements  can  sometimes  be  very  different  from  NASA's  missions: 
DARPA' s  missions  will  not  only  push  the  state-of-the-art  in  materials  and 
structure  but  also  in  various  advanced  control  technologies.  The  most 
important  of  these  are:  structural  modeling  and  verification,  dynamic 
interaction  and  stability  of  dynamic  subsystem,  passive  and  active  damping 
augmentation  in  space,  distributed  actuation  and  sensing,  ground  predicted 
performance  verification,  and  advanced  control.  The  last  includes  preci¬ 
sion  pointing,  vibration  isolation,  figure  control,  slew  control,  and 
off-line  controller  tuning. 

The  set  of  common  technologies  that  must  be  demonstrated  by  DARPA  and 
NASA  are  therefore: 

•  structural  modeling  and  verification 

•  dynamic  interaction  and  stability 

•  active  and  passive  damping 

•  distributed  actuation/sensing 

•  predicted  performance  verification. 

The  objective  of  structural  modeling  and  verification  tests  should  be 
to  demonstrate  that  current  analytical  models  are  sufficient  to  predict, 
at  least  grossly,  structural  behavior  in  space.  Secondly,  the  demonstra¬ 
tion  should  provide  the  experimental  data  required  to  permit  the  confident 
use  of  models  to  predict  the  behavior  of  systems  that  have  not  been  built 
yet. 

The  dynamic  interaction  and  stability  within  the  system  is  also  important 
to  understand  and  demonstrate.  Classically,  weak  dynamic  interaction  is  re¬ 
duced  by  using  local  controls  with  large  loop  gains.  This  is  difficult  to 
do  with  flexible  structures,  however,  because  the  structure  theoretically 
strong-couples  the  system  at  resonance  and  stability  concerns  preclude  the 
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use  of  high  gains.  It  is  therefore  necessary  to  validate  the  current 
understanding  of  the  problem  and  verify  that  proposed  modern  control  con¬ 
cepts  are  a  solution. 

Next,  the  integration  of  active  and  passive  damping  into  a  structure 
must  be  demonstrated  in  space  if  it  is  to  be  a  viable  approach  to  solving 
various  LSS  problems.  Current  ignorance  of  the  relationships  between  ground 
and  on-orbit  passive  damping  is  currently  too  high  and  could  lead  later 
to  confusing  and  costly  trades  that  are  false. 

Distributed  sensing  and  actuation  are  also  commonly  accepted  concepts 
planned  for  solving  various  LSS  problems.  Yet,  because  of  the  large  band¬ 
width  and  dynamic  range  encountered  in  LSS,  the  integration  of  coarse  and 
vernier  devices  will  be  non-standard  and  must  be  val  i dated. 

Finally,  what  is  the  real  current  performance  1 'mi t  of  the  integrated 
LSS  technology?  Nobody  seems  to  know.  Moreover,  it  is  not  likely  that 
paper  analyses  are  going  to  be  credible  until  current  uncertainties  are 
resolved  in  an  integrated  test. 

2.2  Generic  Set  of  Experiments 

Two  generic  sets  of  experiments  have  been  identified  in  response  to 
the  technology  demonstrations  outlined  in  Section  2.1. 

The  primary  set  includes: 

•  identification  of  structural  dynamic  model 

a  identification  of  structural  damping 

t  passive  damping  augmentation 

0  dynamic  interaction  with  the  orbiter 

0  active  damping  augmentation/control  and  verification 
0  pointing  performance  verification. 

These  experiments  are  interrelated  and  often  one  experiment  is  the  logical 
step  to  the  next.  Most  of  the  tests  can  be  readily  conducted  on  the  ground. 
This  is  desirable  as  the  difference  between  ground  and  on-orbit  results, 
especially  for  performance  verification  and  determining  existing  and  aug¬ 
mented  damping,  is  what  one  needs  to  know  to  improve  modeling  abilities. 
Current  vibration  and  pointing  control  capability  is  also  to  be  verified. 
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The  secondary  set  includes 

•  figure  control  and  associated  performance  verification 

•  slew  maneuvers 

•  control  law  tuning 

and  have  been  chosen  to  add  realism.  The  figure  control  experiment  would 
show  the  interaction  of  the  vibration  control  system  with  the  figure  control 
system.  The  interaction  of  the  attitude  control  system  and  the  vibration 
control  system  during  a  slew  maneuver  adds  realism  to  the  acquisition  and 
tracking  problem.  A  control  law  tuning  experiment  demonstrates  indirect 
adaptive  control . 

2.3  Generic  Set  of  Procedures  for  the  Experiments 

2.3.1  Assumptions  about  the  Generic  Procedures 

The  generic  experiments  were  assumed  conducted  from  the  Space  Shuttle 
Orbiter  and  therefore  hypothesized  to  be  automated.  However,  human  assistance 
with  the  Remote  Manipulator  System  (RMS),  or  through  Extra  Vehicular  Activi¬ 
ties  (EVA),  were  not  initially  excluded  when  the  generic  procedures  were 
defined.  Some  other  assumptions  initially  used  were: 

i)  The  experiment  structure  can  be  deployed  incrementally  or 

sub-structures  can  be  uncaged  to  change  the  dynamic  character¬ 
istics  of  the  structure. 

ii)  The  experiment  structure's  critical  modes  are  lower  than  the 
Orbiter  digital  autopilot  (DAP)  so  that  the  structure  can  be 
excited  by  the  DAP. 

iii)  A  conservative,  active  controller  is  used  to  stabi'Hze  the 
structure  and  to  damp  out  vibrations  between  experiments. 

iv)  The  Orbiter  DAP  signals  are  available  to  the  LSS  control 
system. 

v)  Actuators  installed  on  the  structure  are  used  to  generate 
disturbance  for  evaluating  control  system  disturbance  re¬ 
jection. 

vi )  The  controller  performance  measurement  is  sensed  using  the 

1 ine-of-sight  error  from  a  dummy  payload  to  the  Orbiter  cargo 
bay  or  to  an  inertial  reference  point.  This  information  did 
not  necessarily  have  to  be  available  to  the  control  system. 
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and 


vii)  All  the  actuator/sensor  data,  including  DAP  signals,  are 
recorded  on  tape  for  post-flight  analysis. 

2.3.2  Test  Procedure  Concepts 

Conceptual  procedures  to  the  experiments  identified  in  Section  2.2 
were  derived  next.  The  results  are  summarized  in  Table  2-1.  Listed  are 
the  candidate  experiments,  their  procedures,  and  required  ground/on-orbit 
activities.  The  table  is  subsequently  used  in  this  report  to  guide  the 
selection  of  a  suitable  test  structure  and  as  the  framework  for  the  detailed 
experiment  procedures  described  in  Appendix  A  for  the  baseline  test  structure. 

Referring  to  Table  2-1,  the  experiments  described  follow  the  ACOSS-8 
approach  to  LSS  control:  first,  a  robust,  stable  control  system  is  designed; 
then,  after  one  identifies  the  on-orbit  dynamical  behavior  more  accurately, 
the  control  system  is  tuned  to  obtain  the  best  achievable  performance.  From 
this  approach,  the  following  experiments  are  derived: 

1)  Identification  of  Structure  Dynamic  Model  -  A  more  accurate 
structural  dynamic  model  is  obtained  from  actuator/sensor 
measurements  by  using  parameter  estimation  techniques.  The 
preferred  techniques  are  the  maximum  likelihood  method  and 
the  transfer  function  method. 

2)  Structural  Damping  -  Ground  measurement  of  structural^  damping 
does  not  accurately  estimate  the  on-orbit  damping.  Thus, 
structural  damping  is  determined  on-orbit  to  ensure  the  control 
system  stability  and  performance.  Damping  is  measured  from 
free  vibration  decay  or  from  frequency  of  peak  coincident  re¬ 
sponse. 

3)  Component  and  Substructure  Testing  -  One  of  the  lessons  that 
will  be  learned  from  the  experiment  is  how  one  can  predict  the 
behavior  of  a  larger  structure  by  testing  an  individual  component 
or  substructure.  The  component  and  substructure  can  be  tested 
using  the  same  technique  as  the  full  scale  structure. 

4)  Passive  Damping  Augmentation  -  Effective  use  of  passive  damping 
augmentation  concepts  depends  on  their  being  designed  into  the 
structure  rather  than  added  as  an  after-thought.  The  design 
process  includes  detailed  modal  analysis  of  the  entire  structure 
to  find  those  members  possessing  significant  amounts  of  elastic 
energy  and  therefore  suitable  for  damping  augmentation.  .  The 
selection  of  damping  material  should  consider  the  material ^1 oss 
factor  as  a  function  of  frequency  and  temperature  range.  The 
experiment  structure  can  have  cageable/deployable  appendages 
which  have  the  ability  to  change  the  structure's  modal  char¬ 
acteristics  and  damping. 


2-5 


Table  2-1.  Generic  Experiment  Procedures 


Test  Objective 

Procedure 

EVA/RMS  Activities 

Post-Flight  Analysis 

Remarks 

1 )  Identification 
of  structure! 
dynamic  model 

e  Excite  the  experiment  structure  with. 

-  Actuetion  on  the  test  structure 

-  Orbiter  thrusters 

e  Aunt  the  deployment 
of  experiment 
structure: 

-  Change  the  charac¬ 
teristics  of  the 
experiment 
structure 

a  Maximum  likelihood  parameter 
estimations: 

-  Physical  parameters 

-  Modal  parameters 

a  Transfer  function^ frequency 
response  method 

e  Experiments  we  conducted 
during  and  after  the 
deployment  of  ten  structure 

e  Structural  model  may 
include  Orbiter 

2)  Structural 

Damping 

e  Excite  the  structure  with. 

-  Sinusoidal  signals 

-  frequency  sweep  sinusoidal  signals 

Same  as  above 

a  Determine  damping  from 
free  vibration  decay 

a  Determine  damping  from 
frequency  of  paek  coincident 

response 

e  May  require  long  time  of  no 
OrtMter  thruster  activity 

e  Relatively  "clean"  data  are 
required 

3)  umponent 
•sting 

a  Sat  up  special  test  bench 

a  Test  the  dynamic  properties  of  individual 
component  'truss  cumper  actuator  sensor) 
by  observing  the  response  to  an  input  signal 

e  Place  test  article  on 
test  bench 

e  Determine  component 
physical  properties  (stiffness, 
damping)  end  dynamic 
properties  (frequency  response) 

e  Test  can  be  conducted 
independently 

4 )  Substructure 
Testing 

e  Test  the  dynamic  properties  of  substructure 

Same  as  above 

a  Determine  the  dynamic 
properties  of  each  individual 
substructure  end  combine  the 
results  to  predict  total 
structure  response 

e  Provides  data  on  modal 
synthesis  techniques 

5 )  Pstsiv#  Dempmg 

Augmentation 

a  Excite  the  structure  before  and  after  the 
replacement  of  en  existing  truss  with  a 
damping-enhanced  one 

e  Excite  the  structure  before  and  after  the 
addition  of  passive  damper 

e  Replace  or  install 
passive  damping 
elements 

e  Determine  structure!  damping 
before  end  after  passive 
damping  augmentation 

•  Evaluate  the  effects  of 
damping-enhanced  devices 

6}  Control  lew  and 
pointing 
performance 
verification 

•  )  Select  controller 
e  Open-loop  response 

-  Increase  input  signal  mayutude  to  test 
linearity 

•  Cloaed-loop  response 

-  Increase  delay  to  test  phase  margin 

-  Increase  gem  to  tett  gem  margin 

•  Simulate  disturbance 

-  Measure  disturbance  rejection 

-  Tune  controller 

•  Cloaed-loop  responee  for  tuned  controller 

•  Simulate  actuator  /sensor  malfunction 

b)  Change  structure  characteristics,  repeat  Hem  e) 

c)  Assess  pointing  performance 

None 

e  Compute  open-loop  response 

e  Check  lineenty 

e  Compute  cloaed-loop  response 

-  Compere  with  predictions 
from  open-loop  response 

-  Verify  phase/gsm  margins 

e  Estimate  cloaed-loop  damping 

e  Com  put#  disturbance  refection 

e  LOS  error  from  test 

Itructurs  to  Orbitsr  cargo 
bay  can  be  used  et 
performance  index 

e  Disturbenca  car  be  simu¬ 
lated  by  setueters  installed 
on  the  structure 

e  Disturbenca  signal  can  be 
sinusoidal  or  band-iimlied 
random  noise 

7 )  Dynamic 
Interaction 
with  the 

OrtMter 

•  Firing  orbiter  thrusters  end  observing  the 
structure  responee  with  active  controller  on/off 

•  Using  inertial-  .ype  actuators  to  excite  the 
total  structufJ  end  observing  the  response  with 
combmeticns  of  Orbiter  OAP  on/off  end 
active  controller  on/off 

•  Repeat,  after  changing  the  structure 
characteristics 

•  Change  experiment 
structure 
characteristics 

a  Orbiter  fuel  consumption 

e  Dynamic  coupling  between 
Orbiter  end  experiment 
structure 

8 )  Comparison  of 
Actuators 
(relative  vs 
inertial) 

•  Evaluate  controller  performance  using  relative 
type  actuators  (linear  force  actuators) 

e  Evaluate  controller  performance  using  inertial 
type  actuators  (reaction  wheels  end/or  proof 
mess) 

None 

•  Estimate  cl osed-looQ  damp*  .g 

e  Evaluate  disturbance  faction 
(control  per  forms  nee) 

e  Locations  end  characteris¬ 
tics  of  these  two  types  of 
actuators  ere  different,  only 
qualitative  comparisons  can 

be  made 

e  Reletive  ectuetor  will  not 
effect  the  rigid  body  motion 

a  Rigid  body  will  respond  to 
inertial  actual  v.  This  is 
important  when  LOS  point¬ 
ing  it  to  be  controlled 

t 
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5)  Control  Law  and  Pointing  Performance  Verification  -  Control 
law  characteristics  that  should  always  be  tested  are  stability, 
robustness,  and  meeting  of  performance  design  goals.  Robust¬ 
ness  and  stability  can  be  verified  by  varying  control  parameters 
until  the  system  is  on  the  verge  of  instability,  provided,  of 
course,  that  this  can  be  done  safely.  Damping  design  goals  can 
be  directly  measured  from  sensed  responses.  Pointing  control 
performance  can  be  measured  from  the  LOS  error  between  a  dummy 
payload  to  Orbiter  cargo  bay.  This  information  is  intended  for 
assessing  pointing  performance  only;  it  is  not  available  to  the 
control  system.  The  vibration  control  or  disturbance  rejection 
can  be  tested  by  installing  special  actuators  on  the  structure 
that  have  output  characteristics  resembling  the  on-orbit  dis¬ 
turbance  source,  e.g.,  cryocoolers. 

6)  Dynamic  Interaction  with  the  Orbiter  -  The  dynamic  interaction 
with  the  Orbiter  must  be  well  understood  for  the  following 
reasons:  1)  safety  of  the  Orbiter  and  its  crew,  2)  excessive 
fuel  consumption  of  the  Orbiter  DAP,  3)  effects  of  the  dynamic 
interaction  on  the  test  structure  and  the  outcome  of  the  ex¬ 
periments,  and  4)  being  attached  to  the  Orbiter  will  be  the 
normal  mode  of  operation  for  some  future  LSS  missions.  Of 
particular  interest  is  the  relations  among  VRCS  thruster  firing, 
Orbiter  DAP,  and  LSS  control  system. 

7)  Comparison  of  Actuators  -  There  are  two  types  of  actuators  that 
can  be  used  for  LSS  control:  relative  (e.g.,  linear  force 
actuator)  and  inertial  (e.g.,  reaction  wheel,  proof  mass). 

The  major  difference  between  these  two  types  of  actuators  in 
theory  is  that  the  relative  type  actuator  will  not  affect  the 
center  of  mass  rigid  body  motion  while  the  rigid  body  will 
respond  to  inertial  actuators. 

Tests  of  components,  substructures,  and  actuators  have  been  included 
in  the  list  because  of  their  unique  modeling  importance.  It  is  precisely 
the  relation  between  the  characteristics  of  the  components,  the  substructure, 
and  the  total  structure  that  is  needed  to  be  able  to  develop  analytical 
models  in  LSS.  Also,  some  of  these  experiments  can  be  performed  on  the 
ground  or  in  orbit  but  the  results  can  differ  significantly  in  some  cases 
(e.g.,  damping).  Therefore,  these  tests  should  be  performed  in  both  loca¬ 
tions  so  that  analytical  prediction  methods  can  be  improved  and  future  on- 
orbit  predictions  based  on  ground  data  can  be  more  accurate. 
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3.0  SELECTION  OF  SUITABLE  TEST  STRUCTURE 

Having  defined  the  generic  experiments  to  be  conducted,  the  next  task 
in  the  study  was  to  select  a  structure  to  accommodate  them.  Existing, 
planned  and  new  designs  were  considered.  The  approach  taken  was  to  first 
obtain  a  list  of  the  planned  Shuttle  payloads.  This  list  was  compiled 
from  information  contained  in  NASA  STAR  from  1977  to  1980,  from  NASA 
literature  search  #45835  ("Space  Shuttle  Payloads")  [1,  3,  4,  20,  21] 
and  from  various  DoD  documents  [2,  4].  The  last  was  necessarily  limited 
due  to  "need  to  know"  constraints.  From  this  base  list,  structures  which 
could  accommodate  the  experiments  were  next  identified.  The  criteria 
was  to  select  systems  which  are  consistent  with  the  NASA  and  DoD  Technology 
Models  [1,  2]  and  which  have  reasonable  time  frames  and  level  of  commit¬ 
ment.  At  the  e-d,  new  possibilities  consistent  with  NASA  and  DoD  were 

added  and  considered. 


3.1  NASA  Missions 


3.1.1  Mission  Definition  Process  for  NASA 

NASA  missions  were  examined  first.  As  one  might  have  guessed,  the 
"wish  list"  was  very  large  and  it  became  important  to  try  to  rank  missions 
according  to  credibility.  The  block  diagram  in  Figure  3-1  was  found  useful 
in  performing  this  operation.  It  models  the  NASA  mission  definition  process, 
It  is  based  on  sheer  perception  and  consistency  with  available  data. 


Referring  to  Figure  3-1,  the  process  starts  at  NASA  and  DoD  think  tanks 
Their  organizations  look  at  current  objectives,  what  has  been  successful 
in  the  past,  and  speculate  about  the  future.  The  output  of  the  think  tanks 
are  termed  unvalidated  missions  and  there  is  no  time  frame  associated  with 
them.  Then,  depending  on  the  "ideas",  unvalidated  models  are  discarded  or 
validated  by  NASA  major  offices.  The  offices  also  determine,  depending  on 
usefulness  and  cost,  whether  the  validated  models  are  0-models  or  C-models. 
0-models  are  potential/speculative  missions  with  possible  starts  beyond 
10  years.  C-models  are  candidates  for  possible  start  within  10  years. 

The  validated  models,  together  with  a  speculated  DoD  coordinating  input, 
are  next  processed  by  NASA  to  determine  P-models.  These  are  missions 
planned  to  start  within  5  years.  The  final  NASA  process  takes  P-models, 
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SPECULATIVE  DOD  INTERACTION  PROCESS 


0  MODELS:  POTENTIAL/SPECULATIVE,  POSSIBLE  START  >10  YEARS 


DoD  coordination  and  congressional  approval /funds  into  A-models.  .hese 
are  approved  and  funded  systems,  scheduled  for  launch  within  5  years. 

3.1.2  Feasible  NASA  Missions 

A  preliminary  set  of  criteria  used  to  “sift"  the  NASA  models  of  interest 
was  defined  to  determine  which  missions  are  feasible.  It  was  determined  that 
politically,  only  USA,  NASA  missions  could  be  realistically  considered.  This 
would  exclude,  for  example,  an  International  Solar  Polar  Mission  (ISPM)  like 
mission.  The  system  would  also  have  to  be  Shuttle  launched  and  earth  orbit¬ 
ing  to  permit  monitoring  and  exclude  planetary  spacecraft  which  could  not 
be  retrieved.  Finally,  the  structure  in  the  mission  would  have  to  have 
the  LSS-'i  j  features  of  the  "true"  system  of  interest. 

Validated  Systems.  The  set  of  approved  and  planned  NASA  missions  that 
satisfy  the  feasibility  criteria  above  is  shown  on  Table  3-1.  The  raw  data 
is  taken  from  [1].  The  alpha-numeric  code  designates  the  NASA  mission 
number.  It  is  interesting  that  only  one  mission,  the  Solar  Electric  Power 
Array  test  flight,  had  been  approved  (is  funded  for  launch)  at  the  time 
this  research  was  made,  1981.*  The  remaining  missions  were  "planned". 

The  more  interesting  missions  identified  are  the  radiometers  ( RI 12 ,  RI13, 
and  EI67) ,  LSST  (Large  Space  Structure  Technology  program),  and  30/20  GHz. 
However,  none  of  these  are  expected  to  require  ultra-precision  structural 

control . 

Continuing  the  search  through  the  remaining  validated,  feasible  missions, 
the  systems  become  more  and  more  LSS-like.  Table  3-2  shows  the  "candidate 
and  "opportunity"  missions.  It  is  evident  that  these  systems  are  more 
speculative  than  the  ones  in  Table  3-1.  Nevertheless,  it  is  important 
to  observe  that  the  Space  Operation  Center,  which  had  no  set  date  in  1980, 
is  the  current  Space  Station  for  which  NASA  has  recently  initiated  Phase  A 
study.  Hence,  abrupt  status  changes  are  very  possible. 


*  Many  of  the  NASA  missions  described  have  since  been  delayed  or  eliminated. 
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Table  3-1.  Feasible  NASA  Missions  (1980  Data) 


0112 

SEP  Solar  Array  Shuttle  Flight 

1983 

Launch 

Planned: 

L2 

Space  Science  Platform 

1982 

Start 

R 1 12 

Passive  Microwave  Radiometer 

Mid  1980's 

R 1 13 

Active  Microwave  Sensor 

Mid  1980's 

E 167 

Pushbroom  Radiometer 

<1985 

Start 

C4 

30/20  GHz  Antenna  Wideband  Communications 
Satellite  Program 

1982 

Start 

C5 

Geostationary  Platform  Demonstration 

1983 

Start 

T9 

Solar  Electric  Propulsion  Stage  (SEPS) 

1981 

Start 

U4 

Tethered  Satellite  System 

1981 

Start 

U5 

LSST 

1981 

Start 

U6 

Power  Extension  Package  (PEP) 

1981 

Start 

U7 

Science  and  Applications  Space  Platform 

1982 

Start 

118 

25  KW  Power  System 

1983 

Start 

U9 

Material  Experiment  Carrier 

1983 

Start 

U10 

Material  Experiment  Carrier  II 

1986 

Start 

Ull 

Large  Power  Module 

1990 

Start 
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Table  3-2.  Speculative,  Feasible  NASA  Missions  (1980  Data) 


Candidates : 

A15 
A16 
A 17 
S12 
R8 


Very  Long  Baseline  Radio  Interferometer 
Orbiting  Submillimeter  Telescope 
Large  Ambient  Deployable  IR  Telescope 
Solar  Terrestial  Observatory 
Soil  Moisture  Research  and  Assessment 


<1990 

1985-90 

1985-90 

1985 

Mid  1980 's 


Opportunity : 

A18 

A19 

A20 

A21 

A22 

C6 

C8 

T14 

U13 

U14 


Infrared  Interferometer 
Gravity  Wave  Interferometer 
Coherent  Optical  System  of  Modular  Imaging 
Collectors 

Large  Optical /UV  Telescope,  Diffraction  Limited 
100-Meter  Thinned  Aperture  Telescope 
Multi-Service  Thin  Route  Narrowband  Program 
Orbiting  Deep  Space  Relay  Station 
Advanced  Electric  Propulsion 
Materials  Experimentation  Module 
Space  Operation  Center 


>1990 

>1990 

>1990 

>1990 

>1990 

1984 

>1990 

1990 

1990 

No  Date  Set 


Start 

Start 

Start 

Start 

Start 


Start 

Start 

Start 

Start 

Start 

Start 

Start 

Start 

Start 
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The  final  set  of  possibly  feasible,  validated  missions  are  shown  in 
Table  3-3.  The  available  descriptions  of  these  systems  were  too  ambiguous 
to  determine  whether  they  were  feasible  or  not.  The  list  consists  primarily 
of  instruments  (as  opposed  to  spacecraft)  and  a  few  Earth  resources  space¬ 
craft. 

'-Invalidated  Systems.  Unvalidated  missions  were  also  considered  in  the 
search  for  a  suitable  test  structure.  These  missions  are  primarily  proposed 
and  evaluated  by  "think  tank"  organizations  through  various  funded  studies. 

Some  of  the  more  extensive  studies  include:  "Advanced  Space  Concepts  and  . 
jheir  Orbital  Support  Needs"  [4],  "Space  Resources  and  Space  Settlements"  [5], 
"Space  Horizons"  [6],  and '"Post  LANDSAT-D,  Advanced  Concept  Evaluation"  [7]. 

It  is  evident  from  the  "NASA  Technology  Models"  [  1]  document  that  NASA  used 
many  of  the  missions  described  in  these  documents,  but  only  after  they  had 
been  further  studied  by  two  think  tanks:  General  Research  Corporation  and 
Planning  Research  Corporation.  The  number  of  identified  unvalidated  systems 
is  also  rather  large  (<v.  50)  but  NASA  OAST  goes  through  a  ranking  exercise 
to  determine  which  should  become  validated  missions.  The  missions  listed 
in  Table  3-4  meet  the  feasibility  criteria  and  are  in  the  top  30  of  the 
rankings.  The  Permanent  Orbiting  Outpost  is  number  one  and  the  "Geosyn¬ 
chronous"  spacecraft  are  numbers  three  and  four. 

3.1.3  Most  Promising  NASA,  Feasible  Systems 

All  of  the  feasible  NASA  systems  described  above  were  carefully  examined 
and  only  the  most  promising  were  selected  for  further  evaluation.  The  following 
twelve  candidates  were  judged  to  be  the  most  promising  with  respect  to  near- 

term  starts  (<  5  years): 

0112  SEP  Solar  Array  Shuttle  Flight 

U5  LSST 

T9  SEPS 

U6  PEP 

U8  25  KW 

C4  30/20  GHz 

RI12  Passive  Microwave  Radiometer 

RI13  Active  Microwave  Radiometer 

/\  16  Orbiting  submillimeter  Telescope 

A 1 7  Large  Ambient  Deployable  IR  Telescope 
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Table  3-3.  Possibly*  Feasible  NASA  Missions 


RI6 

pi anned : 

A14 

R 1 10 

Rill 

R 1 14 

R 1 1 5 

R 1 16 

E6 

E7 

E8 

E9 

E10 

EI25 

EI26 

EI27 


Shuttle  Imaging  Radar 


EI31-EI41 

EI45 


Candidates : 


Large  Area  Modular  Array  of  Reflectors 

Gravity  Gradiometer 

Multi  spectral  Mid-IR  Imager 

Thermal  IR  Sensor 

Visible  IR  Sen., or 

Earth  Resources  Synthetic  Aperture  Radar  Unc 

National  Oceanic  Satellite  System  (NQSS) 

Upper  Atmospheric  Research  Satellite 
NOSS  Research  Program 

Ocean  Circulation  Mission  Topography  Experiment 

Operational  Meteorology  Satellite 

Large  Antenna  Mul ti frequency  Microwave  Radiometer 

Scatterometer 

Scanning  Multichannel  -Microwavt  Radiometer 
Instruments  on  E7 

Cryogenic  Limb  Scanning  Interferometer  Radiometer 
Salinity  Radiometer 
Ocean  Synthetic  Aperture  Radar 
Pushbroom  Radiometer 

Large  Area  Modular  Array  Reflectors 
Ocean  Research 


1981  Launch 

<1990  Start 

1982  Start 

1983  Start 
Mid  1980' s  Start 
Mid  1980's  Start 

Undetermined  Start  Date 


1981 

Start 

1982 

Start 

1982 

Start 

1983 

Start 

1980 ' s 

Start 

1981 

Start 

1981 

Start 

1981 

Start 

1982 

Start 

1984 

Start 

<1985 

Start 

<1985 

Start 

<1985 

Start 

1990 

Start 

1985 

Start 

Search  for  Extraterrestrial  Intelligence 
Pinhole  Satellite 

Tracking  and  Data  Relay  Satellite  II 
Unconventional  Orbital  Transfer  Vehicle 


Undetermined  Start  Date 
1984  Start 

1986  Start 

1990  Start 


*  Description  in  documents  too  vague  to  make  judgement. 
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Table  3-4.  Feasible,  NASA  [Invalidated  Missions 


Permanent  Orbiting  Outpost 
Geosynchronous  Telephone  Switching  Utility 
Geosynchronous  Broadband  Computer  Links 
Global  Search  and  Rescue  Location 
Orbiting  Civil  Command  and  Control  Station 
High  Resolution  Earth  Mapping  Radar 
10-m  Optical  Interferometer 
Advanced  Resources/Pollution  Observatory 
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All  of  the  others  were  judged  either  too  far  into  the  future,  too  speculative, 
or  possessed  less  than  average  desirability  with  respect  to  the  "envisioned" 
POC  test. 

3  2  AF  Missions 

Representative  DoD  missions  were  examined  next  to  determine  their 
requirements  and  also  their  suitability  as  test  structures.  It  was  found 
that  many  of  the  systems  are  classified  and  that  one  realistically  has  to 
settle  for  some  very  limited,  available  information  from  the  Air  Force  [2,  4]. 

table  3-5  shows  the  set  of  AF  missions  which  were  determined  to  be 
feasible  (in  the  same  sense  as  the  NASA  missions).  The  missions  are  listed 
according  to  underlying  functions:  radar,  communication,  passive  and  active 
optics.  The  time  frames  for  these  missions  are  unavailable  and  therefore 
could  not  be  assessed.  In  general,  all  of  these  systems  are  excellent 
test  beds.  The  radar  and  communication  systems  listed  are  LSS  having 
stringent  structural  control  requirements.  Moreover,  the  passive  and  active 
optical  missions  probably  represent  the  most  needful  of  the  technology 
that  is  proposed  to  be  demonstrated  by  this  study. 

The  document  "Advanced  Space  System  Concepts  and  Their  Orbital  Needs"  [4] 
proposes  an  extensive  set  of  potential  AF  missions.  Those  that  can  serve 
as  a  test  bed  are  listed  in  Table  3-6.  These  systems  must  be  viewed,  however, 
as  much  more  speculative  than  the  missions  listed  in  Table  3-5. 

3 . 3  Preliminary  Assessment  of  Existing  or  Planned  Systems 

The  systems  listed  in  Subsections  3.1  and  3.2  were  all  feasible.  The 
next  step  taken  in  the  task  was  to  determine  a  criterion  with  which  to 
select  the  baseline  POC  structure  and  choose  one. 

The  desirable  properties  for  a  POC  experiment  are  that  it  be  of  interest 
to  NASA  and  DoD,  yield  as  much  information  about  LSS  as  is  possible,  and 
be  low  cost  and  simple.  In  order  to  rank  the  identified  systems  according 
to  these  attributes,  each  system  was  to  be  graded  in  each  of  the  following 
nine  categories: 
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Table  3-5.  Feasible  AF  Mission  Models 


Radar: 

Mechanically  Steered  Space  Based  Radar  Surveillance  (U) 
Medium  Altitude  Parabolic  Cylinder  Radar  Surveillance  (U) 
Synchronous  Altitude  Active  Lens  Radar  Surveillance  (U) 
Medium  Altitude  Tactical  Surface  Reconnaissance  Radar  (U) 
Distributed  Array  Radiometer  Surveillance  (U) 

Medium  Altitude  Surveillance  Radar  (U) 

Communications : 

Maneuverable  Strategic  MILSATCOM  (U) 

Passive  Optics: 

Staring  IR  Missile  Tracking  (U) 

IR  Step-Stare  Mosaic  Surveillance  (U) 

IR  Panoramic  Pushbroom  Space  Surveillance  (U) 

Active  Optics: 

Space  Based  Laser  for  ASAT/DSAT  (U) 

Space  Based  Laser  Relay  (U) 

Space  Based  Laser  for  Boost  Phase  ABM  (U) 
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1)  Interest  to  OoD  and  NASA  objectives 

2)  Representative  of  LSS  -  many  closed  spaced,  low  frequency, 
low  damped  modes 

3)  Acceptable  dynamic  interaction  with  orbiter 

4)  Low  cost  and  simplicity 

5)  Ground  testability 

6)  Low  impact  on  intended  mission 

7)  Potential  for  future  growth 

8)  Low  weight  and  small  Shuttle  stowed  volume 

and 

9)  Easy  restowabil i ty 

The  total  score  was  then  to  be  collected  and  compared. 

It  became  evident  very  quickly,  however,  that  the  stated  approach  for 
choosing  a  structure  would  not  succeed.  The  problem  is  that  there  are 
really  no  candidates.  More  specifically,  hard  timeframes  and  commitments 
for  the  systems  identified  are  lacking.  A  multitude  of  competing  NASA 
missions  are  being  advanced,  yet  there  is  a  great  uncertainty  as  to  which 
ones,  if  any,  will  fly.  The  DoD  were  highly  suitable  candidates,  never¬ 
theless  the  timeframe  and  commitment  needed  to  judge  them  is  not  accessible. 
Another  challenging  problem  was  to  determine  if  the  sponsor  of  a  multi- 
million  dollar  spacecraft  could  he  convinced  that  a  set  of  experiments, 
deliberately  exciting  his  structure,  is  risk-free  and  to  his  benefit. 
Finally,  a  majority  of  the  missions  required  major,  expensive  modification 
and  added  instrumentation  to  permit  the  experiment.  Hence,  the  preliminary 
conclusion  was  that  a  dedicated  POC  structure  was  too  competitive  an 
alternative  to  ignore  and,  therefore,  had  to  be  considered  before  final 

trades  are  made. 

3.4  Dedicated  POC  Structure  Alternatives 

Dedicated  structures  were  considered  next  as  means  of  getting  more 
control  on  commitment  and  timeframe  and  also  to  improve  the  support  of 
the  sponsor.  Two  different  viewpoints  were  taken  in  generating  these 
dedicated  POC  structure  alternatives.  One  is  a  simple,  modular  approach 
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aimed  at  demonstrating  technology  through  emulations  of  a  real  system. 

The  other  is  geared  at  demonstrating  technology  on  pieces  of  a  potentially 

real  system. 


The  idea  of  the  modular  approach  is  to  split  the  structure  into  two 
pieces,  the  support  structure  and  the  payload.  One  then,  depending  on  the 
desired  technology  demonstration,  interchanges  various  support  structures 
with  the  alternative  payloads.  As  will  be  shown  later,  the  approach  offers 
simplicity  and  at  potentially  low  cost. 


The  idea  of  the  second  approach  is  to  use  pieces  of  existing  or  funded 
hardware  which  are  likely  to  be  used  in  future  systems.  The  advantage  of 
this  is  that  pieces  of  real  hardware  can  be  tested  and  later  be  integrated 
into  the  real  missions.  The  disadvantage  is  that  it  is  limited  to  the 
demonstration  of  system  unique  technologies  and  is  potentially  more  complex 
and  expensive  than  the  modular  approach. 


3.4.1  Modular,  Emulation  Concepts 

Two  support  structures  were  identified  for  use  in  the  modular  concept. 
One,  shown  in  Figure  3-2a,  is  an  Astromast  like  structure  that  has  one  end 
attached  to  a  shuttle  pallet  and  the  other  end  has  a  gimbal .  The  gimbal 
is  free  to  accept  one  of  several  payloads  to  be  identified  next.  The  mast 
itself  can  be  used  to  demonstrate  vibration  suppression  and  figure  control. 

The  second  support  structure,  shown  in  Figure  3-2b,  consists  of  two  con¬ 
centric  volumes.  The  inner  volume  contains  disturbance  generating  devices. 

The  outer  volume  is  attached  to  the  inner  volume  through  isolation  mechanisms, 
One  end  of  the  other  volume  can  be  attached  to  the  Shuttle  through  another 
isolation  mechanism  or  the  whole  structure  can  be  tethered.  At  the  other 
end  of  the  outer  volume  is  a  gimbal  free  to  accept  one  of  several  payloads. 
The  second  support  structure  can  be  used  to  demonstrate  vibration  isolation, 
suppression,  and  to  some  extent,  figure  control. 

One  method  conceived  for  adding  vibration  suppression  and  figure  control 
to  the  mast  support  structure  is  to  attach  collars.  Referring  to  Figure  3  3, 
force  actuators  connecting  the  collars  can  be  used  to  apply  a  controllable 
bending  moment.  Alternatively,  one  may  add  momentum  exchange  devices  that 
are  linear  or  rotational.  A  similar  approach  may  be  taken  with  the  con¬ 
centric  volume  structure. 
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Figure  3-2.  Modular,  emulation  Support  Structures 
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Figure  3-3.  Concept  for  Evaluating  Active 
Structural  Control 


Several  modular  payload  structures  were  identified  next  for  use  with 
either  of  the  described  support  structures.  These  are  shown  in  Figures  3-4 
and  3-5.  They  are  the  thin  rectangular  plate,  the  thin  planar  structure 
with  inexpensive  segmented  mirrors,  the  thin  rectangular  plate  with  tetra¬ 
hedron,  and  the  thin  twin  beam.  They  all  attach  directly  to  the  support 

structure  gimbal . 

The  three  structures,  shown  in  Figure  3-4,  can  be  used  to  emulate  the 
dynamics  of  antennas,  segmented  optics,  and,  with  some  imagination,  antennas 
or  mirrors  with  reflectors.  Bending  moment  and  momentum  exchange  actuation 
attached  on  the  planar  part  of  the  structure  can  be  used  to  demonstrate 
vibration  suppression,  figure  control,  and,  when  the  gimbal  is  rotated, 
post  slew  damping.  The  addition  of  the  tetrahedron  in  Figure  3-4b  adds 
another  level  of  emulation  and  can  be  used  to  demonstrate  isolation  and 
the  control  of  weakly  controllable  modes.  The  tetrahedron  can  also  be 
added  to  the  structure  with  segmented  mirrors.  The  picture  of  the  Sun¬ 
flower  antenna  in  Figure  3-4d  is  shown  to  demonstrate  what  is  meant  by 
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Figure  3-4.  Modular,  Emulation  Payload  Structures 


emulation.  The  topological  resemblance  between  the  plate  with  tetrahedron 
and  the  Sunflower  is  self-evident. 

The  thin  beam  payload  shown  in  Figure  3-5  also  attaches  to  the  support 
structure  gimbal.  It  represents  a  concept  capable  of  demonstrating  optimum 
slew  of  LSS  with  large  solar  panels  but  it  can  be  configured  to  demonstrate 
vibration  suppression,  figure  control,  rotational  vibration  isolation  and 
multibody  alignment.  The  structure  consists  of  two  parallel,  very  thin 
beams.  The  beams  are  hinged  (in-plane)  to  each  other  to  reduce  the  effects 
of  torsional  modes.  The  two  beams  attach  to  the  support  structure  gimbal 
through  an  interface  consisting  of  two  concentric  motors.  The  inner  motor 
is  used  to  disturb  the  structure,  the  outer  motor  is  used  to  rotationally 
isolate  the  disturbance  and  also  effect  the  slew.  Bending  moment  and/or 
momentum  exchange  devices  are  used  to  suppress  "leaked"  vibration  and  pre¬ 
serve  figure.  The  LED  and  receiver  are  used  to  monitor  the  test.  The 
multibody  alignment  experiment  adds  additional  LED  transmitters,  light 
receiver  and  a  motorized  gimbal.  The  motor  is  used  to  effect  alignment  of 
one  end  of  the  beam  to  body  center.  The  motorized  gimbal  would  be  attached 
at  one  end  of  the  beam  and  the  rotational  momentum  exchange  would  be  con¬ 
trolled  to  perform  the  alignment. 

Figure  3-6  illustrates  how  the  shuttle  attached  and  tethered  configura¬ 
tion  would  look  using  the  modular  approach.  The  motion  of  the  structure 
would  be  monitored  by  attaching  LED  transmitters  to  selected  points  and  tian 
observing  the  motion  from  a  shuttle  based  receiver.  Pointing  experiments 
would  reverse  this:  a  receiver  could  be  added  to  the  structure  and  an  LED 
transmitter  would  be  placed  in  the  shuttle  cargo  bay. 

3.4.2  Semi -Real  Concepts 

Alternative  concepts  incorporating  pieces  of  planned,  funded,  or  existing 
hardware  were  also  considered.  With  few  exceptions,  the  majority  are  large 
antennas:  e.g.,  the  General  Dynamics  "Precision  Erectable  Truss  Antenna" 

PETA),  the  Lockheed  "Wrap-Rib  Antenna",  the  Harris  "Hoop-and-Col umn" ,  and 
the  TRW  "Sunflower".  These  systems  were  invariably  rejected  because  of 
their  cost  and  complexity  (for  an  experiment)  and  because  it  is  not  clear 
if  active  control  is  even  required  in  these  systems.  For  this  reason,  the 
decision  to  depart  from  antennas  and  into  optical  systems  was  made. 
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Figure  3-6.  Illustration  of  Attached  and  Tethered  Configuration 
Using  the  Modular  Approach 


Two  optical  structures  were  conceptualized.  The  first,  sketched  in 
Figure  3-7,  is  of  an  Offset  Cassegrain  structure.  The  primary  is  parabolic, 
the  secondary  is  hyperbolic  and  the  focal  plane/detector  is  sideways  behind 
the  primary  mirror.  The  entire  structure  fits  in  one  shuttle  pallet.  The 
concept  was  equipped  with  actuators,  sensors,  and  an  optical  train  to  enable 
one  to  test  vibration  control,  isolation,  figure  control  and,  via  monitored 
shuttle  motion,  active  damping.  The  optical  f iel d-of-view  of  this  was  ex¬ 
pected  to  be  small,  perhaps  0.5° x 0.5°.  The  main  experiment  therefore 
envisioned  to  image  a  star  as  a  point  light  source  at  infinity.  The  quality 
of  the  image  in  the  presence  of  various  disturbance  acting  on  the  optical 
train  would  then  be  used  to  assess  the  effectiveness  of  active  structural 
control .  Earth  imaging  was  a  secondary  experiment. 

The  second  optical  structure  conceptualized  is  sketched  in  Figure  3-8. 
This  configuration  retains  the  essence  of  the  two  mirror  system  but  eliminates 
the  secondary  mirror  and  relocates  the  focal  plane/detector  to  where  the 
secondary  used  to  be.  It  was  initially  believed  that  this  approach  would 
result  in  a  more  reasonable  optical  geometry  while  retaining  the  geometric/ 
jitter  sensitivities  projected  for  future  systems.  It  will  be  shown  in 
Section  4.0  that  this  is  not  true. 

3 . 5  Test  Structure  Candidates 

The  addition  of  the  dedicated  systems  to  the  NASA  and  AF  missions  made 
the  trades  described  in  Section  3.3  possible.  The  results  of  these  trades 
are:  no  near  term  NASA  structures  suitable  to  meet  DARPA  POC  objectives 
could  be  identified.  Those  that  are  of  interest  are  in  the  future  and  are 
too  uncertain  to  commit  a  POC  experiment  design.  Many  of  the  DoD  missions 
examined  could  be  used  but  the  data  required  to  decide  on  whether  the 
structure  is  really  usable  is  unavailable.  Hence,  reliance  on  projected 
missions  to  design  a  POC  test  is  too  risky.  Strictly  generic  structures 
were  also  considered  but  these  are  hard  to  justify  for  a  costly  flight  ex¬ 
periment.  The  POC  test  has  to  simultaneously  demonstrate  a  technology  to 
make  it  worthwhile.  Existing  and  funded  systems  were  additionally  considered 
to  meet  this  later  requirement  but  these  systems  were  found  to  be  antennas; 
too  complex,  too  expensive,  and  it  was  not  even  clear  that  they  required  the 
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Figure  3-7.  Offset  Cassegrain,  Two  Mirror  Telescope  Concept 
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Figure  3-8.  Single  Mirror  Telescope  Concept 


technology.  Therefore,  since  DARPA's  interests  seem  to  lie  in  large  optical 
systems,  we  recommended  the  study  of  an  on-orbit  POC  structure  supporting 
this  technology.  The  two  telescope  structures  incorporating  a  lightweight 
mirror  were  therefore  chosen  for  '  ^ther,  more  detailed  analysis. 


3-23 


4.0  TEST  STRUCTURE  DESIGN  AND  ANALYSIS 


Two  concepts  for  the  optical  system  configuration  were  considered: 
the  single  mirror  concept  and  the  two  mirror  Cassegrain  configuration. 

In  the  following  the  two  concepts  are  discussed  and  the  advantages  and 
drawbacks  of  each  approach  are  explained.  The  reasons  for  the  selec¬ 
tion  of  the  two  mirror  approach  and  the  details  of  the  chosen  configura¬ 
tion  are  given.  The  model  data  obtained  from  finite  element  program  is 

pres  ented . 


4.1  Optical  Prescription  for  Single  Mirror  Concept 

In  designing  the  optical  system  we  are  trying  to  achieve  two  goals: 

(i)  To  arrive  at  an  optical  system  exhibiting  similar  sensitivity  to 
geometric  distortions  and  jitter  anticipated  in  optical  trains  of  future 

systems. 

(ii)  To  make  use  of  existing  optical  hardware,  with  possibly  some 
modifications. 


The  single  mirror  configuration  was  considered  as  a  possible 
candidate  for  the  optical  system  because  its  simple  geometry  results 
in  a  less  complicated  structure.  The  simplicity  of  the  single  mirror 
design  is  apparent  in  the  ray  traces  of  its  optical  prescription  as 
shown  in  Figure  4-1.  The  effectiveness  of  the  active  structural  control 
system  will  be  assessed  by  imaging  a  star  as  a  point  light  source  at 
infinity.  The  quality  of  the  image  in  the  presence  of  various  distur¬ 
bances  acting  on  the  optical  train  will  serve  as  the  performance  measure 

of  the  control  system. 


Although  the  single  mirror  configuration  has  a  simple  structure, 
it  suffers  from  an  inherent  nonuniformity  of  the  image  resolution  over 
the  field  of  view  (field  angle)  as  can  be  seen  in  Figure  4-2.  This 

inherent  drawback  of  the  single  mirror  optical  system  makes  it  un¬ 

acceptable  as  a  POC  test  structure:  The  vernier  thruster  (VRCS)  shuttle 
limit  cycle  has  excursions  of  +0.1°.  Thus  at  least  a  0.3°x0.3  field 
of  view  will  be  required  to  image  a  star.  The  motion  of  the  star 

image  in  the  focal  plane  due  to  the  shuttle  limit  cycle  occurs  at  a 

very  slow  (0.01  deg/sec)  and  predictable  rate  and  can  be  easily  removed 
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post  facto  from  the  observation  data.  However,  it  is  extremely  difficult, 
if  not  impossible  to  decide  how  much  of  the  loss  of  resolution  is  caused 
by  vibrations  in  the  optical  structures  which  are  to  be  controlled  or 
by  the  resolution  deterioration  incurred  as  the  star  image  moves  away 
from  the  center  of  the  field  of  view  (FOV). 

4.2  The  Two  Mirror,  Offset  Cassegrain  Structure 

Because  of  the  drawbacks  of  the  single  mirror  concept  explained 
in  Section  4.1,  it  was  decided  to  go  with  a  two  mirror,  offset  Cassegrain 
system.  This  configuration  allows  for  an  optical  design  which  provides 
a  more  uniform  resolution  over  the  FOV  as  compared  to  the  single  mirror 
approach.  The  structure  for  the  Cassegrain  system  is  designed  to 
have  the  focal  plane  assembly  in  between  the  primary  and  the  secondary 
mirrors.  The  total  distance  along  the  optical  center  between  the  primary 
and  secondary  surfaces  is  only  1.77  m,  assuming  a  mast  not  exceeding 
2.5  meters.  Thus,  the  entire  structure  fits  within  the  Shuttle  cargo 
bay  envelope  with  cargo  bay  doors  closed  and  therefore  no  deployment 
mechanisms  are  needed. 

The  optical  prescription  of  the  two  mirror,  offset  Cassegrain 
structure  is  shown  in  Figure  4-3.  It  represents  a  preliminary  optical 
design,  which  is  more  than  adequate  for  our  present  purpose  of 
defining  an  opto-mechanical  test  structure.  The  optical  configuration 
of  the  Cassegrain  structure  consists  of  two  mirrors:  a  parabolic 
primary  mirror  and  a  hyperbolic  secondary.  The  salient  data  of  the 
system  is  given  in  the  following: 


System 

Effective  Focal  Length:  3.5m 
Fiel d-of-View:  0.01x0.01  rad 
Focal  Plane:  3. 5x3. 5  cm 


Primary:  Parabolic  Surface 
Diameter  of  Off-Axis  Segment:  1.0m 
Radius  of  Curvature:  -4.1464m 
Focal  Length:  2.0732m 


Secondary:  Hyperbolic  Surface 
Aperture  Diameter:  0.15m 
Radius  of  Curvature:  -1.4718m 
Focal  Length:  0.4032m 
Eccentricity:  3.906 
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Figure  4-3.  Optical  Prescription 


The  FOV  in  this  configuration  is  0.01x0.01  rad  (0.57x0.57  deg). 

It  allows  imaging  a  star  while  the  Shuttle  undergoes  +.1  degree  limit 
cycle  motions.  As  mentioned  before,  this  highly  predictable  and  very 
slow  Shuttle  rate  of  0.01  deg/sec  can  be  removed  from  the  data  post 
facto.  Another  advantage  is  that  the  reasonable  FOV  size  also  makes 
earth  viewing  a  meaningful  experiment. 

The  performance  of  this  preliminary  optical  design  is  shown  in 
Figure  4-4.  By  shifting  the  focal  plane  by  +1.8mm  out  of  focus  a 
relatively  uniform  resolution  of  100  yrad  over  the  FOV  can  be  achieved. 
This  resolution  is  not  sufficient  for  an  actual  on-orbit  test  structure 
and  a  fine  tuned  optical  design.  Improvement  of  resolution  and  the 
uniformity  of  resolution  are  needed  before  such  a  system  could  be 
built.  This,  however,  is  not  the  intent  at  this  point.  At  this  stage 
it  is  merely  important  that  the  system  concept  selected  and  the  result¬ 
ing  structure  allow  for  a  high  quality  optical  design.  This  can  be 
accomplished  by  either  figuring  the  primary  and  secondary  surfaces  as 
aspherics,  or  by  adding  a  third  reflective  element.  A  very  desirable 
approach  for  an  on-orbit  P0C  structure  would  be  to  leave  the  large 
primary  surface  as  a  parabola,  add  a  refractive  element  between  the 
secondary  and  the  focal  plane,  aspherize  the  secondary  surface  and 
design  the  refractive  element  to  optimize  the  optical  performance. 

Of  key  interest  here  is  the  fact  that  such  optical  design  improvements 
will  not  affect  the  structural  concept  nor  its  physical  dimensions. 

4.3  Baseline  Structure  Design  of  the  Two-Mirror  Offset  Cassegrain  Concept 

The  baseline  test  structure  was  designed  to  withstand  launch  and 
landing  loads  while  being  extremely  complient  in  space.  To  accomplish 
this  task,  flexible  structures  are  used  to  support  the  mirrors  and  the 
sensor  package.  These  structures  are  in  turn  supported  by  caging 
devices  during  period  of  critical  loading. 

The  design  of  the  test  structure  is  shown  in  Figure  4-5.  It 
consists  of  an  L  shaped  structure  which  supports  the  primary  and 
secondary  mirrors,  a  mast  and  primary  mirror  caging  mechanisms, 
primary  mirror's  coarse  and  fine  actuators,  and  secondary  mirror  control 
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Figure  4-4.  Performance  of  the  Offset  Cassegrain  Concept 
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actuators.  The  structure  is  mounted  on  a  base  which  is  connected  to 
the  Shuttle  pallet  at  ten  attachment  points. 

The  coarse  actuators  of  the  primary  mirror  system  are  mounted  on 
a  graphite-epoxy  honeycomb  panel.  This  panel  is  attached  to  the  base 
by  two  paths.  When  the  experiment  is  not  being  used,  a  set  of  caging 
devices  is  closed,  rigidly  connecting  the  mirror  mount  to  the  base. 

When  testing,  these  caging  devices  are  opened  to  allow  limited  free 
movement  of  the  mirror  (+3mm).  However,  movement  large  enough  to 
damage  the  structure  is  prevented  by  a  pair  of  low  thermal  expansion 
graphite  epoxy  flexural  hinges  that  support  the  primary  mirror  system. 
The  fundamental  frequency  of  the  hinge  mirror  system  is  ~2  Hz. 

The  primary  mirror  caging  mechanism  is  shown  in  Figure  4-6a.  It 
consists  of  a  motor  and  gearhead  combination  which  automates  the 
opening  and  closing  of  jaws  through  a  rack  and  pinion  gear  train.  The 
jaws  close  on  a  T-bar  which  is  connected  to  the  mirror  support.  The 
major  forces  exerted  by  tne  T-bar  on  the  jaws  are  perpendicular  to  the 
allowed  movement  of  the  jaws,  so  the  motor  does  not  carry  the  load 
where  the  jaws  are  open.  The  primary  mirror  is  allowed  limited  free 
motion  in  all  axes. 

The  support  mast  for  the  secondary  mirror  and  sensor  package  has 
a  fundamental  resonance  frequency  of  2.7  Hz.  The  mast  is  caged  and  is 
also  made  of  low  thermal  expansion  graphite-epoxy.  It  is  a  tube 
approximately  6  cm  in  diameter  by  2mm  thick.  The  secondary  mirror  is 
attached  to  the  top  of  the  mast  and  has  variable  height  and  angle 
actuators.  The  sensor  package  is  one  third  of  the  way  down,  see 
Figure  4-5. 

The  caging  devices  for  the  mast  are  shown  in  Figure  4-6b.  These 
devices  have  two  movable,  beveled  jaws  on  each  side  of  the  mast.  The 
jaws  fit  into  matching  "V"  slots  in  a  bar  through  the  mast.  This 
provides  caging  in  all  axes  when  the  jaws  are  closed.  The  caging 
devices  are  attached  to  the  pallet  by  a  truss  of  tubes  from  the  pallet. 
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MIRROR  SUPPORT 


(a)  Primary  Mirror  Caging  (b)  Mast  Caging 


4.4  facial  Analysis  of  the  Test  Structure 

A  finite  element  model  of  the  structure  with  the  two  mirrors 
attached  to  it  was  developed.  The  structural  model  included  43  truss 
elements  and  31  two  dimensional  elements.  The  resultant  dynamic  model 
included  50  nodes  each  one  having  six  degrees  of  freedom  in  space, 
giving  a  total  of  300  degrees  of  freedom.  In  developing  the  test 
structure's  model,  it  was  assumed  that  both  mirrors  were  very  stiff, 
the  fundamental  frequency  of  the  primary  mirror  was  assumed  to  be  '400  Hz. 
The  modal  analysis  was  performed  using  the  TRW  finite  element  program 
TRWSAP. 

Two  sets  of  modal  frequencies  and  mode  shapes  were  generated: 

(i)  The  free- free  modal  representation  of  the  test  structure 
including  mirrors. 

(ii)  The  modal  representation  of  the  test  structure  (mirror 
included)  assuming  that  the  base  is  attached  to  the  ground.  These 
modes  will  be  referred  to  as  the  cantilevered  modes. 

Tables  4-1  and  4-2  show  some  of  the  modal  frequencies  for  both  the 
free-free  and  cantilevered  boundary  conditions.  The  first  four  mode 
shapes  of  each  case  are  shown  in  Figures  4-7  and  4-8.  As  it  can  be 
seen  from  the  tables,  there  is  a  significant  jump  in  free- free  model 
frequencies  after  the  first  two  modes.  However,  when  the  structure 
is  cantilevered  two  more  modal  frequencies  are  lowered.  In  general, 
it  can  be  said  that  the  structure  behaves  like  two  cantilevered  beams. 


4-11 


Table  4-1.  Free-free  Modal  Frequencies  of  Test  Structure 


MODE 

NUMBER 

FREQUENCY 

W 

DESCRIPTION 

1-6 

0 

Rigid  body  modes 

7 

2  896 

Secondary  mirror  moves  along  x,  primary 
moves  along  z 

8 

3.767 

Primary  mirror  rotates  around  x. 

Secondary  mirror  moves  along  y 

9 

17.860 

Secondary  moves  along  x,  primary  rotates 
around  y 

10 

20  41  3 

Secondary  mirror  moves  in  y  direction 

11 

40.644 

Mast  motion  in  yz  plane 

12 

45.148 

Mast  bending  in  xz  plane,  primary  mirror  motion 
in  z  direction 

13 

127.831 

Rotation  of  secondary  mirror  about  the  y-axis 

14 

147.467 

Rotation  of  secondary  mirror  around  x-axis 

15 

157.887 

Rotation  of  secondary  mirror  around  x-axis  and 
motions  of  primary  mirrors  support 

Table  4-2.  Cantilever  Modal  Frequencies  of  Test  Structure 


MODE 

NUMBER 

FREQUENCY 

fs(Hz) 

DESCRIPTION 

1 

1  .970 

Bending  of  Primary  mirror  in  xz  plane 

2 

2.637 

Bending  of  mast  (and  secondary  mirror)  in 
zy  plane 

3 

2.640 

Bending  of  mast  (and  secondary  mirror)  in 
xz  plane 

4 

4.872 

Torsion  of  primary  mirror  around  the  y-axis 

5 

22.360 

Bending  of  mast  (and  secondary  mirror)  in 
xz  plane 

6 

23.014 

Bending  of  mast  in  yz  plane 

7 

40.881 

Bending  of  mast  in  xz  plane 

8 

48.444 

Torsion  of  primary  mirror  around  y-axis 

9 

50.654 

Motion  of  primary  mirror  in  z  direction 

10 

124.505 

Motion  of  secondary  mirror  in  t  direction 
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5.0  DYNAMIC  INTERACTION  WITH  THE  ORBITER 

The  motivation  for  the  dynamic  interaction  study  is  to  assess  the 
effect  of  mounting  the  structure  on  the  normal  operation  of  the  Shuttle 
Orbiter.  In  addition,  one  would  like  to  determine  the  influence  of  the 
Shuttle  thrusters  firings  on  the  structure's  oscillation  so  that  they  car 
be  taken  into  account  in  the  design  of  the  control  law  demonstration  ex¬ 
periments. 

Areas  of  special  interest  are  the  effect  of  added  structural  flexibility 
on  the  system's  dynamic  behavior  and  the  influence  of  excitation  forces 
applied  to  the  structure  during  experiments  on  the  normal  operation  of  the 
Digital  Autopilot  (DAP).  It  will  be  shown  that  the  additional  vibrational 
modes  of  the  test  structure  mounted  on  the  Orbiter  do  not  couple  with  the 
Shuttle's  modes  and  therefore  do  not  result  in  high  amplitude  vibrations. 

It  will  be  also  shown  that  the  excitation  forces  which  act  on  the  structure 
during  experiments  do  not  affect  the  Shuttle's  RCS  system  operation.  The 
simulation  study  shows  that  the  number  of  thrusters  firings  does  not  change 
significantly  with  the  addition  of  the  test  structure  and  thus  the  fuel 
consumption  remains  virtually  unchanged. 

Another  area  of  interest  investigated  was  the  amount  of  excitation 
exerted  on  the  structure  by  the  RCS  firings.  The  level  of  this  excitation 
must  be  assessed  in  order  to  correctly  design  the  experiments  and  to  deter¬ 
mine  the  additional  excitations  needed  to  demonstrate  the  control  law  per¬ 
formance. 

A  dynamic  simulation  study  was  performed  using  a  TRW  developed  digital 
simulation  program.  The  dynamic  model  of  the  coupled  system  was  obtained 
by  combining  the  modal  representations  of  the  Orbiter  and  the  test  structure 
using  a  TRW  modal  synthesis  program,  KQMBINE.  A  detailed  model  of  the  DAP 
was  implemented  using  the  up-to-date  parameters  obtained  from  JSC.  Special 
features  were  added  to  the  program  to  simulate  the  conditions  during  the 
experiments.  The  simulation  program  is  highly  modular  and  can  be  easily 
modified  by  adding  and/or  changing  specific  subroutines  which  might  be 
needed.  This  feature  allows  one  to  perform  simulation  studies  on  difierent 
structures  and  various  flight  conditions  with  minimum  additional  irogramming 

effort. 
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5.1  Systems 's  Dynamic  Representation 

The  system's  dynamic  model  was  obtai  led  by  combining  the  modal  data 
of  the  Orbiter  and  the  test  structure.  The  modal  synthesis  approach  allows 
one  to  make  changes  in  the  test  structure  without  running  an  expensive 
finite  element  program  each  time  a  dynamic  model  of  the  total  system  is 
required.  The  combination  procedure  was  performed  by  a  TRW  developed 
program  KOMBINE  which  will  be  described  in  the  following  section.  Since 
the  dynamic  model  of  the  Shuttle  obtained  from  Rockwell  International  was 
given  in  terms  of  the  Orbiter's  structural  coordinate  system  [8],  this  coordi¬ 
nate  system  was  adopted  as  the  reference  system  for  the  current  study  (see 
Figure  5-1).  In  order  to  avoid  confusion,  all  quantities  and  locations 
were  transformed  into  this  coordinate  system.  This  includes  DAP  parameters 
and  thruster  firing  directions  expressed  in  different  coordinate  systems  in 
Shuttle's  original  documents. 

5.1.1  Modal  Synthesis  Program  KOMBINE 

Most  of  the  existing  structure  modeling  computer  programs  are  general 
purpose  programs  and  therefore  are  very  expensive  to  run.  The  modal  synthesis 
program  developed  in  TRW  permits  fast  and  efficient  modeling  of  a  flexible 
assembly  which  consists  of  any  number  of  interconnected  flexible  and  rigid 
bodies.  The  program  reads  the  free-free  modal  representation  of  each  body 
in  the  assembly  and  combines  them  sequentially  into  one  system.  The  pro¬ 
gram  can  accommodate  very  general  interfaces  between  any  two  adjacent  bodies, 
including  gimballed  and  flexible  connections.  Moreover,  any  two  bodies  can 
be  connected  at  more  than  one  point  which  allows  for  closed  topological 
loops.  The  resultant  dynamic  model  is  linear  since  the  bodies  modal 
representation  is  used  as  input.  However,  the  program  allows  large  angle 
relative  rotations  between  any  two  adjacent  bodies.  Therefore,  it  is 
possible  to  conduct  sensitivity  studies  of  the  dynamic  representation  as 
a  function  of  system's  configuration.  The  program  has  also  an  option  of 
calculating  the  mass  properties  of  the  combined  body  and  of  any  of  the  given 
substructures.  This  feature  is  useful  in  propulsion  system  and  fuel  con¬ 
sumption  studies.  The  program  accepts  modal  data  calculated  by  standard 
finite  element  programs  such  as  NASTRAN  and  TRWSAP.  Moreover,  KOMBINE  is 


5-2 


ORIGIN:  In  the  orbiter  plane  of  symmetry,  400  inches 

below  the  center  line  of  the  payload  bay  and 
235  inches  forward  of  the  tip  of  the  orbiter 
nose  fairing. 

ORIENTATION:  The  XQ  axis  is  in  the  vehicle  plane  of  symmetry, 

parallel  to  and  400  inches  below  the  payload 
bay  centerline.  Positive  sense  is  from  the 
nose  of  the  vehicle  toward  the  sail. 

The  Z  axis  is  in  the  vehicle  plane  of  symmetry, 
perpendicular  to  the  XQ  axis.  Positive  upward 
in  landing  attitude. 

The  Y  axis  completes  a  right-handed  system, 
o 

Figure  5-1.  The  System's  Reference  Coordinate  System 
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compatible  with  existing  control  system  analysis  and  design  programs  such 
as  FLYNSYS.  The  flow  chart  diagram  of  the  program  is  shown  in  Figure  5-2. 

KOMBINE  is  as  accurate  as  a  finite  element  program  but  its  run  costs 
are  1-2  orders  of  magnitude  less.  This  is  achieved  by  using  numerically 
stable  and  efficient  algorithms.  The  procedure  avoids  inversion  of  ill- 
conditioned  matrices  by  eliminating  the  constrained  modes  using  singular 
value  decomposition  which  does  not  require  matrix  inversion.  The  program 
can  be  run  in  either  batch  or  timeshare  modes,  although  for  most  applications 
the  memory  available  in  the  timeshare  mode  is  sufficient.  The  average  cost 
for  combining  2  to  4  bodies  is  $l-$5. 

5.1.2  Combined  Systems  Modal  Representation 

The  dynamic  representation  of  the  orbiter  was  obtained  from  Rockwell 
International  in  the  form  of  free-free  modal  frequencies  and  displacements. 

The  given  model  included  200  modes  and  567  modal  displacements.  The  structure 
was  modeled  as  being  attached  to  one  of  the  existing  cargo  ports.  The  cargo 
port  chosen  is  one  of  the  right-hand  side  attach  points  at  XQ  =  1017.870, 

Y  =  94,  Z  =  424  (node  number  369  in  the  Orbiter1 s  model).  The  attachment 

0  0  .  j 

to  this  node  is  representative  and  is  used  to  generate  a  typical  system  s 
behavior.  Attachment  to  any  other  point  would  not  alter  the  conclusions 
of  the  study. 

The  model  of  the  complete  system,  Orbiter  and  test  structure  was  obtained 
by  combining  30  free-free  modes  of  the  Orbiter  with  30  free-free  modes  of  the 
structure  by  using  the  modal  synthesis  program  KOMBINE.  This  particular 
number  of  modes  was  chosen  after  concluding  that  a  higher  number  of  modes 
does  not  change  the  dynamic  representation  of  the  complete  system.  In  order 
to  keep  the  combined  system's  model  at  a  manageable  size  for  simulation  and 
analysis  studies,  only  selected  nodes  were  retained  in  the  final  model. 

Table  5-1  shows  the  numbering  of  the  selected  nodes  and  their  locations. 

A  schematic  drawing  of  the  structure  indicating  nodes'  locations  is  given 
in  Figure  5-3.  The  CPU  run  time  of  the  modal  synthesis  process  was  ^  2  sec 
and  the  total  run  cost  on  the  CDC  CYBER175  was  $3.00.  Thus,  model  generation 
is  easy  and  inexpensive  and,  if  future  changes  are  needed  either  in  the 
structural  model  or  in  the  number  of  modes  used,  they  can  be  easily  imple¬ 
mented  in  the  simulation. 
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HEAD  NODAL  DATA 
OF  A.  CALCULATE 

NASS  MOWTTES 
CALI  NASSP 

HEAD  MODAL  BATA 
|OF  A.  CALCULATE 
AI,  II.  Cl 
(CALL  SELECT) 


Figure  5-2.  Flowchart  of  Program 
KOMBINE 


Table  5-1.  System's  Node  Location 
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Table  5-2  summarizes  the  modal  representation  of  the  combined  structure. 
In  addition  to  the  system's  modal  frequencies,  the  free-free  natural  fre¬ 
quencies  of  the  Orbiter  and  of  the  structure  are  given.  The  modal  frequencies 
of  the  structure  rigidly  attached  to  the  ground  are  given  for  comparison. 

It  can  be  se  n  that  the  modes  of  the  combined  system  are  very  close  to  the 
modes  of  the  cantilevered  structure.  This  is  to  be  expected  since  the  mass 
of  the  Shuttle-Orbiter  is  very  large  compared  to  the  structure  and,  hence, 
for  all  practical  purposes  the  structure  behaves  as  if  it  were  cantilevered. 

It  can  also  be  observed  that  the  Shuttle's  modal  frequencies  are  not  affected 
by  the  addition  of  the  test  structure.  This  is  due  to  the  large  mass  of 
the  Orbiter  and  also  because  of  the  relatively  low  flexibility  at  the  attach¬ 
ment  point. 

The  vibrational  modes  of  the  structure  as  attached  to  the  Orbiter  can 
be  divided  into  three  distinct  groups: 

(i)  Modes  with  modal  frequencies  1.970,  2.639,  2.645,  4.879  Hz 

(ii)  Modes  with  modal  frequencies  22.732,  23.040  Hz 

(iii)  Modes  with  modal  frequencies  40.884,  51.33  Hz. 

In  designing  the  vibration  suppression  control  system,  this  natural  division 
of  vibrational  modes  is  taken  into  account.  The  control  system  will  be  de¬ 
signed  to  control  the  first  four  vibrational  modes.  As  it  will  be  shown 
later,  two  additional  modes  have  to  be  taken  into  account  during  the  identi¬ 
fication  study. 

5.2  Simulation  Program  Description 

5.2.1  Simulation  Program  Functional  Description 

A  detailed  simulation  program  based  on  the  TRW's  general  simulation 
program  XSIM  was  developed.  The  program  was  written  in  a  modular  form 
which  facilitates  additions  of  new  options  into  the  program  without  alter¬ 
ing  its  general  structure.  Moreover,  if  changes  have  to  be  made,  only 
the  modified  part  of  the  program  is  affected.  The  program  was  divided 
into  seven  autonomous  sections,  each  one  of  them  implemented  by  a  separate 
subroutine.  Figure  5-4  shows  the  functional  description  of  the  simulation 
program.  The  function  of  each  part  of  the  program  is  written  inside  each 
block  and  the  name  of  the  subroutine  implementing  it  is  denoted  above  the 
corresponding  block. 
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Table  5-2.  Natural  Frequencies 


f 


SYNTVESIS 


The  different  functional  parts  of  the  program  and  the  subroutines 
which  implement  each  part  are  as  follows:  (I)  The  plant  (DYNAMIC)  -  The 
state  space  model  of  the  Orbiter  and  the  test  structure  as  obtained  from 
the  model  synthesis  program  KOMBINE.  The  inputs  to  DYNAMIC  are  the  thrust 
forces  of  the  VRCS  system  and  disturbance  forces  acting  on  the  structure. 

The  outputs  are  the  three  attitude  angles  at  the  I MU  sensor  location  and 
the  deflections  of  the  nodal  points  on  the  test  structure.  (II)  Phase 
plane  logic  (PPLANE)  -  calculates  the  required  control  torques  needed  to 
maintain  a  given  attitude.  Subroutine  PPLANE  accepts  the  estimated  attitude 
error,  (0  ),  attitude  rate  error  (0g),  and  acceleration  caused  by  disturbances 
(Ur))  and  computes  the  control  torque  (C)  needed  for  attitude  correction. 

The  objective  is  to  limit  the  attitude  error  to  less  than  0.1  degree  and 
attitude  rate  error  to  less  than  0.01  deg/sec.  (Ill)  Jet  selection  logic 
(JETSEL)  selects  the  vernier  jets  that  must  fire  to  maintain  the  system's 
attitude.  JETSEL  obtains  the  control  torques  command  C  and  calculates 
the  jet  numbers  J  which  a^e  then  supplied  as  inputs  to  subroutine  INSUB. 

(IV)  Conversion  of  jet  numbers  to  corresponding  forces  (INSUB)  -  This  sub¬ 
routine  converts  the  numbers  of  the  jets  that  fire  at  any  given  time  into 
corresponding  force  components  acting  on  the  plant.  (V)  State  estimator 
model  (FILTER)  -  Two  filters  which  estimate  and  interpolate  the  Orbiter's 
attitude  rate  and  disturbance  acceleration.  Subroutine  FILTER  accepts 
the  attitude  parameters  of  the  plant  and  performs  the  following  operations: 

(a)  Estimates  the  attitude  parameters  and  disturbance  acceleration  given 
their  measurements  (sample  at  6.25  Hz)  implemented  by  subroutine  FILTER1 ; 

(b)  Extrapolates  attitude  parameters  over  half  the  sampling  period  of 
FILTER1,  implemented  by  subroutine  FILTER2.  (VI)  Disturbance  generator 
(DISTURB)  -  Simulates  the  disturbance  forces  acting  on  the  structure  during 
experiments.  It  models  the  disturbance  forces  (D)  acting  on  the  structure 
during  control  law  performance  demonstration  (band  limited  gaussian  noise) 
and  during  the  parameter  identification  experiment  (pulse  sequence). 

(VII)  Line-of-sight  calculations  (MIRSENS)  -  calculates  LOS  motion,  given 
the  motion  of  various  mirrors'  support  points. 

Subroutines  PPLANE,  JETSEL  and  FILTER  constitute  the  Digital  Autopilot 
(DAP)  of  the  Orbiter.  The  various  parts  of  the  DAP  operate  at  a  sampling 
rate  of  12.5  Hz,  except  the  estimation  part  of  subroutine  FILTER  (FILTER1) 
which  operates  at  6.25  Hz.  The  parameters  used  in  the  implementation  of 
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the  DAP  were  obtained  recently  from  JSFC  and  are  up  to  date.  The  detailed 
description  of  the  DAP  is  given  in  various  Shuttle  documents  [9].  A  brief 
description  of  the  DAP  is  given  in  Appendix  C.  In  the  following,  the  sub¬ 
routines  specific  to  the  current  simulation  study  are  described. 

5.2.2  Disturbance  Generating  Subroutine  DISTURB 


DISTURB  is  a  subroutine  which  generates  the  required  disturbance  forces 
as  a  function  of  time  and  applies  them  to  various  points  on  the  structure. 

In  the  current  study  two  types  of  functions  were  simulated: 

(i)  A  band  limited  Gaussian  process  which  simulates  artificial 

disturbance  forces  acting  on  the  structure.  These  disturbance 
forces  (Shakers)  induce  the  vibration  used  to  assess  the  per¬ 
formance  of  the  proposed  control  law.  The  shakers  themselves 
simulate  expected  environmental  disturbances,  e.g.,  pumps, 
CMG's,  etc.,  which  act  on  a  spacecraft  during  normal  operation. 
The  subroutine  which  generates  the  random  process  permits 
specification  of  the  standard  deviation,  the  cutoff  frequency, 
and  the  mean  of  the  process.  An  important  feature  of  the 
subroutine  is  the  fact  that  the  noise  bandwidth  is  preserved 
when  the  process  is  sampled  at  non-increasing  instances  of 
time.  This  feature  is  important  in  simulation  studies,  for 
example,  when  Runge-Kutta  integration  algorithm  is  used. 

(ii)  A  pulse  sequence  simulating  the  excitation  signals  used  in 
the  parameter  estimation  studies.  The  subroutine  allows 
adjustment  of  various  pulse  sequence  parameters  such  as 
amplitude,  pulsewidth,  period  and  start  time. 


5.2.3  Line-of-Sight  Calculations 

A  ray  tracing  program  for  the  two  mirror  Cassegrain  configuration  was 
developed.  The  program  computes  the  mirror  misalignments  as  function  of 
the  mirrors'  motion.  Two  alternative  ways  of  measuring  misalignments  are 
possible.  One  is  in  terms  of  the  ray  scatter  (r,w)  and  the  location  of 
trace's  center  on  the  focal  plane  (y  ,  zm).  The  other  method  of  measuring 
is  in  terms  of  focal  plane  perturbation  (az)  and  the  change  of  LOS  direc¬ 
tion  (ey,  ez) .  The  two  alternative  ways  of  LOS  measurement  are  illustrated 
in  Figure  5-5. 


In  order  to  avoid  time  consuming  calculations  of  LOS  motion  in  the 
simulation  program,  a  linear  approximation  was  performed.  It  was  assumed 
that  the  mirrors  motion  is  small  and  only  the  linear  terms  of  Taylor  series 
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are  needed  to  evaluate  LOS  motion.  Therefore,  the  various  perturbations 
can  be  written  as: 


aF 


[^]  ax1  +  [$21  ax2 


(5.1) 


and 


aL 


az 

9y 

0. 


AX|  +  C d; 2 3  ^x2 


(5.2) 


where: 


Af,  aL :  Perturbed  LOS  vectors 

$  ^ :  Sensitivities  of  LOS  to  mirror  1  motion 

$  Sensitivites  of  LOS  to  mirror  2  motion 

Ax  ax~:  Vector  9  x  1  expressing  the  three  dimensional 
2  motion  of  three  points  on  mirrors  1  and  2. 


The  sensitivity  matrices  were  obtained  by  perturbing  each  mirror  support 
points  separately  and  computing  the  trace  motion  on  the  focal  plane  using 
the  ray  tracing  program.  The  various  sensitivity  matrices  are  given  in 

Table  5-3. 


5.3  Simulation  Results 


5.3.1  Setting  the  Initial  Conditions 

In  order  to  create  the  worst  case  structure  excitation  caused  by  VRCS 
thruster  firings,  the  initial  conditions  of  the  Orbiter  ought  to  be  set  in 
a  way  which  will  cause  frequent  firings.  Moreover,  in  order  to  reduce  simula¬ 
tion  costs,  the  firings  should  occur  in  the  initial  stage  of  the  simulation. 
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Table  5-3.  LOS  Sensitivity  Matrices 
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Since  the  system  dynamics  is  represented  in  terms  of  modal  coordinates, 
a  procedure  was  developed  to  obtain  the  initial  conditions  of  the  modal 
coordinates  given  a  set  of  physical  initial  states.  The  procedure  assumes 
that  only  the  rigid  body  initial  conditions  are  different ‘from  zero.  Let 
q  be  the  physical  coordinates  and  nR  and  ng  be  the  rigid  body  and  flexible 
modal  coordinates,  respectively.  One  can  write: 


then 


q 


[*R  !  *ei 

I 


q(o) 


(5.3) 


(5.4) 


or 


(5.5) 


In  general,  q(o)  is  a  6  dimensional  vector  expressing  initial  displacements 
and  rotations  at  a  given  location.  However,  if  only  a  subset  of  the  physical 
coordinates  is  initialized,  only  the  corresponding  number  of  rigid  body 
modes  is  calculated.  In  our  case  only  three  modes  are  given  initial  con¬ 
ditions  because  only  the  three  angular  attitudes  are  initialized.  The 
initial  rate  is  obtained  directly  from  Equation  (5.5). 

nR(o)  =  V1  >i(0)  <5'6> 


It  was  decided  to  impose  the  largest  attitude  deviation  on  the  roll 
axis.  This  is  because  the  roll  moment  of  inertia  is  the  smallest  and  ro¬ 
tation  around  this  axis  will  conceivably  cause  the  largest  disturbances. 

The  initial  conditions  on  the  other  axes  were  set  to  also  violate  the  limit 
cycle  limits,  thus  causing  frequent  thruster  firings.  The  initial  conditions 

are  summarized  in  Table  5-4. 


Table  5-4.  Attitude  Initial  Conditions 


ROLL 

PITCH 

YAW 

ATTITUDE 

0.128 

0.11 

0.11 

(deg) 

RATE 

0.01 

0.01 

0.01 

(deg/sec ) 
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As  expected,  the  above  initial  conditions  caused  frequent  firings  in 
the  initial  stage  of  the  simulation.  The  various  vernier  thruster  firings 
durinq  the  first  20  seconds  are  shown  in  Figure  5-6.  Figure  5-7  shows  the 
attitude  of  the  Shuttle  during  the  first  50  seconds  of  a  simulation  run. 

The  behavior  shown  is  typical  of  the  Shuttle  when  disturbances  cause  the 
attitude  to  drift  outside  the  bounds  of  the  limit  cycle.  Obviously,  the 
exact  behavior  of  the  Orbiter  is  highly  dependent  on  the  initial  conditions 
and  the  disturbance^  acting  on  it  during  normal  operation.  However,  the 
behavior  presented  here  is  representati ve  and,  since  it  is  a  "worst  case" 
study,  it  gives  an  upper  limit  on  the  motion  of  the  test  structure.  As 

it  can  be  seen,  the  limit  cycle  is  irregular  and  a  definite  time  period 

is  difficult  to  observe.  This  is  due  to  the  fact  that  the  various  thrusters 
do  not  produce  a  pure  moment  around  one  axis  only.  A  compromise  is  made  by 
the  jet  selection  logic  algorithm  which  produces  a  torque  that  is  not  exactly 
in  the  direction  of  the  commanded  torque.  After  the  initial  stage  of  fre¬ 
quent  transfer  firings,  the  Orbiter  enters  the  limit  cycle  deadband  and 
coasts  without  VRCS  activity  for  a  relatively  long  period  of  time.  Therefore, 
during  normal  operation  of  the  Shuttle,  there  will  be  long  periods  of  time 
in  which  there  will  not  be  structural  excitation  caused  by  thruster  firings. 

Figure  5-8  shows  the  vibrations  induced  by  system's  flexibility  at  the 

location  of  the  I MU  sensor.  As  it  can  be  seen,  the  dynamic  interaction  of 

the  test  structure  with  the  Orbiter  does  not  cause  noticeable  vibrations. 

The  change  of  attitude  due  to  the  vibrational  motion  is  of  the  order  of 

_g 

10  rad  which  cannot  be  even  detected  by  the  IMU  sensor.  Moreover,  this 
level  of  vibration  is  within  the  modeling  error  of  the  finite  element 
model . 

Figures  5-9  and  5-10  show  the  vibrations  induced  by  the  VRCS  firings 

on  the  primary  and  secondary  mirrors.  The  motion  of  the  structure  due  to 

_2 

the  rigid  body  rotation  of  the  Orbiter  is  of  the  order  of  10  m.  The 
amplitude  of  the  vibrational  motion  excited  by  the  VRCS  thrusters  is  of 

_5 

the  order  of  10  meters.  This  is  a  relatively  low  level  of  disturbance 
which  will  not  interfere  with  the  experiments  performed  on  the  structure. 
Moreover,  for  long  periods  of  time  the  thrusters  are  not  firing  because  the 
Orbiter  is  within  the  deadband  of  the  limit  cycle.  Hence,  during  periods  of 
inactivity  the  vibrations  are  damped  out  by  the  existing  structural  damping. 
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As  can  be  observed  in  the  figures,  four  vibrational  modes  are  important: 

i)  The  bending  mode  of  the  secondary  mirror  in  the  xz  plane. 

The  vibrations  due  to  the  VRCS  limit  cycle  are  of  the  order 
10  urn. 

ii)  The  bending  mode  of  the  secondary  mirror  in  the  yz  plane. 

The  VRCS  thruster  firings  caused  vibration  with  amplitude 
up  to  65  um .  The  vibrations  in  this  direction  are  higher 
than  for  the  previous  mode  because  the  induced  roll  motion 
is  larger  than  the  pitch  motion. 

iii)  The  bending  mode  of  the  primary  mirror  in  the  xz  plane.  This 
causes  vibrations  with  an  amplitude  of  up  to  40  um. 

iv)  Combined  bending  and  torsion  of  the  primary  mirror,  which 
gives  deflections  of  up  to  10  um. 

In  conclusion,  as  it  was  predicted,  the  test  structure  behaves  like  two 
cantilevered  beams  with  concentrated  masses  at  their  ends. 

The  qualitative  behavior  of  the  system  will  not  change  for  different 
sets  of  initial  conditions.  The  change  that  might  occur  for  different 
thruster  firing  sequence' is  the  relative  deflection  of  the  mirrors.  In 
particular,  if  the  yaw  motion  is  larger,  the  vibration  mode  (iv)  will 
have  a  larger  amplitude. 

In  addition  to  the  interaction  study,  another  goal  of  the  simulation 
was  to  determine  the  level  of  excitation  needed  for  demonstrating  the 
control  system  performance.  Preliminary  calculations  indicated  that 
deflections  of  0.2  um  of  the  mirrors  will  cause  1 ine-of-sight  (LOS)  error 
of  0.8  urad.  Therefore,  in  order  to  demonstrate  vibration  suppression  by 
a  factor  of  1000,  the  imposed  excitation  should  cause  deflections  of  ^  300  um. 
Based  on  previous  results,  it  was  determined  that,  in  crder  to  excite  the 
significant  structural  modes,  the  disturbance  forces  should  be  applied  at 
the  top  of  the  mast  and  on  one  of  the  side  support  points  of  the  primary. 

The  simulation  study  determined  that  the  following  set  of  Gaussian,  band- 
limited  (25  Hz)  forces  causes  the  required  deflections: 

a)  Force  in  x  direction  acting  on  top  of  the  mast  (node  3), 

0-mean,  a£(0)  =  10"3(Hz_1). 

b)  Force  in  y  direction  acting  on  top  of  the  mast,  0-mean 
o2(0)  =  10-3(Hz-1). 
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c)  Force  in  z  direction  acting  on  one  of  the  side  supports  of  the 
primary  mirror  (node  14),  0-mean,  <j2(0)  =  10_3(Hz"l). 

Figures  5-11  and  5-12  show  the  vibrations  of  the  two  mirrors  as  the 
result  of  the  above  disturbance  forces.  The  primary  mirror's  support 
(node  14)  undergoes  vibrations  with  amplitude  up  to  100  pm  in  the  z  direc¬ 
tion,  40  urn  in  the  y  direction,  and  ^  10  pm  in  the  x  directions.  The 
secondary  mirror's  deflections  are  ^  50  m  in  the  x  direction,  <v.  40  pm  in 
the  y  direction  and  ^  10  pm  in  the  z  direction. 

The  resulting  LOS  perturbation  is  shown  in  Figure  5-13.  It  can  be 
seen  that  the  LOS  error  is  in  the  range  which  will  permit  control  law 
demonstration  of  vibration  suppression  by  a  factor  of  a  1000.  The  open 
loop  LOS  error  is  of  the  order  of  500  prad  which  is  within  the  desired 
range. 
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Secondary  Mirror  in  Presence  of  Disturbance  Forces 


Table  6-1.  Candidate  List  of  Experiments 


I.  COMPONENT  TESTS  (GROUND) 

OBJECTIVE:  TO  DETERMINE  ANALYTICAL  MOOELS  AT  THE  POC  STRUCTURE  COMPONENT  LEVEL. 

SET  OF  TESTS; 

A.  DETERMINE  TRUSS  MATERIAL  PROPERTIES  AND  GEOMETRY  NECESSARY  TD  DEVELOP  FINITE  ELEMENT  .MODELS  (FEM)  AT 
EXPECTED  TEST  AND  OPERATING  ENVIRONMENT  (E.G.,  TEMPERATURE,  PRESSURE,  RADIATION,  ETC.). 

B.  DETERMINE  ACTUATOR  AND  SENSOR  DYNAMICS  OVER  OPERATING  RANGE 

1.  LINEARITY 

2.  RESOLUTION 

3.  DYNAMIC  RANGE 

4.  NOISE  CHARACTERISTICS 

5.  BANDWIDTH  ANO  FREQUENCY  RESPONSE  AS  A  FUNCTION  OF  SIGNAL  AMPLITUDE  IN  THE  OPERATING  RANGE 

C.  OETERMINE  ENVIRONMENTAL  DISTURBANCES  THAT  WILL  NOT  PARTICIPATE  IN  PDC  TEST  AND  "RECORD" 

1.  CRYOCDOLER  "NOISE"  AND  COOLANT  MASS  SHIFTS 

2.  CMG  IMBALANCE  AND  NOISE 

3.  RW  IMBALANCE  AND  NOISE 

4.  SOLAR  PANEL  MOTION  AND  SLIPRING  NOISE 

5.  IMPULSIVE  DISTURBANCES 


SUBSTRUCTURE  TESTS  (GROUND) 


OBJECTIVE:  TO  DETERMINE  ANALYTICAL  MODELS  AT  THE  SUBSTRUCTURE  LEVEL. 


SET  OF  TESTS: 

A.  SHUTTLE  PALLET 

1.  I.D.  MODAL  DATA 

2.  I.O.  OAMPING 

3.  COMPARE  "REAL"  VERSUS  FEM  MOOEL 

B.  PRIMARY  MIRROR  +  SUPPORT 

1,  2,  AND  3.  SAME  AS  II. A. 

C.  TEST  STRUCTURE 

I,  2,  AND  3.  SAME  AS  II. A. 

D.  MISCELLANEOUS  FLEXIBLE  (<  XIO  CONTROL  BW)  COMPONENTS  TO  BE  ATTACHED  TO  TEST  STRUCTURE 
I,  2,  AND  3.  SAME  AS  II. A. 


III.  SYSTEM  TESTS  (GROUNO,  NO  SHUTTLE  ANO  ON-ORBIT) 

OBJECTIVE:  TO  CALIBRATE  THE  SYSTEM  AND  TO  VALIDATE  OPEN  ANO  CLOSED  LOOP  SYSTEM  MOOELING. 

A.  CALIBRATE  OPTICAL  MEASUREMENT  ANO  ALIGNMENT  SYSTEM 

OBJECTIVE:  ESTABLISH  NOMINAL  CAPABILITY  AND  OPERATING  POINT  FOR  DISTURBANCE  FREE  SYSTEM 

METHOD:  APPLY  KNOWN  ROTATIONAL  AND  TRANSLATIONAL  PERTURBATIONS  TO  ELEMENTS  IN  OPTICAL  LOS  PATH,  TUEN 
COMPARE  MEASUREO  COMPONENT  MOTION  VS  REAL.  USE  TRW  SURFACE  ACCURACY  MEASURMEENT  SENSOR  (SAMS)  IN 
TRIANGULATION  MODE.  THIS  SENSOR  WILL  MONITOR  SIX-OEGREE-OF-FREEDOM  MOTION  OF  RIGID  POINTS  ON  PRIMARY 
AND  SECONDARY  MIRROR  STRUCTURE  (AND  ANY  OTHER  ELEMENTS  ALONG  THE  OPTICAL  PATH)  IN  FOCAL  PLANE  COORDINATES. 


6-2 


Table  6-1.  Candidate  List  of  Experiments  (Continued) 


CALIBRATE  IMAGE  MOTION  SENSDR  SYSTEM 

OBJECTIVE:  ESTABLISH  NOMINAL  CAPABILITY  AND  OPERATING  POINT  FOR  DISTURBANCE  FREE  SYSTEM.  ALSO,  TEST  AND 
CALIBRATE  OPTICAL  ELE^NT  MOTION  TO  LOS  MOTION  TRANSFORMATION. 

METHDD-  USE  IMAGING  CHARGE  COUPLED  DEVICE  LIKE  USED  IN  TRW'S  MADAN  AND  FOCUS  CENTROID  OF  “TRUTH"  TARGET; 
THEN  APPLY  KNOWN  ROTATIONAL  AND  TRANSLATIONAL  PERTURBATIONS  TO  ELEMENTS  IN  OPTICAL  LOS  PATH  AND  COMPARE 
PREDICTED  LOS  MDTION  VERSUS  REAL. 

CALIBRATE  ACTUATOR  TO  SENSDR  "DC"  TRANSFORMATION 

OBJECTIVE:  CALIBRATE  AND  VALIDATE  THE  INFLUENCE  THE  ACTUATORS  (E.G.,  CONTROL,  DISTURBANCE,  I.D.)  HAVE  ON 
LOS,  OPTICAL  ELEMENT  MDTIDN,  AND  STRUCTURE. 

METHDO:  APPLY  KNOWN  STATIC  PERTURBATIONS  TO  EACH  ACTUATOR  IN  A  SEQUENTIAL  MANNER  AND  MONITOR  DC  CHANGE 
ON  ALL  SENSORS.  COMPARE  PROJECTED  VERSUS  REAL. 

MODAL  SYNTHESIS  EXPERIMENT,  PHASE  I 

OBJECTIVE:  VALIOATE  ANALYTICAL  MODAL  SYNTHESIS  PROCEDURES  AND  ESTABLISH  A  REFERENCE  FOR  THE  METHOD'S 
MOOELING  ERROR  PROPAGATION. 

CONDITIONS:  GROUND,  NO  SHUTTLE,  CONTROL  OFF,  PASSIVE  CONTROL  OFF  (MEMBER  DAMPERS  AND  ISOLATORS  "LOCKEO"), 
DISTURBANCES  OFF. 

METHOD: 

1.  PLACE  ASSEMBLED  STRUCTURE  PLUS  PALLET  (IF  POSSIBLE)  ON  TEST  BED,  AND  IDENTIFY  MODAL  DATA  AND  DAMPING 

2.  COMPARE  WITH  FEM  MODEL,  EXPERIMENT  I. A. 

3.  COMPARE  WITH  MODAL  SYNTHESIS  FEM  MODEL 

4.  COMPARE  WITH  MODAL  SYNTHESIS  USING  "REAL"  SUBSTRUCTURE  DATA  -  EXPERIMENT  II. 

DISTURBANCE  RECORDING  VALIDATION 

OBJECTIVE:  VERIFY  THAT  OISTURBANCE  "RECORDING"  IS  REPRESENTATIVE  OF  EXPECTED  DISTURBANCES. 

CONDITIONS:  SAME  AS  III.2.D  BUT  WITH  DISTURBANCES  ON  AS  REQUIRED. 

METHOD:  TURN  ON  DISTURBANCE  ACTUATORS  IN  A  SEQUENTIAL  MANNER  AND  MONITOR  EFFECT  USING  ALL  SENSORS. 

COMPARE  WITH  ANALYTICAL  PREDICTION. 

VERIFICATION  OF  LUMPED  PASSIVE  ELEMENT  MODELING/PASSIVE  CONTROL  EXPERIMENT,  PHASE  I 

OBJECTIVE:  VERIFY  ANALYTICAL  METHODS  USED  TO  MODEL  STRUCTURES  WITH  INTEGRATED,  LUMPED,  PASSIVE  ELEMENTS. 
CONDITIONS:  GROUND,  NO  SHUTTLE,  CONTROL  OFF,  PASSIVE  CONTROL  ON  AS  REQUIRED,  DISTURBANCES  ON  AS  REQUIRED. 
METHOD: 

1.  VALIDATE  EFFECT  OF  PASSIVE  ISOLATOR  MODELING;  MEMBER  DAMPERS  LOCKED,  ISOLATORS  UNLOCKED. 

ID  WITH  AND  WITHOUT  DISTURBANCES.  MONITOR  LOS  AND  STRUCTURAL  DEFLECTIONS  AND  COMPARE  WITH  PREDICTED. 

2.  VALIOATE  EFFECT  OF  PASSIVE  MEMBER  DAMPERS  MODELING:  MEMBER  DAMPERS  UNLOCKED,  ISOLATORS  LOCKED. 

ID  WITH  AND  WITHOUT  DISTURBANCES.  MONITOR  LOS  AND  STRUCTURAL  DEFLECTIONS  AND  COMPARE  WITH  PREDICTEO. 

3.  VALIDATE  PASSIVE  CONTROL  SYSTEM  PERFORMANCE  PREDICTION:  MEMBER  DAMPERS  AND  ISOLATORS  UNLOCKED.  ID 
WITH  AND  WITHOUT  DISTURBANCES.  MONITOR  LOS  AND  STRUCTURAL  DEFLECTIONS  AND  COMPARE  WITH  PREDICTED. 
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Table  6-1.  Candidate  List  of  Experiments  (Continued) 

G.  VERIFICATION  OF  ACTIVE  STRUCTURAL  CONTROL  MODELING 

OBJECTIVE:  VERIFY  ANALYTICAL  METHODS  USED  TO  MODEL  ACTIVE  CONTROL 

CONDITIONS:  GROUND,  NO  SHUTTLE,  PASSIVE  CONTROL  ON,  ACTIVE  STRUCTURAL  CONTROL  ON  AS  REQUIRED,  DISTURBANCES 
ON  AS  REQUIRED. 

METHOD: 

1.  COMPUTE  SIMPLE  ACTIVE  LAWS  (E.G.,  PIB)  FOR  FEM  MODEL  (I. A.),  “REAL"  MODAL  SYNTHESIZED  MODEL  (II. D), 

ANO  "REAL"  SYSTEM  (III.F.3). 

2.  FOR  EACH  SET  OF  GAINS  ABOVE:  ID  WITHOUT  DISTURBANCES,  MEASURE  STABILITY  MARGINS,  MONITOR  LOS  AND 
STRUCTURAL  DEFLECTIONS  WITH  DISTURBANCES. 

3.  COMPARE  WITH  ANALYTICAL  PREDICTION  AND  III.F.3. 

IV.  TECHNOLOGY  DEMONSTRATIONS  (GROUND  WITH  NO  SHUTTLE  AND/OR  ON-ORBIT  WITH  SHUTTLE  AS  REQUIRED! 

A.  ADVANCED  CONTROL  LAWS 

OBJECTIVE:  DETERMINE  STABILITY  AND  PERFORMANCE  LIMITATIONS  OF  ADVANCED  CONTROL  LAWS 

CONDITIONS:  PASSIVE  CONTRDL  OH  AS  REQUIRED,  ACTIVE  CONTROL  ON  AS  REQUIRED,  DISTURBANCES  ON  AS  REQUIRED. 

METHOD: 

1.  IDENTIFY  MODAL  DATA  AND  DAMPING  OF  UNCONTROLLED  STRUCTURE  WITH  NO  DISTURBANCES. 

2.  MEASURE  UNCONTROLLED  SYSTEM  RESPONSE  TO  DISTURBANCES. 

3.  FOR  COMPETING  FIXED  GAIN  COLOCATED  CONTROL  LAWS,  NON-COLOCATED  CONTROL  LAWS,  DIRECT  ADAPTIVE, 

AND  FIXED  GAIN  LAW  BASED  ON  “REAL"  MODAL  DATA  (INDIRECT  ADAPTIVE): 

A.  DETERMINE  STABILITY  ROBUSTNESS  BY  INSERTING  GAIN  AND  PHASE  PERTURBATIONS  INTO  THE  CONTROL  LOOP. 

3.  DETERMINE  PERFORMANCE  BY  MEASURING  LOS  AND  OPTICAL  TRAIN  RESPONSE  TO  DISTURBANCE  ENVIRONMENT. 

4.  DETERMINE  PERFORMANCE  (AS  ABOVE)  USING  RELATIVE  AND  INERTIAL  CONTROL  LAWS  DESIGNED  USING  THE  SAME 
CONTROL  PHILOSOPHY  AND  MODEL  INFORMATION. 

B.  WAVE  PROPAGATION  EXPERIMENT 

OBJECTIVE:  EVALUATE  PERFORMANCE  AND  RELATIVE  MERITS  BETWEEN  CONTROLLING  AT  SOURCE,  CONTROLLING  ALONG 
PROPAGATION  PATH  (STRUCTURAL  CONTROL)  AND  CONTROLLING  AT  DESTINATION  (HIGH  8ANDWIDTH  OPTICAL  AUTOALIGNMENT 
SYSTEM).  ALSO,  ESTABLISH  PERFORMANCE  LIMITATIONS. 

METHOD: 

THE  PERFORMANCE  OF  THE  THREE  CONTROL  APPROACHES  WILL  BE  COMPAREO  INDIVIDUALLY  AND  IN  COMBINATION.  LOS 
ANO  OPTICAL  TRAIN  RESPONSE  TO  DISTURBANCES  WILL  BE  THE  CRITERIA. 

C.  IMAGING  EXPERIMENT 

OBJECTIVE:  EVALUATE  PERFORMANCE  AND  RELATIVE  MERITS  BETWEEN  DITHERED  ACTIVE  OPTICS,  WAVEFRONT  COMPENSA¬ 
TION,  AND  ACTIVE  STRUCTURAL  CONTROL.  ALSO,  ESTABLISH  PERFORMANCE  LIMITATIONS. 

METHOD: 

IMAGE  QUALITY  WILL  BE  EVALUATED  USING  A  TRUTH  TARGET  AND  EACH  OF  THE  THREE  IMAGE  CONTROL  TECHNIQUES  ABOVE, 
ALSO  WITH  NO  CDNTROL .  ACTIVE  STRUCTURAL  CONTROL  WILL  THEN  BE  USED  IN  COMBINATION  WITH  ONE  OF  THE  OPTICAL 
METHOOS  TO  ESTABLISH  PERFORMANCE  LIMITS. 
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Table  6-1.  Candidate  List  of  Experiments  (Continued) 


•LLANEOUS  (ON-ORBIT  WITH  SHUTTLE) 

POC  STRUCTURE  DYNAMICAL  INTERACTION  WITH  ORBITER 

OBJECTIVE:  VALIDATE  ANALYTICAL  MOOELS  FOR  OYNAMICAL  INTERACTION  BETWEEN  PRECISION  STRUCTURE  ANO  HOST 
SPACECRAFT.  ALSO.  ESTABLISH  THE  SHUTTLE  ENVIRONMENT  MOOEL. 

METHOO: 

1.  MONITOR  OPEN  LOOP,  DISTURBANCE  FREE  OE FLECTIONS  OF  POC  STRUCTURE  AS  SHUTTLE  IS  IN  NORMAL  OPERATION. 
ALSO  MONITOR  SHUTTLE  CONTROL  SIGNALS. 

2.  MONITOR  CLOSEO  LOOP,  DISTURBANCE  OEFLECT IONS  OF  POC  STRUCTURE  AS  SHUTTLE  IS  IN  NORMAL  OPERATION. 
MONITOR  SHUTTLE  CONTROL  SIGNALS. 

RAP  10  SLEW-TRACK-SETTLE 

OBJECTIVE:  VALIOATE  RAPIO  SLEW-TRACK-SETTLE  ALGORITHMS  REQUIREO  IN  SOME  MISSIONS. 

METHOO: 

1.  USE  SHUTTLE  TO  PERFORM  OPTIMUM  OPEN  LOOP  SLEW. 

2.  COUPLE  SIGNALS  FROM  SHUTTLE  AUTOPILOT  INTO  POC  OPTICAL  ALIGNMENT  SYSTEM  ANO  EXECUTE  FINE  SLEW  CONTROL 
AND  TRACK  WITH  ACTUATEO  OPTICAL  ELEMENTS.  SETTLE  BY  COUPLING  SLEW  ANO  TRACK  SIGNALS  INTO  STRUCTURAL 
CONTROLLER. 


The  objective  of  these  sequence  of  tests  is  to  validate  open  and 
closed  loop  structural  control  models. 

The  need  to  validate  open-loop  structural  models  is  self  evident. 
To  do  it,  all  actuators  and  sensors  are  first  calibrated,  then  the 
test  structure  is  identified.  The  resulting  test  data  is  compared 
with  the  predictions  made  with  the  two  modal  synthesis  models  from  the 
substructure  tests. 


The  need  to  validate  closed-loop  structured  models  at  sub  pm  and 
yrad  levels  is  also  of  utmost  importance.  Up  until  the  present  time 
only  the  rigid  body  modes  of  a  spacecraft  have  actively  been  controlled 
All  of  the  flexible  modes  have  only  been  stabilized  (mainly  gain 
stabilized).  Hence  it  hasn't  been  too  important  whether  closed- 
loop  flexible  modes  act  according  to  linear,  lumped  parameter  feed¬ 
back  theory  or  not.  With  the  current  emphasis  on  precision  structural 
control  however,  there  is  a  need  to  determine  this  information.  It 
is  expected  that  the  baseline,  truss  like  system  will  act  as  the 
mechanical  analog  of  a  transmission  line  with  wave  delays  and  energy 
absorbing  "stubs"  (active  and  passive  control). 

The  proposed  validation  experiment  is  to  first  verify  passive  con¬ 
trol  modelling.  Here  the  actual  open  and  closed  loop  responses  of 
the  system  to  known  disturbance  inputs  are  compared  with  those 
predicted  by  linear  feedback  models.  Once  the  analytic  feedback  model 
has  been  validated,  active  control  models  will  be  evaluated  in  a 
similar  fashion.  Passive  control  is  described  first  so  that  a  gross 
error  in  the  analytical  model  will  result  in  only  a  benign  miscalcu¬ 
lation  rather  than  a  potential  hardware  malfunction. 

6.1.4  Technology  Demonstrations 

The  technology  demonstration  consists  of  these  parts:  advanced 
control  law  validation  and  trades,  wave  propagation  experiment;  and 
imaging  experiments. 


The  objective  of  the  advanced  control  law  demonstration  is  to  com¬ 
pare  the  various  colocated,  noncolocated,  fixed  gain  and  adaptive 
control  laws  developed  for  active  control.  Performance,  as  measured 
by  LOS  error,  and  stability  margins  are  the  key  criteria.  Comparing 
the  performance  of  relative  (eg.,  bending  moment)  and  inertial 
(eg.,  proof-mass)  control  systems  is  also  an  objective.  It  is  noted 
that  only  one  or  two  of  the  most  promising  approaches  are  to  be  demon¬ 
strated  on-orbit;  the  remaining  will  be  discarded  based  on  the  results 
of  the  ground  test. 

The  objective  of  the  wave  propagation  experiment  is  to  determine 
what  type  or  combination  of  control  approaches  is  best  suited  for 
optical  missions  with  high  disturbance  sources.  Currently,  one  can 
control  at  the  disturbance  source  location  by  isolating  it  before  its 
energy  gets  into  the  structure;  one  can  control  along  the  path  of 
energy  propagation  through  the  structure  with  passive  and  active 
vibration  suppression  techniques;  or  one  can  control  at  the  juncture 
between  the  optical  train  and  the  structure  by  using  a  high  bandwidth 
autoalignment  system.  Trades  between  these  systems  are  possible 
in  terms  of  performance,  cost,  weight,  complexity,  and  power.  Again, 
only  the  most  promising  approach  is  intended  for  on-orbit  test.  The 
performance  criteria  is  optical  train  motion  and  LOS  error. 

The  final  technology  demonstration  and  trade  is  the  imaging 
experiment.  The  objective  here  is  to  experimentally  determine  the 
relative  merits  between  active  optics  (dithered  or  wavefront  compen¬ 
sation)  and  active  structural  control.  The  intent  is  to  determine  the 
most  advantageous,  individual  areas  of  application  and  combination. 
Image  quality  is  to  be  evaluated  by  using  truth  targets.  Only  the 
most  promising  combination  is  to  be  evaluated  on-orbit. 

6.1.5  Miscellaneous 

This  set  of  experiments  addresses  the  dynamical  interaction 
between  a  large  structure  (Shuttle)  and  a  payload,  and,  also 
slew-track-settle  maneuvers. 


The  objectives  of  the  dynamical  interaction  experiments  are  to 
i)  model  the  Shuttle  environment  prior  to  performing  any  of  the 
on-orbit  experiments;  and  ii)  test  the  analytical  prediction 
capability  of  modal  synthesis  between  two,  dynamically  complex 


structures . 


The  objectives  of  the  rapid  slew-track-settle  experiment  are  to 
test  some  control  coordination  algorithms  required  in  some  missions 
The  algorithm  arises  when  one  tries  to  coordinate  a  slew  maneuver 
(large  dynamic  range  and  small  bandwidth}  with  structural  control 
(medium  dynamic  range  and  bandwidth)  and  figure  control  (small 
dynamic  range  and  high  bandwidth).  This  control  algorithm  is 
necessarily  non-linear,  multivariable,  and  distributed  in  nature. 
Hence,  it  needs  to  be  tested.  The  on-orbit  experiment  presents  a 
good  opportunity  to  do  this  since  the  Shuttle  itself  can  be  used  in 
place  of  the  typically  very  expensive  attitude  control  system. 

6.2  Relative  Priority  Between  Experiments 


Table  6-2  check  lists  the  set  of  candidate  experiment.-  according 
to  whether  they  have  primary  or  secondary  importance  in  a  modeling 
intensive  program,  a  technology  verification  intensive  program  or 
both.  The  "recommended"  set  is  based  on  the  analysts  understanding 
of  the  current  state  of  the  technology.  The  actual  experiments  and 
procedures  chosen  from  this  candidate  list  are  described  in  Appendix  A. 
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6.3  The  Control  Subsystem 

The  control  subsystem  was  designed  to  perform  the  various  forms  of 
control  required  by  the  experiments  sequence  described  in  Section  6.1.  In 
order  to  be  able  to  implement  the  required  control  forms,  the  control  sub¬ 
system  was  divided  into  three  functionally  independent  control  modules.  Eac 
module  is  to  be  activated  or  disabled  according  to  the  experiment  to  be 
performed.  This  section  also  establishes  the  baseline  location  of  actuators 
and  sensors  for  each  control  module  as  well  as  hardware  and  software 
requirements  to  implement  the  defined  control  laws. 

6*3.1  Control  System  Functions 

The  functional  diagram  of  the  control  subsystem  is  shown  in  Figure  6-1. 
The  three  parts  of  the  control  systems  are  shown  in  the  figure  as 
fol  lows  : 

(i)  Passive  Control  System:  This  system  consists  of  a  series  of 
dampers  and  isolators  which  can  be  locked  and  unlocked  at  any  given  time. 
These  passive  elements  are  placed  in  locations  at  which  maximum  energy 
dissipation,  and  thus  maximum  vibration  suppressors,  occurs.  The  passive 
control  system  will  be  used  to  demonstrate  the  ability  to  control  the 
structure  at  the  source  of  the  disturbance. 

(ii)  Autoal ignment  System:  The  autoalignment  system  consists  of 
actuators  which  adjust  the  position  of  the  mirrors  to  minimize  the  LOS 
disturbances  as  measured  by  the  LOS  and  figure  sensors.  The  system  has 
low  and  high  bandwidth  operating  modes.  The  high  bandwidth  mode  will  be 
used  to  demonstrate  the  performance  of  the  control  subsystem  while 
controlling  the  system  at  the  destination  during  the  wave  propagation 
experiment. 


The  low  bandwidth  operation  mode  will  be  used  in  conjunction  with  other 
control  modules,  e.g.  passive  control,  to  assess  the  performance  of 
a  control  approach  combination. 

(iii)  Vibration  Suppression  System:  This  module  is  divided  according 


to  the  control  approach  into  either  a  inertial  or  relative  structural  control 
subsystems.  The  inertial  system  measures  the  absolute  deflections  at  various 


Figure  6-1.  Control  Subsystem 


points  and  applies  inertial  forces  to  minimize  the  disturbance.  Whereas 
the  relative  system  measures  the  relative  slopes  between  points  on  the 
structure  and  applies  a  set  of  bending  moments  to  counteract  the  dis¬ 
turbance  forces  acting  on  the  system.  The  vibration  suppression  system 
operating  in  either  the  inertial  or  relative  operating  modes  will  be  used 
to  demonstrate  controlling  the  system  along  the  propagation  path  (structural 
control)  during  the  wave  propagation  experiment. 

6.3.2  Actuator  and  Sensor  Locations 

The  modal  participation  factor  (g-factor)  was  used  as  a  criterion  for 
establishing  the  locations  of  the  various  actuators  and  sensors.  The 
inertial  g-factor  for  mode  i  at  location  j  is  defined  by: 


The  relative  motion  g-factor  between  locations  j  and  k  for  mode  i  is  given 
by: 


_  (*i.i  ~ 


(5.2) 


where  <j>  is  the  mode  shape, 


A  g-factor  analysis  was  performed  on  the  test  structure  using  a  TRW 
developed  program  MDLSFT.  The  program  orders  the  nodes  according  to  their 
g-factors,  one  ordering  for  each  one  of  the  modes  and,  if  required,  orders 
the  modes  according  to  the  magnitude  of  the  g-factor  at  a  given  node.  The 
program  performs  this  task  in  an  interactive  fashion  and  can  calculate  both 
inertial  and  relative  g-factors  as  required. 


The  locations  of  the  actuators  and  sensors  for  the  different  control 
modules  are  given  next.  The  node  locations  are  defined  in  Figure  5-3. 


The  inertial  control  system  consists  of  linear  momentum  exchange  devices 
located  on  the  top  of  the  mast  (node  3)  and  at  two  support  points  of  the 
primary  mirror.  The  actuators  located  at  the  top  of  the  mast  act  in  the 


x  and  y  directions,  and  the  actuators  at  nodes  14  and  15  act  in  the 
z-direction.  Sensors  in  the  form  of  accelerators  are  placed  next  to  the 
actuators,  creating  a  colocated  system. 

The  relative  control  system  consists  of  four  bending  moment  actuators. 

Two  actuators  are  connected  between  nodes  17  and  18  and  supply  torques  T 
and  Tx  to  control  bending  in  the  x  and  y-directians ,  respectively.  The  other 
two  actuators  are  connected  between  nodes  10  and  14,  and  10  and  15,  supplying 
moments  around  the  x  axis,  thus  controlling  the  z  deflection  of  the  primary 

mi rror .  Sensors  are  located  at  the  attachment  points  of  the  bending  actuators, 
creating  a  colocated  system. 

The  autoalignment  system  actuators  are  located  at  the  supports  of  the 
two  mirrors.  Since  each  mirror  has  a  kinematic  mount,  the  mirrors'  orien¬ 
tation  can  be  adjusted  by  moving  the  support  points  without  causing  stresses 
in  the  structure. 


The  passive  system  includes  isolators  and  dampers  connected  at  the 
connection  points  of  the  relative  control  system  i.e.,  nodes  17  and  18  for 
bending  in  the  xz  and  yz  planes  and  between  joints  10  and  15,  and  10  and 
14.  The  passive  system  is  tuned  to  dissipate  energy  from  low  frequency 
modes:  the  first  two  bending  modes  of  the  mast  and  the  bending  and  torsional 
modes  of  the  primary  mirror  and  its  support. 

6.3.3  Hardware  Requirements 

A  preliminary  study  of  the  hardware  requirements  for  the  various  control 
loops  was  also  performed.  The  results  of  this  study  are  as  follows: 

^  Inertial  Control  System:  The  baseline  inertial  control  system's 
actuators  are  four  linear  momentum  exchange  devices.  Each  actuator  is  base¬ 
lined  to  have  a  bandwidth  of  Wg  =  150  Hz  with  dynamic  range  +  2N,  accuracy 
of  +0.02N,  and  resolution  of  0.02N.  Sensors  in  the  form  of  piezoelectric 
accelerometers  will  be  placed  next  to  the  actuators  and  thus  will  create 
a  colocated  control  system.  The  sensors  are  baselined  to  have  a  bandwidth 
of  200  Hz,  range  +50g,  and  noise  range  of~0.0005g  rms. 


(ii)  Relative  Control  System:  The  relative  control  system  will  consist 
of  four  bending  moment  actuators  located  in  the  places  specified  earlier. 

The  actuator  characteristics  are  wg%200  Hz,  dynamic  range  +  2Nm,  and  a 
resolution  of  +0.03  Nm.  The  sensors  colocated  with  the  actuators  will  have 
the  following  characteristics:  bandwidth  a>g=  200  Hz,  dynamic  range  of 

— 1 500  urn,  and  resolution  of~0.1  pm. 

(iii)  Autoal ignment  System:  The  autoalignment  actuating  system  will 
consist  of  a  stepper  motor  operating  in  a  microstepping  operation  mode, 

100:1  harmonic  gear  drive,  and  a  sensitive  ball  screw  system.  An  advan¬ 
tage  of  the  microstepper  is  its  easy  interface  with  a  digital  controller. 

The  autoalignment  assembly  is  baselined  to  have  a  bandwidth  of  100  Hz, 

a  motion  resolution  of  0.1  urn,  and  a  dynamic  range  of  — 1 500  urn.  Since  the 
autoalignment  actuators  are  located  at  the  support  points  of  the  mirrors, 
the  control  system  will  consist  o^  six  control  loops. 

6.3.4  Software  Requirements 

At  the  last  stage  of  the  control  subsystem  design  study,  the  software 
requirements  were  estimated.  These  requirements  are  summarized  in  Table 
6-3.  Software  estimates  are  based  on  the  assumption  that  each  control 
loop  contains  an  n-th  order  filter  with  up  to  n-poles  and  n-zeroes.  It  is 
assumed  that  each  filter  is  implemented  using  parallel  realization  which 
consist  of  second  and  first  order  filters.  The  assumptions  on  the  number  of 
operations  of  each  filter  were  as  follows:  n  first  order  and  second  order 
filters  are  as  follows: 

n  -  first  order  filters 


3n+l  Multiplications 

3n+l  Additions 

8n  Logical  instructions 

4n+2  Storage  locations  for  constants 

8n  Storage  locations  for  instructions 


n  -  second  order  filters 


5n+l  Multiplications 

5n+l  Additions 

16n  Logical  instructions 

9n+2  Storage  locations  for  constants 

16n  Storage  locations  for  instructions 

Assuming  fixed  point  operations,  the  following  assumptions  are  made 
concerning  the  timing  of  each  operation: 

Addition  1  cycle 

Multiply  4  cycles 

Logical  instruction  1  cycle 
Division  30  cycles 

Trigonometric  function  50  cycles 

It  is  assumed  that  each  cycle  takes  1  usee.  The  A/D  converter's  time 
is  taken  as  8  usee  and  the  ready  and  transfer  time  is  assumed  to  be  2  usee 
The  D/A  converters  delay  time  is  assumed  to  be  3  usee. 
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The  parameter  estimation  subsystem  serves  many  functions  of  the 
POC  experiment.  They  are: 

1)  To  provide  an  accurate  on-orbit  test  structure  model  which 
can  be  related  to  analytical  and  ground  test  results. 

These  relationships  can  help  one  learn  how  to  model  future 
large  space  structure  more  accurately. 

2)  To  predict  the  performance  and  stability  of  high  technology 
demonstration  experiments  before  conducting  them,  thus  ensuring 
the  success  and  safety  of  the  tests. 

3)  To  validate  the  analytical  method  used  to  model  passive  and 
active  control  systems  by  considering  the  control  systems  as 
integral  part  of  the  structure  whose  damping  and  other  dynamic 
properties  are  being  determined. 

Parameter  estimation  algorithms  for  LSS  applications  are  already  analyzec 
in  our  previous  (ACOSS-8)  study  [19],  which  identifies  the  maximum  likeli¬ 
hood  method  as  the  baseline  approach  and  the  frequency  response  method  as 
a  way  of  providing  redundant  information  for  cross-checking.  Both  methods 
are  reviewed  briefly  in  this  report.  Two  important  advances  in  the  maximum 
likelihood  method  are  made  during  the  current  study.  One  is  the  develop¬ 
ment  of  a  substructure  identification  technique  which  allows  one  to  identify 
the  test  structure  in  spite  of  Orbiter  attachment  and  movement  .  The  other 
is  the  utilization  of  a  general  purpose  finite  element  analysis  program  in 
the  identification  loop;  this  allows  the  parameter  estimation  software  to  be 
easily  adapted  to  any  other  structure. 

Simulation  studies  of  applying  these  methods  to  the  ACOSS-14  test  struc¬ 
ture  are  conducted  in  order  to  obtain  design  data.  Results  show  that  the 
Orbiter  dynamic  interaction  has  minimal  effect  on  the  parameter  estimation 
subsystem.  The  maximum  likelihood  method,  as  expected,  gives  very  good 
estimates  and  the  frequency  response  method  yields  good  frequency  estimates 
but  is  coarse  on  damping  and  mode  shapes.  The  selected  actuators  and 
sensors  are  found  to  be  able  to  generate  adequate  test  data  for  parameter 
estimation.  Most  of  the  actuators  and  sensors  are  the  same  ones  used  in  the 


Let  z(t)  be  observed  actual  response  of  a  LSS  during  time  interval. 
[0,T].  The  closeness  of  model  predicted  response  to  the  actual  response 
is  measured  by  the  negative  log  likelihood  function: 

J<9’  T)  =  J  /TCz(t)  -  y(t,  0)]T  G'1  [ z ( t )  -  y(t,  e)]  dt 

The  unknown  parameter  0  is  adjusted  so  as  to  minimize  the  above  function. 

The  minimization  is  currently  being  performed  by  the  modified  Newton-Raphson 
method. 


Physical  and  Modal  Parameters 

The  unknown  parameters  of  a  LSS  are  directly  related  to  its  model 
representation.  The  LSS  model  is  usually  obtained  through  the  finite 
element  analysis  which  results  in  two  equivalent  ways  of  describing  the 
LSS  dynamics,  each  leading  to  a  different  set  of  parameters. 

One  model  is  the  physical  representation  described  in  terms  of  the  mass 
matrix  M  and  the  stiffness  matrix  K.  The  parameters  of  this  model  is  natural  1 
the  elements  of  the  M  and  K  matrices,  e.g..  Young's  modulus,  Poisson  ratio 
and  lumped  masses.  The  advantage  of  physical  parameter  is  that  one  can 
identify  the  whole  structure  by  estimating  only  a  small  number  of  parameters. 
But  the  estimation  process  is  rather  complicated  since  it  iterates  upon  a 
finite  element  model  of  the  LSS.  In  the  current  study,  computer  softwares 
are  developed  to  incorporate  existing  general  purpose  finite  element  analysis 
program  in  the  identification  loop.  This  reduces  some  of  the  difficulties 
in  estimating  the  physical  parameters. 

The  other  model  is  the  modal  representation  obtained  through,  an  eigen¬ 
value  transformation  from  the  physical  representation.  The  parameters 
associated  with  this  model  are  modal  frequencies  and  mode  shapes,  called 
modal  parameters.  The  estimation  procedure  for  this  set  of  parameters  is 
much  simpler  but  the  large  number  of  parameters  restricts  one  to  identify 
only  a  few  critical  modes. 


Combining  above  two  equations,  the  dynamic  equation  for  the  test  structure 
in  terms  of  the  new  variables  becomes 


M2  x2  +  K22  x2  =  F2  +  M2  K-l  K2j  Xj 

The  above  equation  is  independent  of  the  Orbiter  (xQ).  If  the  accelerations 
at  the  connecting  point  (x^)  are  measured,  one  can  consider  x^  as  an  input 
to  the  test  structure  and  use  the  above  equation  to  identify  the  test 
structure  parameters. 


6. 4. 1.3  Frequency  Response  Method 


Many  versions  of  frequency  response  method  exist  in  the  literature. 
The  one  used  in  this  study  is  as  follows: 

1)  Establish  the  frequency  response  curve  of  an  actuator-sensor 
pair  by  using  Fast  Fourier  Transforms  (FFT) 


H(oj) 


FFT  (Output  Signal) 
FFT  (Input  Signal) 


2)  Assume  the  transfer  function  of  a  particular  bending  mode  is 
of  the  form  (rate  sensor) 


joj  G. 

H.M  -  — . ,  1 - ? 

-w  +  J  2%.  OJ.  (1)  +  CJ.j 


where  =  transfer  function  gain;  ^  =  damping  ratio; 
co^  =  modal  frequency 

3)  For  each  sufficiently  separated  spike  in  the  response  curve, 
measure  the  frequency  and  the  magnitude  of  the  spike  M.. . 
Then 

“1  =  “mi  and  Mi  '  V2?i  “1 


4) 


Substitute  above  equations  into  H^(oj)  and  determine  the  damping 
ratio  c-  by  minimizing  the  cost  function 

J(?i}  =/W)|-|Hi(.)|]^ 

“1 


where  oi,  and  oi9  are  chosen  so  that  only  one  mode  is  involved. 


purpose  finite-element  analysis  program  in  the  identification  algorithm. 

The  block  diagram  in  Figure  6-2  shows  how  such  a  program  (TRWSAP,  TRW 
Structural  Analysis  Program)  is  used. 

® ^ • 2 • 2  Selection  of  Actuators  and  Sensors  for  Parameter  Estimation 

The  parameter  estimation  subsystem  should  make  as  much  use  as  possible 
of  the  actuators  and  sensors  already  selected  for  the  active  control  sub¬ 
system.  The  criterion  of  placing  actuators  and  sensors  for  parameter 
estimation  is  very  close  to  the  one  for  active  control,  i.e.,  they  should 
be  placed  at  locations  where  all  critical  modes  can  be  excited  and  measured. 
Additional  actuators  or  sensors,  however,  may  be  required  to  enhance  the 
identifiabil ity  of  certain  parameters. 

The  locations  of  actuators  and  sensors  for  parameter  estimation  are 
shown  in  Figure  6-3  and  are  also  listed  in  Table  6-4.  Other  than  the  follow¬ 
ing  two  exceptions,  all  actuators  and  sensors  are  the  same  ones  used  in 
the  active  control  subsystem. 

1)  Two  additional  linear  momentum  exchange  devices  are  mounted 

at  the  mid-mast  in  the  XQ  and  Y0  directions.  These  two  actuators 
are  used  to  excite  the  second-order  mast  bending  modes  (22.4 
and  23.0  Hz),  which  are  critical  in  distinguishing  the  masses 
at  the  secondary  mirror  and  mid-mast. 

2)  Accelerometers  and  gyros  are  mounted  at  the  base  connecting 
the  structure  to  the  Orbiter.  When  this  acceleration  infor¬ 
mation  is  available,  one  can  isolate  the  test  structure  dynamically 
from  the  Orbiter  motions,  permitting  the  use  of  the  substructure 
identification  technique  described  earlier. 

6-4-2-3  Physical  Parameter  Estimation  Results 

Five  structural  parameters  are  determined  as  critical  physical  parameters 
that  must  be  estimated.  They  are: 

^  E  Young's  modulus  of  the  graphy  epoxy  material  used  to 

make  the  mast  and  most  of  the  primary  mirror  support 
structure 
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Figure  6-2.  Simulation  Software  for  Physical  Parameter 
Estimation 
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Figure  6-3.  Actuator  and  Sensor  Placement 
for  Parameter  Estimation 
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The  next  step  is  to  select  a  test  signal  that  will  make  the  unknown 
parameters  more  identifiable  from  structural  responses.  Table  6-5  lists 
the  actuators,  and  the  major  structural  modes  and  related  physical  parameters 
they  make  more  identifiable. 

There  are  procedures  for  computing  the  "optimal"  test  signal  which 
maximizes  the  Fisher  information  matrix  (high  estimation  accuracy),  while 
satisfying  some  energy  constraint.  But  optimal  input  design,  in  general, 
is  a  difficult  problem  to  solve.  The  results  depend  on  the  unknown  parameters 
and  may  violate  the  test  structure  load  constraint  or  linearity  assumption. 

The  approach  taken  in  this  study  is  to  design  a  test  signal  which 
produces  reasonably  good  estimates.  Several  test  signals  were  tried.  It 
was  found  that  the  masses  of  the  secondary  mirror  and  the  mid-mast  are 
difficult  to  distinguish  if  the  second  mast  bending  modes  (modes  5  and  6) 
are  not  sufficiently  excited.  Also,  the  Poisson  ratio  and  the  mass  of  the 
primary  mirror  are  highly  correlated  if  the  primary  mirror  torsional  mode 
(mode  4)  is  not  excited.  A  test  signal  that  can  excite  modes  1-6  is, 
therefore,  crucial.  An  example  of  such  signals  is  shown  in  Figure  6-4, 

Structural  responses  are  measured  using  the  SAMS  sensor  (Section  7), 
which  measures  the  relative  range  from  the  mid-mast  (focal  point),  to  nodes 
8,  14,  and  15  of  the  primary  mirror;  these  measurements  are  denoted  as 
Yl>  y2>  and  y3>  respectively.  Gaussian  white  noise  of  1  pm  (la)  is  added 
to  the  measurement  signal.  The  accelerations  at  the  node  connecting  to  the 
Orbiter  are  also  measured;  this  permits  one  to  Identify  the  test  structure 
parameters  in  spite  of  Orbiter  structural  vibrations  and  VRCS  firings. 


The  physical  parameter  estimation  technique  is  applied  to  1  second 
of  simulated  data.  Figure  6-5  shows  the  comparison  of  the  measured  data 


Actuators 


Major  Modes 


Physical  Parameters 


Secondary  Mirror 

2,  3 

E,  M 

x2,  y2 

’  top 

Primary  Mirror 

1,  4 

E,  v,  M  . 

V  h 

prim 

MEASURED  DATA 
PREDICTED  DATA 


o  *1 

Y1  ~  I- 

(ICTER) 

— /V 

0 

o 

i2  3 

(METER)  ~  °  ' 

0  0  1  0  *  L  •  Of  •  o 

T«  'UCO.JI 

• 

f*  », 

o 

Y3  3 

(hCTER) 

9  o  /  0  «  Oft  6. a  t  o 

»tn  umoi 

0  1 

sO 

a 

O 

K  ft- 

residual" 

Y1  -Iff 

(*TER)  M 

>  tf  1  0*  P  ft  Of  If 

T|«t  lUCMIt 

^  a  c 

2., 
RESIDUAL—  ; 

Yz 

(METER)  [r 

1  0.1  0  «  Of  0,0  1.0 

'IK  iStCOBOi 

o.c 

so 

2*1 

RESIDUAL— 

y3  • yi 

(METER)  ’ 

-*■  — 

1  0  1  0  •  0.4  0.0  it 

'IK  'tlDMl 

l^y 

•  o.i  g«  ss  o.i  i ,  o 

TltC  IIIOMI 


Figure  6-5.  Comparison  of  the  Measured  Response 
with  the  Predicted  Response  for 
Physical  Parameter  Estimation 


with  the  data  predicted  using  the  final  estimated  parameters.  The  differences 
are  so  small  one  can  hardly  tell  them  apart.  The  residuals  which  are  the 
differences  between  measured  signal  and  predicted  signal  are  also  shown 
in  Figure  6-5.  The  estimated  parameters  in  each  of  the  Newton-Raphson 
iterations  are  shown  in  Figure  6-6. 

The  shape  of  the  cost  function  which  must  be  minimized  to  yield  the 
maximum  likelihood  estimates  depends  on  both  the  input  signals  and  the 
physical  parameters.  For  the  input  signal  presented  before,  the  cost  function 
for  some  of  the  physical  parameters  are  shown  in  Figure  6-7.  In  each  case 
only  one  parameter  is  allowed  to  vary,  all  other  parameters  ar.e  fixed  at 
their  true  values.  As  we  h"ve  studied  before  in  ACOSS-8,  the  cost  function 
of  longer  data  (2  seconds)  has  a  sharper  curve  at  the  minimum  point  (which 
implies  greater  estimation  accuracies)  than  that  of  shorter  data  (1  second). 
However,  the  region  of  convergence  for  2  seconds  of  data  is  still  very 
wide,  so  one  can  start  the  estimation  algorithm  with  poor  initial  estimates 
and  it  will  still  converge. 

6*4. 2.4  Results  of  Modal  Parameter  Estimation 

The  modal  parameters  of  the  first  four  structural  modes  are  obtained 
using  3  separate  tests.  In  each  test  only  one  actuator  is  energized  and 
the  corresponding  critical  modes  are  estimated. 

The  time  histories  of  data  involved  in  a  typical  modal  parameter 
estimation  is  shown  in  Figure  6-8.  The  input  signal  is  a  rectangular  pulse 
with  an  amplitude  of  1  Newton  and  duration  of  0.3  second.  The  measured 
output  is  the  relative  rate  of  deflection  of  the  secondary  mirror  in  the 
x  direction.  The  output  measurement  contains  the  responses  due  to  Orbiter 
VRCS  firings  and  the  dynamic  interaction  with  Orbiter  structural  modes. 

Also,  for  the  particular  case  shown,  the  output  measurement  has  nonzero 
initial  dynamic  conditions  which  have  to  be  estimated  along  with  other 
modal  parameters.  The  predicted  signals  are  computed  using  the  resulting 
parameters  after  5  Newton-Raphson  iterations.  Both  the  input  and  output 
signals  are  sampled  at  50  Hz  and  quantized,  with  proper  scaling,  to  8  bits. 
Since  it  is  planned  that  the  parameter  estimation  be  done  on  ground,  the 
predicted  signal  is  generated  by  a  ground  computer  with  a  much  higher 
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quantization  level.  The  impact  of  sample  rate  and  quantization  level 
on  parameter  estimation  accuracy  will  be  discussed  later. 


The  results  of  the  three  modal  parameter  tests  are  summarized  in 
Table  6-6.  In  some  cases  the  initial  dynamic  conditions  are  also  estimated. 
Since  the  initial  condition  is  a  function  of  the  particular  data  starting 
point  and  not  properties  of  the  test  structure,  they  are  not  shown. 

6. 4. 2. 5  Results  of  Frequency  Response  Method 

Figure  6-9  shows  the  signals  involved  in  estimating  the  second  mode 
of  the  test  structure.  Band-limited  (25  Hz)  Gaussian  random  noise  (a)  is 
applied  to  the  actuator  at  the  secondary  mirror  in  the  X  direction  and  the 
structural  response  at  the  same  location  is  measured  (b).  The  frequency 
contents  in  the  first  20.48  second  data  of  (a)  and  (b)  are  computed  using 
the  Fast  Fourier  Transforms  (c  and  d).  The  frequency  response  curve  (e) 
is  obtained  from  the  ratio  of  curves  (d)  and  (c). 

A  computer  program  is  developed  to  estimate  the  modal  frequency, 
transfer  function  gain  and  damping  ratio  by  fitting  a  second-order  transfer 
function  to  the  frequency  response  curve  in  the  neighborhood  of  each  modal 
spikes.  Figure  6-10  shows  the  result  of  this  frequency  domain  matching  for 
the  second  mode. 

The  estimation  results,  summarized  in  Table  6-7,  show  very  accurate 
modal  frequency  estimates  but  rather  poor  estimates  of  modal  damping  ratios 
and  transfer  function  gains.  These  results  are  expected  from  our  ACOSS-8 
work  and  this  is  one  reason  for  using  this  approach  only  as  a  cross-check. 

6. 4. 2. 6  Quantization  and  Sample  Rate 

Tests  are  made  with  various  sample  rates  and  quantization  levels  for 
the  modal  parameter  estimation.  Figure  6-11  shows  the  average  (over  unknown 
parameters)  parameter  error  as  function  of  sample  rate  for  the  cases  of 
8-bit  and  16-bit  fixed  point  quantization.  The  quantization  level  is  selected 
so  that  all  signal  amplitudes  are  within  its  dynamic  range.  Both  8-bit  and 
16-bit  quantization  produce  almost  identical  results  in  estimation  accuracy 
which,  however,  degrades  significantly  for  sample  rate  below  50  Hz.  Thus, 
the  baseline  quantization  level  is  specified  at  14  bits  and  the  sampling 
rate  at  175  Hz. 
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7.0  DATA  ACQUISITION  AND  INSTRUMENTATION 


W.th  the  establishment  of  the  baseline  concept  and  its  basic 
compatibility  with  the  Shuttle,  finer  details  of  the  experiment  system 
can  be  defined  and  analyzed.  This  section  defines  how  experimental  data 
is  to  be  acquired,  stored,  and  transmitted.  The  system's  electrical  and 
thermal  requirements  are  also  briefly  assessed  and  compared  with  the 
capability  of  the  Shuttle. 

7.1  LOS  Measurement  and  Optical  Alignment 

The  performance  criteria  for  the  technology  demonstrations  will  be 
optical  element  motion  and  LOS  motion/image  quality.  A  concept  for 
measuring  it  is  illustrated  in  Figure  7-1.  The  basic  idea  to  measuring 
optical  element  motion  is  to  deploy  8  Ga  A1  As  laser  diode  (LED) 
transmitters  on  the  perimeter  of  the  primary  mirror,  preferably  at 
rigid  points  subtending  from  the  support.  The  six  degree-of-freedom 
motion  of  these  points  (hence  primary  mirror  rigid  body  motion)  is 
then  directly  measured  with  a  TRW  Surface  Accuracy  Measurement 
Sensor  (SAMS)  rigidly  attached  to  the  focal  plane.  The  motion  of  the 
secondary  mirror  is  indirectly  measured  by  observing  the  LED  reflection 
off  the  secondary  with  CCD  focal  plane  sensor  and  then  subtracting  the 
LED  image  motion  caused  by  the  primary.  In  order  to  make  the  LED 
light  focus  on  the  foca'  olane,  four  of  the  LEDs  are  lensed  to  focus 
at  F. 


The  LOS  and  imaging  performance  of  the  system  will  be  measured 
with  the  same  charge-coupled-device  (CCD)  sensor  system  above,  one 
similar  to  TRW's  MAuAN.  LOS  is  directly  measured  by  el ectronically 
tracking  the  centroid  of  a  stellar  or  other  point  source  truth  target 
image  in  the  sensor.  Imaging  quality  is  deduced  by  observing  and 
el ectronical ly  imaging  a  more  complicated  truth  target  and  ground. 

7.1.1  CCD  Sensor  Requirements 

The  LOS,  imaging  and  secondary  mirror  range  sensing  functions 
defined  above  are  performed  by  a  CCD  array.  TRW  is  currently  developing 
such  a  sensor  (MADAN)  for  other  applications,  a  picture  of  it  is  shown 


in  Figure  7-2.  During  the  course  of  the  study,  the  MADAN  system 
designer  was  consulted  to  outline  the  requirements  and  modification 
required  to  make  this  sensor  work  for  the  present  application.  The 
optical  requirement  is  to  have  a  .Olx.Ol  rad  field  of  view  with  approx¬ 
imately  .2  yrad  resolution.  The  throughput  requirement  is  to  be  able 
to  sample  "figure"  at  175  Hz,  and  sample  "image"  at  two  frames/second. 

The  modification  required  to  make  the  MADAN  sensor  meet  these 
specifications  are  to  replace  the  current  quad  arrangement  of  Hughes 
CCD  arrays  (see  Figure  7-3)  by  four  RCA53612  arrays.  The  RCA  chips  are 
540x512  pixels  each,  so  if  one  optically  shutters  the  set  and  uses  the 
standard  .05  pixel  interpolation  algorithm, .2  yrad  resolution  is 
possible.  One  can  further  improve  this  resolution  by  using  the 
preferred,  .02  pixel  interpolation  algorithm.  The  resolution  with  this 
algorithm  is  .1  yrad. 

In  order  to  meet  the  throughput  requirements,  it  is  assumed  that 
the  laser  diodes  (designating  critical  structure  locations)  are  aimed 
at  pre-speci fied,  nominal  locations  on  the  sensor.  The  area  surrounding 
these  locations  are  used  only  for  structure  measurements  and  not  to 
image.  The  electronics  then  channels  the  pixel  samples  to  either  the 
image  or  structural  measurement  processor,  depending  on  the  location 
of  the  pixel  on  the  chip.  Pixel  windowing  is  specified  at  10x10  with 
all  four  chips  running  iri  parallel  at  .1  ysec/pixel.  The  windowed  data 
i;;  accumulated  using  .0015  seconds  integration  and  then  line  dumped 
at  .4  ysec/line.  The  signal  strength  provided  by  the  laser  diodes  has 
been  analyzed  and  it  will  permit  175  Hz  samples  of  the  structure  with 
good  S/N.  Similarly,  the  light  gathering  capability  of  the  telescope 
has  been  analyzed  and  it  will  permit  two  frame/second  imaging  rates. 

The  power  consumption  of  the  CCD  sensor  is  conservatively  estimated 
at  31  W,  including  data  processing  electronics.  The  total  weight  is 
less  than  20  lbs. 


7.1.2  Surface  Accuracy  Measurement  Sensor  (SAMS)  Requirements 


A  modified  SAMS  has  been  baselined  for  ranging  the  primary  mirror. 
The  SAMS,  shown  in  Figure  7-4  was  originally  developed  by  TRW  for  NASA 
Langley  and  was  delivered  to  them  in  1980.  Since  then,  the  sensor  has 
continually  been  improved  through  better  detectors,  electronics,  and 

element  stabilization.  The  unique  properties  of  the  sensor  are:  its 

4  /  \ 

5x10  dynamic  range  (independent  of  field  of  view);  its  excellent 

linearity  properties;  and  bandwidth  greater  than  1  KHz.  The  sensor 
plus  electronics  consumes  about  20W,  and  the  space  qualified  sensor 
would  weight  10  lbs. 

For  the  present  application  the  field  of  view  of  the  sensor  would 
be  designed  at  5000  yrad.  This  would  allow  a  .1  yrad  measurement 
accuracy.  In  order  to  simultaneously  view  all  of.  the  laser  diodes  on 
the  primary  mirror  support,  the  optics  on  the  sensor  head  will 
spacially  multiplex  the  nominal  range  angles  into  the  FOV.  This  will 
allow  up  to  +400  yrad  excursion  at  each  range  point.  The  approach  has 
been  successfully  used  elsewhere. 

The  SAMS  sensor  can  only  measure  angles  relative  to  its  own 
coordinate  frame.  In  order  for  it  to  range  objects  relative  to  itself, 
some  additional  information  must  be  available.  One  approach  to  obtain 
this  information  is  through  the  stadiametric  ranging  technique  described 
in  Figure  7-5.  The  approach  is  described  for  a  planar  coordinate 
frame  for  simplicity.  The  basic  idea  is  that  if  one  knows  the  relative 
distance  between  at  least  three  points  (four  points  for  three  dimensions) 
on  an  object,  and,  one  knows  the  angles  to  the  targets  on  the  object, 
then  the  object  can  only  lay  in  one  of  two  positions.  These  are  marked 
with  solid  and  dashed  lines  in  the  figure.  But  note  however,  that  in 
on  case  the  angle  to  C  is  greater  than  the  angle  to  B,  and,  in  the 
other  case  the  angle  to  B'  is  greater  than  to  C1.  So  if  one  can 
associate  the  measured  angle  to  the  targets  with  the  targets,  the 
ambiguity  is  resolved,  and  the  range  to  all  three  targets  can  be 
computed.  The  concept  used  in  the  SAMS  is  to  uniquely  tag  the  target 
by  pulsing  them  at  different  rates. 
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7.1.3  Simulated  Performance  of  the  SAMS  and  CCD  Sensor 

Monte  Carlo  runs  were  used  to  determine  the  LOS  and  figure  sensing 
capabilities  of  the  baselined  SAMS  and  CCD  sensor  systems.  The  simula¬ 
tion  results  of  the  SAMS  showed  that  it  will  be  able  to  stadiametrical ly 
range  the  position  of  point  sources  in  the  primary  mirror  to  0.16  pm  (3a) 
in  all  three  axes.  This  result  includes  various  expected  random  mis¬ 
alignments  expected  by  the  SAMS  designer:  referring  to  Figure  7-5,  .1  pm, 
la  uncertainty  in  the  reference  distance  S  and  R;  .1  prad,  la  uncertainty 
in  reference  angle  0;  and  .1  prad,  la  measurement  error  in  a,  6  and  y  . 
This  information  was  incorporated  with  the  expected  capability  of  the 
baselined  CCD  sensor  in  order  to  estimate  the  ranging  sensitivity  to 
the  secondary  mirror.  This  resulted  in  0.2  pm(3a),  3-axis  ranging. 
Integrating  this  information  into  the  baselined  optical  train  sensitivity 
further  resulted  in  a  0.8  prad  closed-loop  vibration-control  reference. 

The  CCD  sensor  itself  is  expected  to  resolve  point-source  LOS 
measurements  to  0.1/0. 2  prad  for  0.02/0.05  pixel  interpolation,  respec¬ 
tively.  Hence,  the  CCD  sensor  will  provide  a  closed-loop  target  tracking 
reference  which  will  easily  be  able  to  measure  the  performance  of  the 
vibration  control  system, 

7.2  Experiment  Control  and  Data  Handling  System  (CDHS) 


All  of  the  electrical,  mechanical  and  optical  hardware  described 
(and  to  be  described)  interacts  and  is  baselined  to  be  controlled  by  the 
control  and  data  handling  system.  Referring  to  Figure  7-6,  the  baseline 
architecture  for  this  system  is  "top-down"  with  a  special  provision  made 
to  automate  the  experiment:  Man  interacts  with  the  experiment  through 
the  Experiment  Monitor  Software  package.  Data  entry  is  made  by  the 
astronaut  through  a  data  terminal  (CRT,  keyboard,  etc.)  or  by  ground 
through  a  radio  link.  The  experiment  monitor  will  enable  either  a 
manually  controlled  experiment  or  a  complete  set  of  pre-defined  experi¬ 
ments.  Manual  control  is  allowed  through  a  data  bus  to  all  experiment 
subsystems.  Automated  control  is  executed  through  an  Automated  Experi¬ 
ment  Sequencer  and  Supervisor  package.  Whether  manual  or  automatic, 
the  experiments  defined  earlier  (system  tests,  technology  demonstrations 


rHf~U  SCF^ARE 
:«0 


WglMT* 

eiperi^wt  kjhitoh  inti^aci  n  *aim*  aw  status  graphics 
EIPERIXNT  >O*»T0 *  INTERFACE  TD  SWTTL£  (DAP) 

AUTO*TlD  EAPERI^IT  SEQUENCER  ENABLE 

upeb:»€*t  T-ams  hjhitoh 

i  1DEQ  CONSOL 

rectwoeb  coktroi 

EXPER  1>CMT  HAXUAl  CONTROL 
FAILURE  DETECTION  EHABLE 
HAS  TIB  PC*ER  EHABLE 

WHM. 

ajC*aonics 

*€  CHAM  CAL 

HANUAL  ElFLRl^WT  SELECT 
Stsri"  TST 

TtC.-MLOGT  XWHSTRATIOH 
MSCElLABEOuS 

EJ^Rl^.HT  TVL^ICC/REPROGflAX 
S,r?<5!-E¥T  ’APVfTR  MOOIPf 


WHrTOH  HAAOUAAE 


CRT 

STB  IP  CHAATS 

tOBOMD 

GRAPHICS 


^ _ 

GOTH®  SOFTWARE 

DATA  PROCESSING 

THERMAL 
ELECTRONICS 
*€C:<AHICAl 
TAPE  i  VIDEO 
3I5TURBANCE  G£B. 


Figure  7-6.  Baseline  Control  and  Data  Handling  System  Architecture 
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and  miscellaneous)  are  stepped  through  by  commanding  the  state  of  the 
hardware  via  a  Hardware  Enable  and  Supervision  System.  The  purpose  of 
the  supervision  system  is  to  flag  hardware  anomalies  to  the  experiment 
monitor  so  that  corrective  action  can  be  taken.  As  the  hardware  is 
stepped  through,  the  Control  and  Identification  Software  will  execute 
according  to  manual  or  automatic  algorithms  stored  in  memory  or  in  the 
Experiment  Sequencer.  Software  failure  detection  with  an  overridable 
automatic  shutoff  is  also  provided  here  as  an  additional  precaution. 

During  the  execution  of  the  experiments,  we  have  identified  31 
signals  to  be  monitored  and  stored  for  ground  processing.  These  include 
six  disturbance/stimuli  generators,  four  active  control  sensors,  twelve 
optical  autoalignment  sensors,  six  inertial  sensors,  three  digital- 
autopilot  signals,  and  time.  In  addition  to  these,  we  have  added 
another  nine,  as  a  pad,  for  a  total  of  40.  The  experiments  command  and 
data  handling  system  (CDHS)  will  multiplex  these  measurements  and 
communicate  the  data  to  the  Shuttle  payload  tape  recorder.  The  tape 
recorder  is  part  of  the  Shuttle  standard  equipment,  and  will  store  all 
the  information  gathered  during  the  duration  of  the  experiment.  After 
the  experiments  are  completed,  the  data  will  be  downlinked  on  the 
Shuttle's  Ku-band  channel.  The  signal  is  received  by  standard  NASA 
links  where  it  is  then  recorded  on  the  ground,  ready  for  processing. 

The  data  recording  requirements  were  computed  to  be  128K  bits/second 
This  assumes  forty  data  signals  sampled  at  200  Hz,  each  16  bits  long. 

The  total  data  record  time  is  upper  bounded  at  111  minutes,  including 
various  crosschecks  and  safety  factors. 

Two  data  recorder  modes  meeting  the  above  specification  were 
considered:  6  inches/second  with  fourteen  parallel  tracks  recording  a 
three-to-one  multiplexed  signal;  and  30  inches/second  with  a  single 
track  recording  a  forty- to-one  multiplexed  signal.  We  are  baselining 
the  second  mode  as  it  appears  to  be  simpler  overall,  and,  it  will  be 
able  to  store  data  for  224  minutes  at  the  required  throughput  rate. 

The  downlink  operation  will  be  automated  and  enabled  by  the  CDHS. 

It  will  replay  the  recorded  data  at  120  inches/second  with  all  fourteen 


Table  7-3.  Summary  of  Electrical  and  Thermal  Requirements 


i 

o 

E 

u 

CD 

jC 

S- 

“O 

"O 

a> 

CD 

a) 

i — 

c 

c 

o 

•r 

•r— 

o 

03 

m 

u 

4-> 

1  1  1  1  1 

c 

c: 

a  i  i  i  i  i 

o 

o 

’r 

CJ 

<j 

S- 

4-> 

4— 

4- 

u 

t— 

r— * 

0} 

CD 

CD 

m 

c/i 

cd 

c/l 

c 

CD 

s- 

E 

CL> 

CD 

5 

S- 

O 

•r— 

Q- 

3 

cr 

(D 

3 

3 

3 

3 

3 

3 

O 

rH 

cn 

LD 

CO 

00 

C\J 

ro 

r— 

CNJ 

CXI 

ro 

CNJ 

r-H 

CO 

r-4 

E 

CD 

CD 

E 

> 

4-> 

CD 

•r— 

E 

l/l 

4-> 

4-> 

CD 

>> 

C/1 

Ql 

4-> 

CO 

>» 

03 

C/1 

CO 

“O 

E 

>> 

CD 

cC  E 

CD 

in 

C 

c 

CD 

-M 

•r- 

o 

C/l  4-> 

C/l 

s- 

j— 

•r— 

3  01 

>> 

E 

o 

■a 

4-> 

•—  >» 

CO 

CD 

■M 

c 

03 

Q-  CO 

4-> 

•r“ 

ra 

E 

CD 

C/1 

c 

31 

E 

•r“ 

CD  »— 

c 

>> 

o 

CD 

4-> 

S_  O 

•r“ 

CO 

2: 

03 

C/1 

3  S- 

(/) 

4-> 

C/1 

LU 

4->  +■> 

c 

S- 

03 

>1 

a  C 

CD 

o 

c 

Q 

CO 

L- 

3  O 

m 

c/l 

CD 

CD 

S-  o 

c 

£= 

“O 

r— 

4-> 

4-> 

CD 

CD 

•r“ 

c 

o 

cd 

CO  C/l 

S- 

CO 

u 

03 

S- 

E 

u 

3 

CD 

E 

4— > 

03 

4—1  -r- 

CD 

Q 

CL 

E 

c 

S- 

C/l  4-> 

•r“ 

O 

X 

O 

o 

03 

CD  CL 

U. 

O 

LU 

o 

o 

Q- 

h-  O 

8.0  DEVELOPMENT  OF  THEORY  FOR  NONCOLOCATED  CONTROL 

As  a  result  of  an  earlier  study  performed  by  TRW  (ACOSS-8),  a  stability 
ensuring  design  methodology  was  developed  for  LSS  systems  permitting  the 
colocation  of  actuators  and  sensors.  The  design  procedure  has  the  advantage 
that  it  is  relatively  insensitive  to  modal  parameters.  It  permits  design 
objectives  to  be  designed  for  and  sought  in  both  the  time  and  frequency 
domains.  And,  the  approach  leads  directly  to  solution  that  can  be  imple¬ 
mented  using  current  microprocessor  technology.  In  a  number  of  applications 
of  interest  however,  colocation  is  not  possible  and  further  theoretical 
development  was  required. 

One  objective  of  the  current  MOD  contract  study  was  to  extend  the 
ACOSS-8  theory  to  include  noncolocated  LSS  systems.  This  section  describes 
the  approach  taken  to  meet  this  objective  and  the  resulting  design  extension. 
The  plan  was  to  re-examine  the  ACOSS-8  design  approach,  determine  the  proper¬ 
ties  that  make  it  work,  explore  ways  of  synthesizing  the  desired  properties, 
and  derive  the  extended  design  procedure. 

8. 1  Review  of  ACOSS-8  Results 

8.1.1  Stability  Ensuring  Design  Method 

The  stability  ensuring  design  methodology  in  ACOSS-8  makes  use  of 
positivity/energy  dissipation  theory  to  ensure  the  stability  of  the  LSS. 

The  key  results  were: 

*  The  negative  feedback  connection  of  a  positive  system  with  a 
strictly  positive  system  is  stable 

•  An  LSS  with  collocated  actuators  and  ideal  rate  sensors  is  a 
positive  system 

e  Embedding  operations  permit  the  treatment  of  non-positive  systems 

The  power  of  the  theory  becomes  apparent  as  one  realizes  that  no  assumptions 
are  made  about  the  dimensionality  or  linearity  of  the  systems  involved.  It 
is  for  this  reason  that  the  theory  permits  the  treatment  of  modal  truncation, 
stability  robustness  and  nonlinearities.  The  price  one  pays  for  the  generality 
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is  however,  that  designs  tend  to  be  conservative.  It  was  therefore  proposed 
in  ACOSS-8  to  use  positivity  in  conjunction  with  on-orbit  identification/ 
control  "tune-up":  Positivity  is  used  in  the  pre-flight  design  where  there 
is  a  high  parameter  uncertainty,  then  as  information  is  gathered  on  orbit, 
multivariable  frequency  domain  techniques  combined  with  positivity  (to 
treat  non-linearities)  is  used  to  tune  the  design  for  high  performance. 

It  was  clear  that  as  system  knowledge  increases,  the  need  for  a  conser¬ 
vative  pre-flight  design  diminished.  In  cases  where  good  pre-flight  informa 
tion  is  available  over  ^  twice  the  control  band,  the  posi ti vi ty/ frequency 
domain  approach  is  to  be  applied  directly  without  on-orbit  identification. 

8.1.2  Summary  of  Positivity  Theory 


ons  are  omy  sufficient  conditions  tor  staointy 
nservative  at  times.  Since  well-known  necessary 
ist  to  establish  the  stability  of  linear  systems 
y  techniques  have  in  the  past  generally  not  been 


M 


SB 


In  order  to  summarize  the  theory,  some  definitions  are  required, 

Let  f(t)  be  a  real,  possibly  vector  valued  function  defined  on 
[0,«).  The  truncation  of  f(t)  at  t=T  is  defined  as 


fT(t) 


f(t),  0  -  t  -  T 

0  ,  t  >  T 


Define  the  extended  Hilbert  Space  Hg  such  that  f(t)  e  H  if  and  only 
if  f  (t)  e  Hilbert  Space,  H. 

An  operator  with  a  domain  and  range  in  Hg  is  then  defined  to  be 
(strictly)  positive  if  for  f(t)  e  H  , 

<fT(t),  HfT(t)>  -  0  (>0) ,  VT  -  0 

Positive  operators  have  the  following  useful  properties: 

1)  The  inverse  of  a  positive  operator  is  positive. 

2)  The  sum  of  positive  operators  is  positive. 

3)  Positivity  Theorem  [11,  12,  13]:  The  negative  feedback  system 
shown  in  Figure  8-1  is  bounded-input-bounded-output  BIBO 

stable  for  all  inputs  in  H  if  both  G  and  H  are  positive  operators 
and  at  least  one  of  them  is  also  strictly  positive. 


Figure  8-1.  Cascade/Unity  Feedback  Interconnection 


3.1.3  Operator  Embeddinc 


The  condition  imposed  by  the  Positivity  Theorem  are  often  too  restricti 
in  practice.  For  example  when  actuator  and  sensor  dynamics  are  accounted 
for,  or,  when  actuators  and  sensors  are  not  colocated  the  LSS  will  no 
longer  be  positive.  Therefore,  the  technique  of  operator  embedding  was 
introduced  in  ACOSS-8  in  order  to  enlarge  the  class  of  plants  that  can 
be  treated  by  the  theorem. 


Embedding  transformations  were  defined  in  ACOSS-8  as  the  set  of  block 
diagram  manipulations  which  do  not  alter  a  system's  stability  property. 

Two  types  of  embedding  transformations  were  considered:  "F"  embedding 
refers  to  the  cascade  transformation  in  Figure  8-2. 


Figure  8-2.  "F"  Embedding 

D"  embedding  refers  to  the  parallel  transformation  in  Figure  8-3 


Figure  8-3.  "D"  Embedding 

The  important  thing  to  realize  about  embedding  transformations  is  they  ar 
only  mathematical  and  are  not  physically  implemented.  The  embedding  may 
impose  however,  conditions  on  the  controller  which  is  to  be  physically 
implemented.  For  example,  constant  "D"  embedding  may  impose  that  the  gain 
of  the  controller  must  be  less  than  some  number  at  all  frequencies. 


4/f’jCv 


Generally,  embedding  requires  some  knowledge  of  the  plant  to  be 
controlled.  For  example,  constant  "D"  embedding  will  require  that  one  has 
estimated  (albeit,  crudely)  the  norm  of  the  mode  shapes  and  damping  ratio 
of  the  high  frequency  modes.  If  one  has  additional  knowledge  of  the  plant, 
more  complex  embedding  can  be  used  to  yield  less  conservative  stability 
conditions.  In  a  gross  manner,  embedding  permits  one  to  quantify  stability 
conditions  as  a  function  of  system  knowledge. 

8.1.4  Linear,  Multivariable  Frequency  Domain  Techniques 

The  linear  multivariable  frequency  domain  design  approach  was  used  in 
ACOSS-8  to  guide  the  control  system  design  whenever  reasonably  accurate 
system  models  are  available.  The  combination  of  it  with  Positivity  allows 
one  to  trade  sensitivity  properties  for  performance  and  vice-versa.  The 
combinations  will  be  summarized  in  Section  8.1.5. 

The  key  to  the  extension  of  the  classical  Nyquist/Bode  approach  to 
the  multivariable  case  are  the  characteristic  gains  of  transfer  function, 
G(S) ,  describing  the  open-loop  plant.  The  characteristic  gain  X ( S)  are 
defined  as  the  solutions  to  the  equation 

DET  [X (S)  ’  I  -  G(S )]  =  0;  Sc  [oj,»j) 


This  equation  arises  from  the  eigenvalue  problem  for  a  transfer  matrix 


G'SW 


G(S)  W.j  (S)  =  X.(S)  W.  (S)  is  evaluated  for  each  Se  Cbj^j) 


where 


A .  ( S )  are  called  the  characteristic  gains  of  G(S) 


W.  (S)  are  called  the  characteristic  vectors  of  G(S) 


V.j(S)  are  called  the  dual  characteristic  vectors  of  G(S) 


The  importance  of  the  characteristic  gains  stem  from  the  fact  that  the 
closed-loop  characteristic  gains  of  a  unity  feedback  system  are  related  to 


the  open-loop  values  in  a  manner  analogous  to  single-input-single-output 
sys  terns .  i .  e . , 


closed-loop  gain  =  open-loop  gain 

*  4<  n  n  _  nn  n  /n 


1  +  open-loop  gain 


This  fact  can  be  demonstrated  by  considering  the  dyadic  expansion  of  G(S) 

m  T 

G(S)  =  z  X.  (S)  W. (S)  V  (S)  =  W(S)  [X.(S)]  V(S) 
i=l  1  1  i  1 


If  one  now  closes  the  loop  around  G  (shown  in  Figure  8-4) 


IB 

m 

a 

Figure  8-4.  Unity  Feedback  System 

and  considers  the  dyadic  expansion  of  the  closed  loop  system  Gc>  then 

GC(S)  =  [i+g(S)J'1G(S)  =  [I+W(S)  X(S)  V(S)]'1  W(S)  X(S)  V(S) 

=  [W(S)  V(.S)  +  W(S)  X(S)  V(S)]'1W(S)  x (S)  V (S ) 

=  W(S)  [I+x(S)]_1  X (s )  v(s) 
m  Xi/-,  T 

-  ^  TTxfsT  w-(s)  v  (s) 

i=l  1  |  i 

Other  key  results  using  this  technique  are: 

1)  Generalized  Nyquist  Criterion:  Closed-loop  stability  occurs  if 
and  only  if  the  net  sum  of  counter-clockwise  encirclements  of  the 
(-1,  0.)  by  the  characteristic  gains  is  equal  to  the  number  of 
open-loop  unstable  poles. 

2)  To  achieve  low  control  interaction  in  a  frequency  band,  it  is 
sufficient  that  either  all  the  characteristic  gains  have  a  large 
magnitude  over  that  frequency  band,  or  the  dual  eigenvectors  are 
nearly  orthogonal  over  that  band. 


3)  In  frequency  bands  of  low  interaction,  the  notions  of  gain  margin 
and  phase  margin  may  be  applied  as  a  qualitative  assessment  of 
performance. 

4)  There  is  a  relationship  between  tracking  accuracy  and  the  magnitude 
of  the  characteristic  gains.  Under  the  appropriate  assumptions 

the  relationship  is  analogous  to  the  single-input-single-output 
case. 


8.1.5  Positivity  Design  Steps 


In  practice,  both  positivity  and  frequency  domain  compensation*  may  be 
needed  in  order  to  meet  stability  and  performance  objectives.  The  two  are 
combined  by  first  determining  the  characteri stic  gains  and  positivity  index 
of  the  plant  and  shaping  them  using  techniques  similar  to  those  of 
Mac  Farlane  [15].  The  shaping  network  is  a  precompensator  M(s)  which  is 
designed  for  performance.  Once  the  design  of  M(s)  is  completed,  one  can 
determine  the  D  embedding  operator  that  makes  the  precompensated  plant 
positive.  This  is  typically  done  by  determining  the  positivity  index  <S(o>) 
of  the  precompensated  plant,  finding  its  most  negative  peak  <5  ._,  and 


min 


choosing  D  =  Sq 1 ,  5q  >  “l5min  >  0.  The  "H"  part  of  the  controller  consists 


then  merely  of  gains  which  must  be  set  to  less  than  or  equal  to  1/<Sq. 

The  actually  implemented  controller  satisfying  the  positivity  design  is  then 


given  by  (1/Sq)M(s)  (See  Figure  8-5). 


PRE-COMPENSATED 
PLANT 


Figure  8-5.  "D"  Embedding  a  Precompensated  Plant 


*  ACOSS-8  also  describes  a  time  domain  alternative  using  optimal  control 
theory.  This  study  considers  only  the  frequency  domain  approach. 


The  actual  steps  are: 


1)  Determine  the  characteristic  gains  x(w)  and  positivity  index  fi(u) 
of  the  plant  +  actuators  +  sensors 

2)  Shape  characteristic  gains  and  positivity  index  (cascade  compensa 
tor  M[sJ)  using  Nyquist/Bode  approach  and'  CAD  filter  design 
program 

-  Design  for  damping 

-  Design  for  non  modal  sensitivity 

-  Design  for  noise  and  disturbance  rejection 

-  Design  high  frequency  roll-off  characteristic 

-  6 (to)  is  shaped  for  modal  sensitivity 

3)  Evaluate  positivity  index  s(ta)  for  entire  pre-cnmpensated  system 

and  determine  S  .  (maximum  negative  peak) 
mm  3  r 

4)  Implement  controller  as  ~  .  M(s)  where  6  >  -6  .  >  0 

o  o  min 

8.2  Generalizations  Reauired  for  Noncolocated  Control 


In  order  to  extend  the  ACOSS-8  approach,  the  distinguishing  features 
of  the  noncolocated  control  problems  were  first  isolated.  It  was  found 
that  the  unique  features  are  that  an  ideal,  noncolocated  LSS  plant  is  not 
necessarily  positive  nor  is  it  necessarily  square  (the  number  of  system 
inputs  and  outputs  are  not  necessarily  equal). 

The  first  possibility  was  considered  benign  with  respect  to  stability 
as  the  ACOSS-8  theory  was  never  restricted  to  positive  systems  at  all.  In 
theory  the  embedding  operation  required  to  make  the  embedded  system  positive 
can  always  be  found,  regardless  of  how  negative  the  plant  is.  This  was 
demonstrated  in  the  ACOSS-8  examples.  A  more  practical  problem  was  to  find 
however,  the  embedding  operations  which  optimize  performance  in  the  noncolo¬ 
cated  case.  Recall  that  embedding  can  be  interpreted  as  the  possible 
excursions  that  the  real  plant  can  take  from  its  nominal  model.  For  colocated 
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systems  both  the  real  plant  and  the  model  were  close  to  being  positive. 
Hence,  one  could  be  fairly  sloppy  in  choosing  "F"  embedding  without  much 
penalty  on  the  performance.  [Recall  that  embedding  imposes  constraints  on 
the  controller  gains,  and  therefore  affects  performance.]  Now,  however,  the 
nominal  and  real  systems  are  possibly  non-positive  and  the  role  of  the  "F" 
embedding  becomes  more  important. 


The  second  possibility  was  potentially  more  serious  as  both  the  posi¬ 
tivity  and  frequency  domain  theory  used  in  ACOSS-8  was  for  square  systems. 
Fortunately,  it  was  found  that  the  frequency  domain  theory  had  been  extended 
recently  [16,17],  but  it  was  not  easy  to  use  yet.  Moreover,  the  theory 
from  which  positivity  was  derived  did  not  require  the  system  to  be  square. 

The  later  finding  was  critical  indeed  as  it  is  what  ensures  robustness.. 

The  theoretical  effort  focused  into  three  areas:  generating  a  stability 
ensuring  design  condition  for  nonsquare  plants;  developing  a  procedure  for 
making  a  nonsquare  plant  square  (squaring  procedure);  and  developing  an 
erabedding/compensation  procedure  for  maximizing  performance. 

8.3  Generalizing  Positivity 

The  first  goal  is  to  develop  a  stabil ity-ensuring  design  criterion  for 
nonsquare  plants.  This  condition  had  to  be  simple  and  compatible  to  the 
conditions  derived  in  ACOSS-8.  The  strategy  to  meet  this  goal  was  to  revisit 
the  derivation  of  the  positivity  condition,  find  the  assumptions  requiring 
square  plants,  remove  them  and  then  continue  through  the  rest  of  the 
original  derivation. 

8.3.1  Conici ty 

It  was  found  that  the  Positivity  Theorem  is  a  special  case  of  the 
Conici ty  Theorem;  and  the  Conici ty  Theorem  can  be  derived  by  combining 
embedding  with  the  Small  Gain  Theorem.  Neither  the  Conici ty  or  Small  Gain 
Theorem  require  a  square  plant.  The  derivation  below  demonstrates  this 
and  follows  the  original  work  of  Zames  [13]. 


The  Small  Gain  Theorem  states  that  the  GH  system  in  Figure  8-6  is 
stable  if  for 


Figure  8-6.  Feedback  Interconnection 


g ( h)  =  sup 


Hf 

f 

,  V f S  Extended  Hilbert  Space 


g(H).g(G)<l.  Note  that  this  is  a  generalization  of  gain  stabilization.  If 
one  now  embeds  the  GH  system,  as  shown  in  Figure  8-7,  and  applies  the  small 
gain  condition,  it  is  implicit  that  if 


I | (G+D)  F- 1 | | <  1 


and 


-1, 


F  (I-H0)"'H|  |  <  1 


(8.1) 

(8.2) 


Then  the  GH  and  GH  systems  are  stable.  This  is  the  linear  version  of  the 
Conicity  Theorem.  The  name  comes  from  the  fact  that  the  inequalities  (8.1) 
and  (8.2)  imply  (if  the  inverses  exist) 


| | (G+D)  f | | 


and 


-D)  f 


IF  f||  Vf 

I  I F  f[ I  V  f 


*  The  Theorem  is  true  whether  or  not  the  loop  is  closed  with  H  in  the  feed¬ 
back  or  cascade  path. 


8.3.2  Sectors 


Since  Zames'  early  work  on  conicity,  many  new  theories  have  been 
proposed.  The  concept  of  sectors  [18]  is  to  our  knowledge  the  most  general 
input-output  stability  result  available.  It  includes  everything  described 
as  a  special  case.  The  embedding  operations  are  even  explicitly  accounted 
f o r  as  S-j-jj  S  -j  2 ,  S^i  and  S  2  2 1 

An  operator  G  is  said  to  be  inside/outside  sector 

S  S 

S  4  11  12 

s  s 

_  21  i22_ 

if  <Sn  Gf  +  S]2  f,  S2]  G  f  +S22  f  >  <0  (>  0 ) . 

The  pertinent  property  of  sectors  that  is  particularly  interesting  is  the 
Sector  Stability  Theorem  -  If  G  is  inside  a  sector  S  and  H"1  is  outside 
the  sector  S,  the  closed  loop  GH  system  is  stable.  Hence,  the  Small  Gain 
Theorem  resul ts  when 

S11  =  S21  =  S22  =  1  and  S1 2  =  ~l 
The  Conicity  Theorem  results  when 

S11  S21  =  I;  S12  =  D-F  and  S22  =  D+F 

The  Positivity  Theorem  results  when 

^11  -  ^22  -  0;  S-j2  =  -I  and  S22  =  I 

Sectors  also  possess  an  algebra  (summarized  in  Table  8-1)  that  makes  their 
manipulation  simpler  than  embedding.  This  makes  it  easy  to  derive  stability 
ensuring  design  conditions  on  a  case  by  case  basis.  Since  generating 
particular  results  is  not  the  intent  of  the  study,  conicity  (or  conic 
sector)  will  be  used  to  generate  the  general  results  sought.  Nevertheless, 
it  is  noted  that  the  sector  concept  may  be  required  in  practical  situations 
wherein  performance  cannot  be  sacrificed  for  an  overly  conservative  stability 
robust  design. 


Table  8-1.  Summary  of  Section  Properties  (Lemma  6.2  in  [18]) 


1.  Complimentary  Sector 
G  inside  sector 


S11  S12 
S21  S22 


G  outside  sector 


I  0 
0  -I 


2.  Multiplier 
sector 


S11  S12 
MS21  MS22 


sector 


aM*Sn  aM*S12 


3.  Inverse  Relation 


G  inside  sector 


4.  Sums  of  Relations 


S11  S12 
S21  S22 


S11  S12 


G”  inside  sector  S  • 


0  I 

I  0 

fl  G, 


G-|  inside  sector  s  s  4=^G^-G2  inside  sectors*  g  j 


5.  Composition  Products  of  Relation 
a.  fsn  S12 


G0Gi  inside  sector  c  c 
2  I  S21  ^22 


g2  0 

G-|  inside  sector  S*  g  j 


Go G-i  inside  sector  ^  <. 

L  1  b21  b2Z 


S11  S12 


I  0 

G2  inside  sector  S*  Q  g-1 


S11  G1  S12 

G0  inside  sector  c  r  c 

L  b21  bl  H2 


S11  S12 

G,  *G~  inside  <- 

<-  ^2i  ^22 


Notes:  a  b  >  0;  M*  is  the  adjoint  of  M;  S..  are  operators  mapping  into 

•  J 

extended  inner-product  spaces;  S- ■  0  =  0;  properties  1-4  hold  if 

*  J 

inside  is  replaced  by  outside. 


The  Conicity  Theorem  will  now  be  used  to  derive  a  stability-ensuring 
design  condition  for  noncolocated  systems. 

First,  the  uncertainty  of  the  plant  G  is  modeled  as 
1 1 G+D| j<  | | F| | 

This  is  the  first  condition  of  the  Conicity  Theorem  and  it  states  that  the 
true  system  G  differs  from  the  nominal  system  -D  by  at  most  F  (how  it  varies 
within  the  limit  F  is  arbitrary).  The  plant  G  is  then  partitioned  into 
[G-j ]  G2 ]  for  Case  1  (more  inputs  than  outputs,  Figure  8-8)  and  into  Gl_  for 
Case  2  (more  outputs  than  inputs.  Figure  8-9).  The  partitioning  L  2_ 
is  not  necessary  but  assigning  G^  to  represent  the  colocated  part  of  the 
plant  and  G^  to  represent  the  noncolocated  part  of  the  plant  will  allow  an 
intuitive  feel  for  the  results.  The  assumption 

||G]  +0^1  <  1 1  F,  1 1  and  ||G2  +  D2||  <  ||F2|| 


can  be  interpreted,  for  example,  as  the  uncertainty  model  for  the  colocated 
and  noncolocated  parts  of  the  system.  Now,  in  order  for  the  Conicity  Theorem 
to  hold,  the  second  condition 


F  ( I-HD) 


-1 


<  1 


% 


must  be  satisfied.  Substituting  the  variables  corresponding  to  Case  1 
yields: 


■1 


LH2J 


<  1 


[pi  !  F2] 


I  + 


C°1  i  °2] 


V 

r1 

lA'J 

1 

<  1 


(F1  H1  +  F 2  H2)(I  +  Dl  Hi  +  d2  H2) 


-1 


<  1 


F1  H1  +  F2  H2 


(I  +  D]  H]  +  D2  H2) 


-1 


Substituting  the  variables  corresponding  to  Case  2  yields: 

?:  ! 1  F  ["i  i  ¥  [-oil  !  1  [H,  i  h2]  I  <  1 


f’-  [Hi  i  V  j 1  *  -o’-  Ch,  j  h2]  [ _1  <i 


F1  H1  F1  H2 
F2  H1  F2  H2 


I  +  H]  0] 

h2 


H1  °2 
i  +  h2  d2 


If  one  now  chooses  an  norm  in  these  expressions,  Parseval's  Theorem 
may  be  used  to  evaluate  them.  This  results  in 

Case  1  : 

(Fi  H1  +  F2  H2}*  (F  H1  +  F2  H2}  j  <  \nin  J(I  +  H1  D1  +- H2  °2}*  (l  +  Hi  Di  +  H2  V  l 


Case  2: 


F1H1iF1H2]  r  F1  H1  '  F1  H2  "I  f  irI+Hl  D1  '  H1  D2l*rI+Hl  “l  '  H1  V 


where  Xmax  and  Xmin  denotes  take  the  maximum  and  minimum  eigenvalues 
of  {•},  respectively. 

Note  that  it  has  been  assumed  that  G] ,  G2,  F] ,  Fz,  D] ,  D2,  H]  and  H2 

have  a  Laplace  or  Z-Transform  representation,  with  the  variable  S  or 
ST 

Z=e  as  the  argument.  The  conditions  are  to  be  evaluated  at  every  fre¬ 
quency  S,  where  S=ju>  and  «  ranges  from  0  to  »  for  continuous  systems  and 
from  0  to  tt/T  for  sample  data  systems. 

Note  also  that  these  conditions  degenerate  to  the  positivity  conditions 
in  ACOSS-8  if  the  system  is  square  and  tiie  conditions  on  page  8-12  are  imposed. 
The  latter  are  not  possible  to  enforce,  however,  in  the  nonsquare  case. 


8.4  Squaring  Down  Procedure 

As  explained  before,  the  design  technique  adopted  by  TRW  combines 
positivity  theory  with  the  multivariable  frequency  domain  approach.  Both 
of  these  techniques  were  originally  defined  for  square  systems,  i.e.,  systems 
for  which  the  number  of  actuators  equals  the  number  of  sensors.  Since  in 
many  large  space  structures  this  is  not  the  case,  the  system  should  be 
"squared  down"  before  the  usual  design  procedure  is  followed.  Thus,  the 
design  procedure  for  rectangular  systems  calls  for  the  following  two  stage 
procedure. 

(i)  Calculation  of  the  Square  Down  Filter.  The  square  down  filter  K(Z) 
is  chosen  to  improve  on  the  performance  of  the  system.  This  is  done  by 
designing  K(Z)  to  discriminate  against  modes  which  have  limited  influence 

on  the  performance  of  the  system,  e.g.,  do  not  affect  the  motion  of  line  of 
sight  (LOS),  in  favor  of  modes  which  affect  the  performance  the  most.  In 
addition,  K(Z)  should  be  constructed  to  enhance  the  stability  and  robustness 
characteristics  of  the  system  which  are  both  related  to  its  positivity 
index.  This  can  be  done  by  designing  K(Z)  in  a  way  which  makes  the  re¬ 
sultant  square  system  as  close  as  possible  to  a  positive  real  system.  As 
usual  in  design  problems,  a  suitable  compromise  must  be  made  between  the 
conflicting  requirements  of  stability  and  performance. 

(ii)  Design  of  Controller.  After  the  system  is  squared  down,  a 
compensator  G c(Z)  is  designed  using  the  positivity  design  approach.  Since 
the  square  down  filter  K(Z)  is  designed  to  enhance  performance,  the  important 
modes  will  be  prominently  represented  in  the  feedback  signal.  As  a  result., 
the  damping  of  these  modes  can  be  increased  using  lower  controller  gains. 

It  should  be  noted  that  dividing  the  control  system  into  two  separate 
parts,  controller  +  square  down  filter,  is  artificial  and  was  introduced 
here  only  for  the  explanation  of  the  proposed  design  procedure.  The  complete 
control  law  which  is  finally  implemented  is  given  by  H(Z)  =  Gc(Z)  K(Z) 
where  H(Z)  is  a  rectangular  digital  filter,  as  shown  in  Figure  8-10. 
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Figure  8-10.  Topology  of  a  "Squared  Down"  System 


It  should  be  also  noted  that  the  squaring  down  procedure  can  also  be 
used  to  augment  the  representation  of  the  important  modes  in  the  feedback 
signal  of  a  square  noncolocated  system.  This  can  be  done  prior  to  designing 
the  controller,  thus  facilitating  the  task  of  damping  augmentation. 

8.4.1  Positivity  Properties  of  Large  Space  Structures 

Before  explaining  how  colocation  of  actuators  and  sensors  makes  LSS 
behave  as  a  positive  operator,  various  dynamic  representations  of  vibrational 
systems  will  be  introduced.  These  dynamic  descriptions  will  be  later  used 
in  defining  various  characteristics  of  the  system  at  hand. 


In  the  following,  three  different  descriptions  of  a  large  space  structure 
will  be  used  alternately. 


(i)  Modal  Representation.  Motion  of  a  linear  vibratory  system  is  given 


by: 


Mq  +  Kq  =  u 


(8.3) 


where  M  is  a  positive  definite  nxn  mass  matrix,  K  is  a  non-negative, 
symmetric  n  xn  stiffness  matrix  and  u  is  a  time  dependent  forcing  function. 


For  a  complicated  structure,  the  matrices  M  and  K  are  defined  using  the 
finite  element  approximation  method.  Modal  analysis  is  then  performed  to 
compute  a  transformation  matrix  $  which  relates  the  physical  coordinates  q 
to  the  modal  coordinates  n — >q  =  <j>n.  The  equations  of  motion  can  be  then 


written  as: 


n  +  fin  =  <p  u 


(8.4) 


where  n  =  diag  (ui^ . 

Let  the  actuators  and  sensors  locations  be:  a, _ a  .  and  s. ...s 

1  m  1  r 

respecti vely.  Then  after  including  modal  damping  and  assuming  rate 
measurements,  the  equations  of  motion  take  the  form: 


n.j  +  2?aii  ni  +  to.*  ni  =  Uj  +  . 

yk  *  +1  (si>  "l  +  (sk>  % 


+  <t> .  (a  )  u 
M  v  nr  m 


(8.5) 

(8.6) 


i  =  1. . .  .n 
k  =  1. . . .  r 

where  ?  is  the  damping  coefficient  and  y  i.s  the  observation  vector, 
(ii)  State  Space  Representation 


Defining  x  4  (n^, 
written  as: 


.nn»  np)  equations  (8.5)  and  (8.6)  can  be 


x  =  Ax  +  Bm 


(8.7) 


I 


C=  0  C.  0  c0..o  c 
12  n 


r  =  1  j- 

i 


♦i  <sr> 


( i  i i )  Transfer  Matrix  Representation 

By  definition  the  transfer  function  of  a  system  defined  by 
equations  (3.7)  and  (8.8)  is  given  by: 


G  (s)  =  C  (sI-A)-1  B 


Defining 


B*=  B1  , 


c*=  [cr  .  .  .cn] 


(8.9) 


equation  (8.9)  can  be  written  as: 

Gp(s)  =  C*  D(.s)  B*,  D(s 


G  (s)  =  C*  D(s)  B*,  D(s)  =  diag  (  - 5 - 

s^+Z^.s+o)/ 


(8.10) 


As  explained  in  Section  8.1.2,  the  positivity  of  a  square  system 
can  be  tested  according  to  the  definition  by  evaluating  the  matrix 
G(jto)  +  G*( joj)  for  continuous  systems  [or  G(e^“T)  +  G*(eJ'“T)  for 
discrete  systems]  for  0  <  to  <  »  ,  and  evaluating  the  positivity  index. 

For  multi-input  multi-output  (MTMO)  systems  the  frequency  domain  approaci 
requires,  however,  extensive  computing  effort  to  calculate  the  eigen¬ 
values  of  the  above  matrices  for  many  frequency  values.  To  overcome  this 
difficulty  R.D.O  Anderson  developed  a  relatively  simple  positivity  test 


in  the  time  domain.  The  theorem  on  which  the  test  is  based  is  a  general 
ization  of  the  Kalman  -  Yakubovich  lemma  for  single-input  single  output 
(SISO)  systems  and  can  be  stated  as  follows: 

Theorem:  Let  G(s)  be  a  square  matrix  of  rational  transfer  functions 
such  that  G(»)  is  finite  and  G  has  poles  which  lie  in  Re  [s']  <  0  or  are 
simple  on  Re  [s]  =  0.  Let  {A,  B,  C,  G(»)}  be  a  minimal  realization  of 
G.  Then  G(s)  is  positive  real  if  and  only  if  there  exists  a  symmetric 
positive  definite  matrix  P  and  matrices  WQ  and  L  such  that 

(1)  PA  +  ATP  =  -LLT 


(ii)  WoTWo  =  G(°°)  +  G^"(°°) 


(iii ) 


PB  +  LW_ 


Further,  G(s)  is  strictly  positive  real  if  LL  is  positive  definite. 

In  the  case  of  a  large  space  structures  G(®)=0,  equation  (ii) 
reduces  to  WQ=0  which  results  in  equation  (iii): 


CT  =  PB 


(8.11) 


Positivity  of  a  system  with  colocated  actuators  and  sensors  can  be 
established  in  the  frequency  domain  by  observing  that  for  a  colocated 
actuators  and  rate  sensors  one  has: 


T  * 
C*  =  B 


Hence  from  equation  (8.10) 


Gp(jco)  +  Gp*( jw )  =  B*T  E(«)B* 


(8.12) 


(8.13) 


where 


E(o»)  =  diag  [Re( — ~ ^ - 7)] 

L  -W  +2?  Juju  .+U)  •  J 


E (a) )  is  positive  definite  for  all  w  and  thus  Gp(joj)  +  Gp*  (ju)  is 
also  positive  definite. 


In  the  time  domain  it  can  be  shown  that: 


where 


2 

,Jji  pi 

0 

’  L.  = 

0 

0 

0 

pi 

i 

0 

^aP  j 

_ 

— 

l 

- 

J 

(8.14) 


(8.15) 


is  a  solution  to  the  equations  of  Anderson's  theorem.  This  is  provided 
the  actuators  and  rate  sensors  are  colocated  and  have  infinite  bandwiths 

The  above  solution  is  valid  for  any  number  of  modes.  Hence,  if  the 
dynamics  of  the  actuators  and  rate  sensors  are  not  included,  the  spill¬ 
over  modes  cannot  destabilize  the  system.  This  is  because  any  number 

of  spillover  modes  form  a  positive  system,  i.e.,  a  system  which  is 
energy  dissipative  and  therefore  stable. 

The  situation  for  noncolocated  systems  is  not  that  fortunate.  For 
SISO  noncolocated  systems  if  sign  (a^  )j  =  sign  [^.(s^J  (8.16) 

then  the  system  is  positive.  This  can  be  shown  by  taking  _  Msl^ 

for  i  =  1, .  This  gives  which  is  positive  1  ‘h ^al^ 

definite  and  a  /pT  which  is  real.  Obviously,  there  are  some  high 
frequency  modes  for  which  this  condition  cannot  be  satisfied.  However, 
with  a  judicious  placement  of  actuators  and  sensors  these  modes  can 
be  placed  outside  the  control  bandwidth  and  therefore  will  not  affect 
system  stability. 

Intuitively,  one  would  like  to  extend  the  above  result  to  MIMO 
noncolocated  systems.  Unfortunately  it  can  be  shown  that  satisfying 
the  equation 

sign  [>.(aj)]  =  sign  (sj)]  i  =  l.tJ 


(8.17) 


for  j=l,...m  does  not  assure  a  positive  system. 

MIMO  structures  which  can  be  shown  to  be  positive  are  classes  of 
structures  with  colocated  actuators  and  sensors.  Although  non- 
colocated  MIMO  structures  are  not  generally  positive,  it  can  be  shown 

that  when  a  loop  is  closed  on  a  noncolocated  system,  the  system  is  stable 

IT  - 

provided  (Bo  +  C  B)  is  a  semi  positive  definite  matrix. 

8.4.2  Squaring  Down:  Time  Domain  Approach 

The  time  domain  square  down  procedure  makes  use  of  Anderson's 
theorem  for  establishing  positivity.  If  the  number  of  controlled  modes 
nc  satisfies  nc=max(r,  m^a  constant  square  down  matrix  K  which  makes 
the  system  positive  can  be  obtained.  If  nc  =  r,  then  using  Anderson's 
lemma  one  can  write: 

P  B  =  CT  K  T  (8.18) 

P  is  chosen  as  in  equations  (8.14)  and  (8.15)  and  because  of  the 
peculiar  form  of  B,  equation  (8.18)  reduces  to: 

KT  =  (C*1)'1  B*  (8.19) 

Similarly  if  nc  =  m,  it  can  be  shown  that  the  a  square  down  matrix  which 
makes  the  system  positive  is  qiven  by: 

K  =  B*"1  CT  (8.20) 

The  requirement  nc  =  max  (r,m)  leads  to  a  high  number  of  sensors 
(or  actuators).  If  one  is  forced  to  limit  the  number  of  sensors  and 
actuators,  the  time  domain  approach  can  still  be  pursued  by  creating  a 
square  system  which  is  "close"  to  a  square  positive  system.  Let 
nc  =  r  then  using  the  Moore  -  Penrose  pseudo  inverse,  the  square  down 
constant  matrix  is  given  by 

K  =  (CV  B*  (8.21) 

where  +  denotes  a  pseudo-inverse.  Using  K  as  given  by  equation 
(8.21)  generates  a  square  system  that  is  closest  to  a  positive  square 
system  in  a  Euclidean  norm  sense.  In  the  program  which  implemented  this 
approach,  the  pseudo  inverse  was  calculated  using  a  singular  value  sub¬ 
routine  taken  from  the  UNPACK  library. 


In  the  time  domain  approach  the  square  down  criterion  is  to  make  the 
system  as  positive  as  possible  without  paying  attention  to  the  perfor¬ 
mance  requirements.  The  performance  characteristics  of  the  system  are 
designed  for  in  the  second  stage  of  the  design  by  calculating  an  app¬ 
roximate  compensation  filter  G  (z). 

8.4.3  Squaring  Down:  Frequency  Domain  Approach 

In  this  approach  the  square  down  matrix  is  allowed  to  depend  on 
frequency.  Although  it  is  possihle  to  calculate  K  for  any  given 
frequency,  the  calculation  will  be  carried  out  only  at  the  natural 
frequencies  because  at  these  frequencies  a  sudden  jump  in  system's  gains 
occurs.  In  the  approach  taken  here,  a  DC  value  of  K,  KQ  =  K(o)  is 
first  calculated  and  then  a  digital  filter  is  added  to  adjust  the 
compensator's  gain  at  modal  frequencies. 

The  DC  value  of  the  squaring  matrix  K  is  calculated  from  the 
equation: 

Kc  Gp{0)  =  Ro  (8.22) 

The  matrix  Gp(o)  is  the  value  of  the  transfer  matrix  at  oi  =  0  and  is 
given  by: 

n 

Gpij^0)  =  X  ^k  (sJMaJ  (8.23) 

k=l  J 

Another  interpretation  of  each  element  of  Gp(o),  G  (o),  is  that  it 
expresses  the  motion  at  location  s.  due  to  a  unit  force  at  location  aj. 

The  purpose  of  this  square  down  scheme  is  to  scale  the  feedback  signals 
according  to  the  motion  of  the  LOS.  This  is  done  by  calculating  a 
diagonal  matrix  Rq  according  to  the  equation: 

Jj  ^  /  2  p  ■  2*™ 

rjj  =  X  ^k(aj}  ^klx  +  *kly'  +  ^klz“  ,  j=l...m  (8.24) 

k=l 


r . . 


.1 .1 


-  is  the  motion  exerted  by  a  unit  force  at  actuator  j  on  the  LOS. 


Since  Kq  is  a  rectangular  matrix  the  system  of  equations  defined  by 
(8.22)  is  underdetermined  and  many  solutions  are  possible.  In  the 
computer  program  based  on  this  approach,  the  matrix  Kq  is  calculated 
using  a  singular  value  decomposition  which  results  in  a  minimum  norm 
sol ution. 

In  the  second  part  of  the  square  down  procedure  the  frequency 
dependence  of  the  matrix  K  is  computed.  This  is  done  by  calculating  the 
norm  of  the  transfer  matrix  at  the  natural  frequencies,  a.  and  equating 
it  to  the  modal  participation  factor  -g^  at  the  LOS  evaluated  at  this 
frequency.  The  norm  of  the  transfer  matrix  ||K  G(oi^)||  at  a  given 
frequency  is  taken  as  the  maximum  singular  value.  Since  the  underlying 
assumption  is  that  the  modal  damping  is  very  low,  there  is  only  one 
singular  value  of  Kq  G(oj^)  which  is  different  from  zero. 

The  modal  participation  factor  is  given  by: 


.  gi 

The  ratio  —  gives  the  gain  required  at  each  frequency.  This  gain  should 

be  looked  upon  as  a  measure  of  the  relative  importance  of  various  modes 
as  reflected  in  the  LOS  motion.  After  this  calculation  is  completed  for 
U  =  to.,  i  =  l...n  a  sequence  of  gains  is  obtained.  The  next  step  is  to 
approximate  this  series  of  gains  a(aj-),  i  =  l...n  by  a  realizable  digital 
filter  a ( z ) .  The  resulting  square  down  matrix  is  given  by  a(z)-KQ.  Since 
frequency  domain  design  approach  results  in  a  diagonal  controller 

Gc(z)  *  8(z)  •  I  (24) 

The  resulting  complete  controller  H(z)  is  given  by 

H(z)  =  a(z)  3(z)  Kq  (25) 

The  main  criterion  of  the  square  down  procedure  based  in  the  fre¬ 
quency  domain  is  to  enhance  system  performance.  The  relative  stability 
characteristics  of  the  closed  loop  system  still  has  to  be  established 
while  designing  the  compensator  G  (z).  Therefore,  this  procedure  can  be 
employed  to  increase  the  representation  of  the  important  modes  for 


square  noncolocated  systems  prior  to  the  usual  compensator  design  pro¬ 
cedure. 

An  important  special  case  of  non- square  systems  is  the  case  of  a 
system  which  possesses  a  colocated  subsystem  of  actuators  and  sensors, 

This  is  a  very  realistic  case  since  the  main  difficulty  in  creating  a 
colocated  control  system  is  in  placing  the  actuators  due  to  size  and 
power  limitations.  For  this  class  of  structures  a  preferred  design  approach 
is  first  to  create  a  colocated  system  by  placing  sensors  at  the  actuators 
locations  and  then  place  additional  sensors  at  critical  locations  (e.g.,  at 
LOS)  to  enhance  system's  closed  loop  performance. 

For  the  class  of  systems  with  partly  colocated  actuators  and  sensors, 
the  structure  of  the  square  down  filter  is  set  to  preserve  the  favorable 
stability  character!- sties  of  the  colocated  subsystem. 

Assume  that  the  complete  system  possesses  r  sensors  of  which  m  are 
colocated  with  the  m  actuators.  Calculating  the  square  down  matrix  as 
in  Section  8.4.3  would  destroy  the  positivity  of  the  colocated  sub¬ 
system.  Instead  let  K  have  the  following  structure 

K(z)  =  [i  ;  L  (zj]  (8.28) 

The  above  structure  will  ensure  that  the  colocated  subsystem  remains 
positive  real.  This  form  of  K  will  be  referred  to  as  "partial  square 
down"  matrix.  Figure  8-11  shows  the  functional  diagrams  of  both  the 
full  and  partially  squared  down  systems. 

If  r<  2m  the  system  of  equations  to  determine  the  DC  value  of  L, 

L  ,  is  underdetermined.  In  this  case  a  least  square  procedure  is  em¬ 
ployed. 

A  digital  program  based  on  the  frequency  domain  approach  square 
down  technique  was  used  on  the  flat  plate  and  tetrahedron  examples.  The 
flowchart  of  this  program  is  shown  in  Figure  8-12.  The  program  is 
interactive  and  allows  for  both  the  general  and  partial  squaring  down 
matrix  calculation.  The  program  reads  the  system  state  space  description 
from  a  preassigned  tape  and  the  user  supplies  the  format  of  the  desired 
square  down  matrix.  The  output  of  the  program  consists  of  a  DC  matrix 
K  (or  L  )  and  the  gains  of  the  modal  frequencies.  In  addition,  the  state 


9.0  DESIGN  EXAMPLES  AND  EVALUATION 


Two  design  examples  were  selected  to  demonstrate  the  approach 
discussed  in  Chapter  8.  The  two  examples  chosen  are  the  TRW  flat  plate 
and  the  Draper  Tetrahedron.  The  design  in  both  examples  use  noncolocated 
control  systems  and  thus  are  suitable  for  demonstrating  the  proposed 
design  procedure.  In  the  following,  the  structures  will  be  described, 
then  the  design  steps  will  be  explained;  robustness  properties  of  the 
designs  will  be  evaluated,  and  finally  the  performance  of  the  proposed 
control  systems  for  each  structure  will  be  evaluated  using  digital 
simulation . 

9.1  The  TRW  Flat  Plate  Example 

The  flat  plate  design  example  is  an  existing  experimental  set  up 
shown  in  Figure  9-1.  The  plate  is  1.73  x  1.22  meters  and  is  made  of 
1.6  mm  aluminum.  The  top  edge  of  the  plate  is  clamped  to  a  supporting 
structure  which  is  rigidly  attached  to  the  wall.  Control  is  applied 
to  the  plate  by  using  the  bending  moment  actuators  shown  in  the  figure. 
Each  actuator  has  a  built  in  magnetic  sensor  which  measures  relative 
deflections  and  rates. 

The  model  used  to  design  the  controller  was  generated  using  TRW's 
finite  element  program  TRWSAP.  The  modal  frequencies  of  this  "nominal 
model"  are  shown  in  Table  9-1.  Next  to  these  theoretical  frequencies 
are  the  true  frequencies  of  the  plate  (our  "perturbed  model").  The 
assumed  and  actual  damping  are  also  shown. 


Table  9-1.  Theoretical  and  Experimental  Model 

Frequencies  and  Damping  Coefficients 


Theoretical 
Freq  (Hz) 

Experimental 
Freq  (Hz) 

.413 

.03 

.67 

1  .27 

.03 

1 .67 

2.83 

.03 

2.75 

4.96 

.03 

5.14 

5.87 

.03 

5.43 

8.08 

.03 

8.0 

10.08 

.03 

9.7 

10.5 

.03 

11.0 

12.3 

.03 

12.25 

.06 

.009 

.003 

.032 

.042 

.023 

.062 

.013 

unknown 


9.1.1  Flat  Plate:  Control  System  Design 

In  this  section  the  designs  using  colocated  and  noncolocated  control 
systems  are  described  and  evaluated.  The  design  is  done  using  the  multi- 
variable  frequency  domain  technique  combined  with  the  positivity  approach. 
The  square  down  procedure  described  in  Section  8  is  also  used  for  the 
noncolocated  design.  The  design  evaluation  criteria  used  are  the 
singular  value  decomposition  (SVD)of  the  return  difference  and  the 
positivity  under  the  system.  Both  measure  the  system's  robustness. 

The  assumption  used  in  the  colocated  design  was  that  only  the  two 
colocated  actuators  and  sensors  shown  in  Figure  9-1  are  to  be  used  for 
controlling  the  X  and  Y  motions  of  the  plate.  In  the  noncolocated  design, 
an  additional  sensor  and  an  inertial  rate  sensor  were  placed  at  the 
bottom  right  corner  of  the  plate  (this  sensor  has  not  been  implemented 
yet  on  the  experimental  set  up).  This  location  was  chosen  for  the 
additional  sensor  because  there,  the  low  frequency  modes  are  more 
observable. 

The  goal  of  both  the  colocated  and  noncolocated  designs  are  to 
augment  the  modal  damping  so  that  the  plate  motion  settles  quickly, 
in  response  to  initial  conditions. 

9.1.2  The  Colocated  Flat  Plate  Controller 

The  colocated  controller  developed  in  TRW  for  the  experimental 
set  up  shown  in  Figure  9-1  is  given  by: 

60.94  +  J2-4Z  -  12.1 - 

r  -  1 .77Z  +  0.884 

The  controller  operates  in  each  of  the  two  system's  channel  with  sampling 
time  Ts  =  0.01  seconds.  The  design  was  performed  using  the  positivity 
and  multivariable  frequency  domain  approach.  The  open  loop  character¬ 
istic  gains  for  the  compensated  and  uncompensated  colocated  systems 
in  the  w-plane  are  shown  in  Figure  9-2.  The  minimum  and  maximum  singular 
values  of  the  return  difference  and  a  plot  of  the  positivity  index  are 
shown  in  Figure  9-3. 
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9.1.3  Noncolocated  Flat  Plate  Desi 


As  explained  in  Section  8.4,  the  first  step  of  the  noncolocated 
design  is  to  square  down  the  system.  This  was  done  by  applying  the 
procedure  described  in  Section  8.4.3  and  using  the  program  based  on  it. 
The  square  down  of  the  sytem  was  performed  us i the  corner  sensor 
response  as  the  performance  criterion.  This  criterion  emphasizes  the 
first  two  vibrational  modes.  The  resultant  square  down  matrix  is  given 
by: 

1.4  0  -30.834 

K  = 

0  10.  -16.977 

The  second  step  of  the  design  is  to  compute  the  characteristic  gains 
X(v)  of  the  resultant  squared  down  open  loop  system  in  the  W-plane. 

They  are  shown  in  Figure  9-4 (a ) .  Note  that  the  characteristic  gains  are 
computed  after  the  squaring  down  procedure  of  step  1  since  they  are 
defined  only  for  square  systems. 

The  third  step  of  the  design  is  to  shape  the  characteristic  gains  and 
the  plant's  positivity  under  to  augment  modal  damping  and  the  robustness 
properties  of  the  system.  The  resultant  compensator  is: 

Rm  -  0-0193Z2  +  0. 032437Z  +  0.016218  T  (9.3) 

OK*-]  O  • 

T  -  1 .64706Z  +  0.716807  ^ 


Hence  the  complete  controller  is  given  by: 

0.0193Z2  +  0.032437Z  +  0.016218  I  1  .4 


-30.834 


1 .64706Z  +  0.716807 


0  10.  -16.977 


(9.4) 


The  compensated  systems  characteristic  gains  are  shown  in  Figure  ’-4(b) 
The  maximum  and  minimum  singular  values  of  the  return  difference  and 
the  positivity  index  plot  of  the  noncol located  system  are  shown  in 
Figure  9-5. 

9.1  .4  Noncolocated  Flat  Plate  Design  Evaluation:  Robustness 
Properties  and  Closed  Loop  Poles 

The  robustness  properties  of  the  closed  loop  system  were  assessed 
by  two  methods:  the  positivity  index  5 (co )  and  by  computing  the  minimum 
and  maximum  singular  values  of  the  return  difference  matrix.  The 
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(b)  X(v)  OF  OPEN-LOOP,  COMPENSATED,  NONCOLOCATED  SYSTEM 


Figure  9-4.  Characteri Stic  Gains  of  the  Uncompensated 
and  Compensated  Systems 
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(a)  POSITIVITY  INDEX 


(b)  MINIMUM  AND  MAXIMUM  SINGULAR  VALUES  OF  THE 
RETURN  DIFFERENCE 


Figure  9-5.  Positivity  Index  and  SVD  of  the  Return  Difference 
for  the  Squared  Noncol ocated  Design 


singular  value  analysis  was  performed  on  the  return  difference  matrix 
since  an  additive  perturbation  was  assumed  (see  [22],  [23]),  i.e.  the 
perturbed  system's  transfer  matrix  is  assumed  to  be  in  the  form: 

G  (s)  =  G  (s)  +  E  (s)  (9.5) 

Where  E(s)  is  the  difference  between  the  nominal  and  perturbed  system 
It  can  be  shown  (see  [22],  [23],)  that  a  sufficient  condition  for 
stability  is  given  by: 

£  [I  +  G(jaj)]  >  a  [E(jio)]  Vu)  >  0  (9.6) 

where  a  and  a  denote  the  minimum  and  maximum  singular  values, 
respectively. 

It  should  be  noted  that  the  singular  value  analysis  is  a  very 
conservative  test  of  robustness  for  oscillatory  systems.  This  is  due 
to  the  existence  of  many  lightly  damped  vibrational  modes  that  result 
in  sharp  peaks  in  the  frequency  plot  of  a  [E(ju)].  Since  c  .  [I  +  G(jw)] 
has  to  be  larger  than  amax  [E(jm)]  for  all  frequencies,  the  resonance 
peaks  impose  very  conservative  bounds  on  robustness. 

As  shown  in  the  plots  of  the  positivity  index  in  Figures  9-3(a) 
and  9-5(a) ,  the  colocated  and  noncolocated  designs  will  operate  at 
approximately  the  same  DC  gain  of  1.  These  figures  show  also  that 
both  systems  tolerate  moderate,  low  frequency  modal  uncertainties. 
Moreover,  both  designs  can  withstand  large,  high  frequency,  modal 
uncertainties.  It  is  also  observed  that  the  low  frequency  robustness 
of  the  noncolocated  design  is  poorer  than  the  colocated  design,  but, 
this  situation  is  reversed  at  mid  frequencies.  Hence  the  positivity 
index  analysis  conclusion  is  that  the  noncolocated  design  has  slightly 
better  robustness  properties  as  it  is  more  tolerant  of  uncertainty  at 
the  mid  and  high  frequencies  where  the  model  i_s  less  certain. 

The  singular  value  plots  of  the  colocated  and  noncolocated  designs 
shown  in  Figures  9- 3 ( b )  and  9- 5 ( b ) ,  reinforces  the  conclusion  of  the 
positivity  index  analysis.  The  plots  show  that  the  colocated  design 
will  tolerate  moderate,  additive,  but  otherwise  unstructured  uncertainties 
over  the  band  v  =  0.  to  0.1.  Whereas  the  noncolocated  system  will 
withstand  moderate  perturbations  over  the  band  v-0.  to  0.4.  For  about 


a  decade  later,  both  control  systems  can  tolerate  small  unstructured 
uncertainties.  Finally,  at  high  frequencies,  the  location  of  the  negative 
peaks  in  the  plots  leads  to  the  same  conclusion  as  the  positivity  under 
analysis,  namely,  that  overall  the  noncolocated  design  is  slightly  better. 

To  complete  the  evaluation  of  the  proposed  control  designs,  the  poles 
of  the  closed  loop  system  were  computed.  The  list  of  the  closed  loop  poles 
for  both  the  colocated  and  noncolocated  designs,  using  the  theoretical 
model,  is  given  in  Table  9-2.  It  can  be  observed  that  the  noncolocated 
design  has  better  low  frequency  damping  characteristics  than  the  colocated 
system  but  the  situation  reverses  at  mid  frequencies.  At  high  frequencies, 
the  nearly  positive  character  of  the  colocated  system  becomes  evident, 
having  the  high  frequency  modes  much  better  damped. 

9.1.5  Simulation  Study  Results  for  the  Flat  Plate  Example 

A  digital  simulation  study  was  conducted  to  evaluate  the  performance 
of  the  colocated  and  noncolocated  control  systems  using  both  the 
theoretical  and  perturbed  {"real")  models.  Figure  9-6  compares  the 
uncompensated  response  of  the  flat  plate  with  those  of  the  colocated 
and  noncolocated  designs  for  the  theoretical  model  of  the  plate.  The 
plots  shown  are  the  modal  position  response  of  the  plate  when  the  initial 
condition  of  the  first  ten  modes  is  set  to  1.  They  verify  the  analysis 
and  clearly  demonstrate  the  improved  damping  character  of  the  noncolocated 
design:  the  first  two  modes  are  decaying  much  faster  than  in  the  colocated 
case.  Figure  9-7  shows  the  corresponding  responses  of  modes  6  to  1 5  for 
the  noncolocated  design  with  the  initial  condition  as  above.  It  is 
clear  that  all  the  high  frequency  modes  are  stable.  The  conclusion 
is,  therefore,  that  the  simulation  study  verifies  the  analysis  results 
for  the  theoretical  model  as  far  as  stability,  and  modal  damping  of  the 
modes  are  concerned. 

A  similar  simulation  study  was  performed  on  the  perturbed  flat  plate 
model.  Figure  9-8  shows  the  comparison  between  the  open  loop  performance 
and  noncolocated  designs  as  applied  to  the  "real"  plate.  It  is  evident 
that  the  qualitative  performance  of  noncolocated  designs  remains  the  same: 
it  has  good  low  frequency  damping.  Also,  it  is  clear  that  all  the  modes 
remain  stable  despite  the  mismodel 1 ing,  as  shown  in  Figure  9-9. 
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MODAL  POSITION  (11-15)  MODAL  POSITION  (6-10) 


Figure  9-9.  Modal  Response  of  the  Perturbed  Model  with 
Noncol ocated  Control  System:  Modes  6-15 
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9.2  The  Draper  Tetrahedral  Truss  Example 

The  Tetrahedral  Truss,  shown  in  Figure  9-10,  was  selected  as  the 
second  test  example  the  noncolocated  design  approach.  This  structure 
was  devised  by  Draper  Labs  [24]  as  one  of  the  simplest  non-planar 
geometries  capable  of  representing  a  large  space  structure.  Despite  its 
apparent  simplicity,  it  models  the  feed-tower  in  a  generic  class  of 
large  antenna  application.  The  tetrahedral  apex  represents  the  antenna 
feed,  members  1-6  are  a  part  of  the  feed  support  structure,  and  bi-peds 
7-8,  9-10,  and  11-12  are  supports/controls  attached  to  an  inertially 
stabilized  (assumed)  antenna  dish. 

The  finite  element  method  was  used  to  model  the  tetrahedron. 

Referring  to  Figure  9-10,  the  model  contains  ten  nodes,  each  with  three 
degrees  of  freedom,  and  twelve  truss  members.  All  truss  members  are 
pin  connected  and  are  flexible  in  the  axial  direction  only.  Nodes  5 
through  10  are  pinned  to  a  rigid  support.  Masses  are  lumped  at  nodes 
1  through  4,  each  with  3  degrees  of  freedom.  Tables  9-3  and  9-4  show 
the  elements  areas  and  the  lumped  mass  for  both  the  nominal  and  perturbed 
models  and  Table  9-5  shows  the  modal  frequencies  for  both  these  models. 
Modal  damping  of  ^0.5%  was  assumed  for  all  modes. 

9.2.1  Control  System  Design  for  the  Tetrahedral  Truss 

The  goal  of  the  control  system  design  was  to  come  up  with  a  controller 
designed  according  to  the  first  8  modes  of  the  nominal  structure.  The 
design  was  evaluated  by  robustness  measures  such  as  positivity  index 
and  SVD  analysis  and  by  the  closed  loop  eigenvalues.  In  addition,  the 
time  history  response  to  the  initial  conditions  given  in  Table  9-6 
was  simulated.  The  objective  is  that  the  X-L0S  error  and  Y-L0S  error 
are  less  than  4.  x  10  ^  and  2.5  x  10  ^  respectively  in  20  seconds.  The 
design  will  be  evaluated  with  respect  to  the  12  mode  full  models  of  the 
nominal  and  perturbed  structures. 

Two  control  system  designs  are  compared: 

(i)  A  colocated  system  design  -  Six  colocated  sensors  and  actuators 
are  assumed  to  act  as  member  dampers  in  parallel  with  members  7  through  12 
The  member  dampers  were  assumed  capable  of  sensing  relative  position  and 
velocity  and  exerting  force  in  the  axial  direction  only. 


LEG  1  STATION  1 


Figure  9-10.  Draper  Tetrahedral  Truss 


Table  9-3.  Tetrahedron  Model's  Element  Areas 


El  ement 

Nominal  structure 
cross-section 

Perturbed  structure 
cross-section 

1 

1000 

1200 

2 

100 

150 

3 

100 

150 

4 

1000 

5 

1000 

6 

1000 

7 

100 

150 

8 

100 

150 

9 

100 

150 

10 

100 

150 

11 

150 

12 

100 

150 

Table  9-4.  Tetrahedron  Model's  Lumped  Masses 


Node 

Nominal  structure 
mass 

Perturbed  structure 
mass 

1 

2 

4 

2 

2 

2 

3 

2 

2 

4 

2 

2 

(ii )  Noncol ocated  system  design  -  In  addition  to  the  colocated 
system,  two  sensors  which  measure  X-LOS  and  Y-LOS  motions  were  assumed. 

This  created  a  noncolocated  control  system  with  6  actuators  and  8  sensors. 

9.2.2  The  Colocated  Controller  for  the  Tetrahedron 

The  colocated  controller  for  the  tetrahedron  was  developed  using 
the  design  methodology  developed  in  Section  8  and  is  given  in  Appendix  D. 
The  sampling  rate  was  chosen  to  be  0.2  seconds.  This  rate  is  approximately 
three  times  faster  than  the  highest  controlled  mode.  In  order  to  achieve 
the  desired  rate  of  decay,  the  colocated  controller  was  designed  using 
an  a-shifted  plane  as  explained  in  Appendix  D.  The  digital  filter  was 
designed  using  the  generalized  Nyquist  Theorem  for  multivariable  systems. 
The  decentralized  controller  resulting  from  the  design  in  the  w-plane  is: 

H(Z)  -  I  11  .492  +  3..-.921162Z  -  3,9863812 
X  Z 2  -  1. 70921 Z  +  .846578 

The  open  loop  characteristic  gains  for  the  compensated,  colocated  system 
are  given  in  Figure  9-11.  The  minimum  and  maximum  singular  values  of 
the  return  difference  matrix  and  a  plot  of  the  positivity  under  are 
given  in  Figure  9-12.  The  closed  loop  eigenvalues  are  shown  in  Table  9-7. 

9.2.3  Design  of  the  Noncolocated  Control  System  for  the  Tetrahedral 
Truss 

As  in  the  flat  plate  example, the  first  step  in  the  design  of  a 
noncolocated  nonsquare  control  system  is  to  square  down  the  plant  by 
calculating  on  appropriate  matrix  KQ.  Both  a  full  order  and  partial 
squaring  down  matrices  were  generated  and  are  given  in  Table  9-8  (i) 
and  (ii).  Appropriate  filters  were  also  generated  to  improve  the  system's 
performance.  As  it  will  be  shown  later,  these  filters  are  not  needed 
in  the  design  and  constant  matrices  yield  good  closed-loop  performance 
which  meets  requirements.  Figure  9-13  compares  the  frequency  responses 
of  the  colocated  system  to  two  squared  down  systems  which  were  squared 
down  as  follows: 
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With  Compensator 


Figure  9-11.  Characteristic  Gains  for  the  Compensated  Colocated  System 


POSITIVITY  6  (v) 


Table 

9-7.  Closed  Loop 

Poles  for  the 

Colocated  Desi 

PEAL  PART 

IMAGINARY 

PART 

DAMPING 

RATIO 

1) 

-8. 26864E-01 , 

4 . 61 434E-01 

.02072 

3) 

-4 . 6031  0E-01 , 

5.71 054F.-01 

.13652 

5) 

-5.49398E-01 , 

6 . 24499E-01 

.08010 

7) 

-3 . 84056E -01 , 

5. 99832E-01 

.1 5656 

9) 

1 . 01  249E-01  , 

4 . 74242E-01 

.46967 

7) 

1 .01  941  E  -  01  , 

4. 78232E-01 

.46536 

3) 

1 .48772E-01  , 

4.39533E-01 

.52509 

S) 

6.07551  E-01  . 

5.86374E-01 

.21521 

7) 

4.4181 7E-ol , 

5.87595E-01 

.31527 

9) 

6 . 67944  F - 01  , 

5.86729E-01 

.16104 

1) 

9. 37237E-01 , 

2 . 60093E-01 

.10188 

3) 

8 . 92957E-01  , 

3. 74850E-01 

.08042 

5) 

8.80506E-01  , 

3 . 04850E-01 

20737 

7) 

8 . 9061  1 E-01  , 

2.64750E-01 

.24655 

9) 

8',  62547  E-01 , 

2. 83442E-01 

.29107 

1) 

8. 65940E-01  , 

2. 78600E-01 

.29104 

3) 

8 . 65500E-01 , 

2 , 79053E-01 

.29137 

5) 

8 . 64333E-01  , 

2.80603E-01 

.29160 
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down  matrix  is  given  by  a(Z)  KQ.  KQ  is  given  in  Table  9-8(i) 
and  the  digital  filter  is  given  by  v 

a(Z)  =  1-82-10  3  Z3  +  5.46-10'3  Z 2  +  5.43-10"3  Z  +  1.82-10"3 
Z3  -  2.267  Z2  +  1.713  Z  -  4.317 

(ii)  Partial  Square  Down  Matrix  -  The  square  down  matrix  is  given 

by 


Kq  =  [I  [  Lq]  as  given  in  Table  9-8 ( i i ) 

As  it  can  be  observed  in  Figure  9-13  the  gains  of  the  two  first 
modes  are  increased  relative  to  the  other  modes  when  the  square  down 
procedure  is  applied.  This  is  true  both  in  the  case  of  a  limited  gain 
matrix  and  the  case  of  a  full  gain  matrix  +  third-order  filter.  The 
filter  gives  more  weight  to  the  important  modes  (modes  1  and  2),  however, 
it  introduces  lag  into  the  system  which  complicates  the  design  of  a 
control . 

It  was  found  that  the  squared  down  noncclocated  system  can  have 
its  poles  acquire  the  needed  amount  of  damping  by  using  a  constant 
diagonal  compensator.  This  is  true  for  both  the  partial  square  down 
matrix  Kq  =  [I  Lq]  and  the  full  matrix  Kq.  Since  the  partial  matrix 
results  in  a  more  simple  control  system,  this  design  will  be  explained 
and  compared  to  the  colocated  design.  The  square  down  matrix  is  as  in 
Table  9-8 ( i i )  and  a  compensator  4’I  was  included  to  improve  performance. 
The  resultant  complete  control  1 eruin  Table  9-8(iii).  The  functional 
diagram  of  the  noncolocated  control  system  is  given  in  Figure  9-14. 

Since  a  constant  gain  controller  yields  a  closed  loop  system  with 
adequately  damped  poles,  no  additional  filtering  is  necessary.  This 
simplification  of  the  controller's  structure  is  due  to  the  squaring 
down  matrix,  it  amplified  the  representation  of  important  modes  in  the 
feedback  signal.  This,  in  turn,  allowed  for  damping  augmentation  with 
a  constant  gain  without  destabilizing  the  high  frequency  modes .  Table  9-9 
shows  the  closed  loop  poles  of  the  noncolocated  system  with  the  controller 
as  given  in  Table  9-8 ( i i i ) .  As  it  can  be  observed,  both  the  nominal 
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Figure  9-14.  Functional  Diagram  of  the  Noncolocated  Control  System 


Table  9-9.  Noncolocated  Design  for  the  Tetrahedron:  Closed  Loop  Eigenval 


and  perturbed  models  acquire  considerable  modal  damping.  It  is 
especially  important  that  modes  1  and  2  have  high  damping  coefficients 
since  they  contribute  the  most  to  the  LOS  motion.  It  should  be  noted, 
a  sensitivity  study  showed  that  the  closed  loop  modal  damping  character¬ 
istics  are  not  sensitive  to  the  compensator's  gain:  a  change  of  the 
nominal  gain  (4*1)  by  6  db  gave  damping  within  specifications.  See 
also  Section  9.2.4. 

Figures  9-15  and  9-16  show  the  characteristic  gains,  the  positivity 
index  and  the  SVD  analysis  for  the  compensated  noncolocated  system.  The 
conclusions  one  can  draw  from  the  comparison  of  the  colocated  design 
given  in  Figures  9-11  and  9-12  and  the  noncolocated  design  shown  in 
Figures  9-15  and  9-16,  reinforces  the  results  obtained  for  the  plate 
exampl  e: 

(i)  The  colocated  and  noncolocated  systems  will  operate  at 
approximately  the  same  DC  gain. 

(ii)  Both  systems  can  tolerate  moderate,  low  frequency  modal 
uncertainties. 

(iii)  Both  systems  can  withstand  large  high  frequency  modal 
uncertainties. 

(iv)  The  colocated  system  has  better  robustness  properties  at 

low  frequencies.  The  situation  is  reversed  in  mid  frequencies. 

9.2.4  Simulation  Study  Results  for  the  Tetrahedral  Truss 

A  simulation  study,  similar  to  the  flat  plate  study  was  conducted. 
The  study  evaluated  the  performance  of  the  colocated  and  noncolocated 
designs  for  both  the  theoretical  and  perturbed  models.  The  evaluation 
was  done  by  studying  the  time  response  to  the  modal  initial  conditions 
given  in  Table  9-6  (see  Section  9.2.1).  The  response  of  the  nominal 
tetrahedral  model  for  both  the  colocated  and  noncolocated  designs  are 
shown  in  Figure  9-17.  As  it  can  be  seen  the  LOS  transient  requirements 
are  u^ed  for  both  designs,  however,  the  noncolocated  design  has  a 
superior  performance.  As  can  be  observed,  the  LOS  displacement  for  the 
noncolocated  design  in  both  the  X  and  Y  directions  is  within  the 
required  band  in  about  6  seconds,  which  is  in  about  30%  of  the  required 
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Figure  9-16.  Noncolocated  Design  for  the  Tetrahedron:  Positivity 
and  Singular  Value  Analysis  of  the  Return  Difference 
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time.  The  transient  of  the  tetrahedron  controlled  by  the  colocated 
system  reaches  the  required  band  only  after  about  15  seconds.  This 
could  have  been  expected  since  the  noncolocated  system  has  more  infor¬ 
mation  in  the  feedback  signal,  namely  the  LOS  x  and  y  motions  are 
measured  directly  by  the  two  additional  sensors. 

Figure  9-18  shows  the  response  of  the  nominal  model's  residual 
modes  (modes  not  included  in  the  design)  and  the  control  effort  of 
the  actua  -or  at  leg  #1 .  It  can  be  seen  that  the  actuator  command 
forces  ai  e  modest  in  both  cases,  with  the  noncolocated  design  having 
a  lower  control  force  by  about  40%.  The  residual  modes  are  only 
slightly  excited  by  both  control  configurations  and  they  converge  to 
zero.  The  noncolocated  control  sytem  excites  the  residual  modes  less, 
but  a  high  frequency  mode  is  lightly  damped. 

Figure  9-19  shows  the  LOS  and  transient  response  of  the  colocated 
and  noncolocated  designs  for  the  perturbed  model.  Both  designs  meet 
specification,  but  the  noncolocated  system  has  the  better  performance. 

The  noncolocated  system's  transient  response  settles  within  the  required 
band  in  about  9.5  seconds  compared  to  .'about  20  seconds  for  the  colocated 
system.  Hence,  as  predicted  by  the  analysis,  the  noncolocated  design 
is  slightly  more  robust  than  the  colocated  system. 

Figure  9-20  shows  the  response  of  the  residual  modes  and  the 
control  effort  of  the  perturbed  tetrahedron's  model.  It  can  be  seen 
that  the  control  effort  in  the  case  of  the  noncolocated  controller  is 
lower  than  in  the  case  of  the  colocated  system.  However,  the  colocated 
system's  residual  modes  are  better  damped  than  those  of  the  noncolocated 
controller.  This  was  anticipated  by  the  analysis  which  showed  more 
uniform  damping  at  high  frequencies. 

A  sensitivity  study  was  performed  to  determine  if  an  increase  of 
the  controller's  gains  causes  a  deterioration  in  system's  performance. 

It  was  found  that  the  nominal  controller's  gain,  4.1,  can  be  increased 
to  about  12.1  without  deterioration  of  the  performance.  Actually,  for 
certain  higher  values  of  controller's  gain,  the  system's  performance 
improves.  . igure  9-21  compares  the  LOS  motions  of  the  nominal  tetrahedron 
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Figure  9-18.  Nominal  Model:  Residual  Modes'  Response  and  Control  Effort 
for  Colocated  and  Noncolocated  Designs 


Control  Systems  for  the  Tetrahedron 


model,  with  the  nominal  controller  4*1,  and  with  controller  8*1.  I 
can  be  seen  that  the  performance  of  the  system  with  the  higher  gain 
is  better  and  the  LOS  transient  motion  reaches  the  required  band  in 
less  than  3  seconds. 

9.3  Conclusions  of  the  Design  Examples  Study 

Section  9  described  the  application  of  the  positivity  theory  and 
multivariable  frequency  domain  techniques  to  the  design  of  noncolocated 
control  systems.  The  study  showed  that  the  developed  squaring  down 
procedure  gives  good  results  and  can  be  incorporated  into  the  positivity/ 
frequency  domain  methodology. 

Comparison  of  colocated  and  noncolocated  systems  led  to  the  following 
concl usions : 

(i)  The  noncolocated  design  provides  better  damping  of  the 

important  modes,  i.e.  modes  which  influence  the  performance 
the  most. 

(ii)  The  colocated  design  provides  more  uniform  damping  of  the 
modes  and  therefore,  modes  not  included  in  the  design  model 
are  better  damped. 

(iii)  The  colocated  and  noncolocated  designs  for  the  examples 

presented  have  similar  robustness  properties  with  a  slight 
advantage  to  the  noncolocated  design.  However,  this  is  not 
a  general  conclusion  and  is  specific  to  the  examples  at  hand. 

(iv)  The  computational  requirements  for  both  systems  are  comparable. 
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APPENDIX  A.  TEST  AND  FLIGHT  EXPERIMENT  PROGRAM  PLAN 


The  objectives  of  the  test  plan  are  to  determine  how  and  in  what 
sequence  the  experiments  are  to  be  conducted.  This  includes  determination 
of  how  the  experiments  are  to  be  controlled  and  what  kind  of  on-orbit 
operations  are  required.  The  test  plan  also  addresses  the  testing  and 
correlation  of  analytical  models  with  the  post-flight  experimental  data. 

A . 1  Baseline  Experiment  Sequence 

The  plan  below  describes  how  and  in  what  sequence  the  experiments  are 
to  be  executed.  Referring  to  the  timeline  in  Figure  A-l,  it  will  take 
approximately  two  hours  to  sequence  through  the  experiments,  including 
triply  redundant  data  gathering. 

POC  Structure  Dynamical  Interaction  with  Orbiter-UQ  Seconds 

Function:  1)  Validate  analytical  method  used  to  moc_,  the  dynamic 
interaction  between  the  POC  structure  a, id  the  shuttle. 

2)  Validate  the  shuttle  environment  model. 


Operations:  1)  Monitor  the  temperature  of  structure,  the  shuttle 
DAP,  the  inertial/base  sensor,  and  control  system 
sensors. 


2)  Record  POC  structural  motion. 

3)  Stimulate  structure  with  pulses;  stimulate  with 
frequency  sweep;  stimulate  with  disturbances; 
threshold  against  ground;  record  POC  structural 
motion,  base  sensor  output  and  DAP  commands. 

4)  Apply  maximum  likelihood  modal  parameter  estimation 
technique  and  crosscheck  with  FFT. 


Optical  Element  Alignment-420  Seconds 

Function:  1)  Center  and  DC  calibrate  ranging  sensors. 

2)  Establish  "Disturbance  Free"  behavior. 


Operations:  1)  Use  six  mirror  alignment  actuators  to  center  laser 
diodes  in  the  FOV  of  the  SAMS  and  CCD  sensor. 


2)  Monitor  the  temperature  of  the  structure,  the 
Shuttle  DAP  and  the  inertial/base  sensors. 
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Figure  A-l.  Baseline  Experiment  Sequence  Timeline 


3)  Individually  perturb  each  mirror  alignment  actuator 
through  10  DC  levels.  Low  pass  and  time  average  to 
establish  linearity  behavior.  Threshold  against 
ground  values. 

4)  Observe  behavior  of  mirrors  and  structure  as  Shuttle 
limit  cycles  (3  cycles). 

LOS  Centroid  Alignment-260  Seconds 


Function:  1)  Calibrate  image  motion  sensing  system. 

2)  Calibrate  optical  element  motion  to  LOS  motion 
transformation. 

Operation:  1)  Acquire  bright  star  and  center  by  translating  CCD 
sensor. 


2)  Monitor  the  temperature  of  the  structure,  the  shuttle 
DAP  and  the  fnertial/base  sensors. 

3)  Individually  perturb  each  optical  alignment  actuator 
through  6  DC  levels.  Low  pass  and  time  average 
centroid  motion.  Threshold  against  ground  values. 


Structural  Control  Actuator  arid  Sensor  Alignment-180  Seconds 
Function:  Calibrate  actuator  to  sensor  transformation 


Operations:  1)  Fix  on  star  acquired  earlier. 

2)  Monitor  the  temperature  of  the  structure,  the 
Shuttle  DAP,  the  inertial/base  sensors,  LOS  motion 
and  control  sensor  signals. 

3)  Individually  step  each  structural  control  actuator 
through  6  DC  levels.  Low  pass  and  time  average  all 
sensor  signals.  Subtract  Shuttle  motion  to  get 

DC  gains.  Threshold  against  ground  values. 

Modal  Synthesis  Experiment-210  Seconds 

Function:  Validate  model  synthesis  as  a  technique  for  modeling  LSS 
from  smaller  substructures, 


Operations:  1)  Monitor  the  temperature  of  the  structure,  the  Shuttle  DAP,  the 
inertial/base  sensors,  mirror  motion  and  control 
system  sensors. 


2)  Stimulate  structure  with  pulses  and  store  response 
for  ground  processing. 
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3)  Stimulate  structure  with  sinusoidal  frequency  sweep 
and  mark  peaks.  Individually  stimulate  with  sinusoids 
at  critical  frequencies.  Store  free-vibration  decay 
for  ground  processing. 

4)  Apply  maximum  likelihood  (ML)  modal  parameter  estimation 
technique  and  crosscheck  against  free-vibration 

decay.  Apply  ML,  physical  parameter  estimation  and 
compare  with  assumed  physical  values. 

Disturbance  Recording  Validation-180  Seconds 

Function:  Verify  that  disturbance  "Recording"  is  representati ve 

of  expected  disturbances. 

Operations:  1)  Monitor  the  temperature  of  the  structure,  the 

shuttle  DAP,  the  inertial/base  sensors,  and  control 
system  sensors. 

2)  Turn  on  each  disturbance  actuator  (shaker)  individually. 
Cross-correl ate  with  ground  value  and  threshold. 

Lurcpgd  Passive  Element  Modeling  and  Control  Experiment-960  Seconds 

Function:  Validate  analytical  methods  used  to  model  structures  with 

integrated,  lumped,  passive  elements. 

Operations:  1)  Fix  on  star  acquired  earlier. 

2)  Monitor  the  temperature  of  the  structure,  the 
Shuttle  DAP  motion,  the  inertial/base  sensors, 
control  system  sensors,  and  LOS  sensor. 

3)  Turn  on  disturbances  and  record  response. 

4)  Individually  unlock  each  passive  isolator,  then  all 
simultaneously.  Record  response. 

5)  Lock  all  isolators. 

6)  Turn  off  disturbances. 

7)  Unlock  all  passive  elements.  Stimulate  with  pulse. 
Stimulate  with  sinusoidal  frequency  sweep  and  mark 
peaks.  Stimulate  with  sinusoids  at  critical  modes, 
record  free  vibration  decay  and  threshold  aaainst 
ground  values. 

8)  Turn  on  disturbances  and  record  data. 

9)  Unlock  all  passive  elements  and  record  disturbance  respons 
Threshold  against  ground  values. 

10)  Same  as  step  4  in  Modal  Synthesis  experiment. 
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Active  Structural  Control  System  Model inq  Experiment-960  Seconds 


Function:  Validate  analytical  methods  used  to  model  active  control 
Operations:  1)  Fix  on  star  acquired  earlier. 

2)  Monitor  the  temperature  of  structure,  the  Shuttle 
DAP,  the  inertial/base  sensors,  control  system 
sensors,  and  LOS  sensor, 

3)  Individually  turn  on  insensitive  inertial  and 
relative  control  systems.  For  each  one:  measure 
diagonal  gain  and  phase  margin,  and  threshold  against 
ground  values;  stimulate  with  pulse;  stimulate  with 
sinusoidal  frequency  sweep  and  mark  peaks;  stimulate 
with  sinusoids  at  critical  modes,  record  free- 
vibration  decay  and  threshold  against  ground  values. 

4)  Turn  on  both  inertial  and  relative  control  systems. 
Repeat  procedure  in  step  3. 

5)  Turn  on  disturbances. 

6)  Individually  turn  on  inertial  and  relative  control 
system.  Measure  LOS  error  and  figure. 

7)  Turn  on  both  inertial  and  relative  control  systems. 
Measure  LOS  and  figure. 

8)  Same  as  step  4  in  Modal  Synthesis  experiment. 


Advanced  Control  Law  Demonstration-780  Seconds 


Function:  Determine  stability  and  performance  limitations  of  the 
two  most  promising  control  laws. 

Operations:  1)  Calibrate  SAMS  and  LOS  motion  sensors. 

2)  Acquire  bright  star  and  center  unto  CCD  array. 

3)  Monitor  the  temperature  of  structure,  the  Shuttle 
DAP,  the  inertial/base  sensor,  control  system 
sensors,  and  LOS  sensor. 

4)  For  each  of  two  control  laws:  measure  diagonal 
gain  and  phase  margin,  and  threshold  against  ground 
values;  stimulate  with  pulse;  stimulate  with 
sinusoidal  frequency  sweep  and  mark  peaks;  stimu^fe 
with  sinusoids  at  critical  modes,  record  free 
vibration  decay  and  threshold  against  ground  values. 


5)  Turn  on  disturbances. 

6)  Record  LOS  and  mirror  motions  and  threshold  against 
ground  values. 

7)  Same  as  step  4  in  Modal  Synthesis  experiment. 

Have  Propagation  Experiment-180  Seconds 


Function:  Determine  performance  of  high  bandwidth  autoalignment  system. 
Operations:  1)  Fix  on  star  acquired  earlier. 

2)  Monitor  the  temperature  of  the  structure,  the  Shuttle 
DAP,  i nerti al /base  sensors,  control  system  sensors 
and  LOS  sensor. 


3)  Turn  on  high  bandwidth  autoalignment  system  and 
record  nominal  figure  and  LOS  motion.  Threshold 
against  ground  values. 

4) .  Turn  on  disturbances  and  record  figure  and  LOS  motion. 

Threshold  against  ground  values. 

Imaging  Experiment-600  Seconds 


Function:  Establish  performance  limitations  and  relative  merits  of 
two  competing  technologies. 


Operations:  1)  Point  telescope  at  space, 

2)  Monitor  the  temperature  of  structure,  the  Shuttle 
DAP,  the  inertial/base  sensor,  control  system 
sensors  and  LOS  sensor. 


3)  For  each  of  two  competing  technologies:  store  image 
for  one  Shuttle  limit  cycle;  turn  on  disturbances; 
and  store  image  for  one  Shuttle  limit  cycle. 

4)  Maneuver  Shuttle  to  point  telescope  at  earth. 

5)  Repeat  step  3. 


The  baselined  experiment  sequence  and  timeline  will  not  require  any 
unusual  time  slot  requirements  as  the  total  experiment  lasts  111  minutes. 
Seventy-eight  continuous  minutes  one  day  and  33  continuous  minutes  three* 
days  later.  The  primary  reason  for  the  break  is  to  a  How  for  ground 
validation  of  structural  modeling,  low  gain  feedback  behavior  and  predicted 
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subsystem  performance  prior  to  conducting  the  high  technology  demonstration 
The  approach  ensures  that  all  the  bugs  are  out  and  that  the  demonstration 
is  successful  the  first  time. 


A.  2  Control  of  the  Experiments 

The  baseline  test  structure  fits  comfortably  in  the  Shuttle  cargo 
bay  and  its  dynamic  interaction  with  the  Orbiter  is  immeasurable  to 
the  Shuttle  attitude  control  system.  For  this  reason  the  baseline 
design  also  does  not  require  any  deployment/erection  activities,  extra 
vehicular  activity  (EVA),  use  of  the  remote  manipulator  arm,  restowage, 
nor  an  emergency  jettison  procedure.  The  role  of  the  astronauts  during 
the  experiments  will  therefore  be  relatively  minor:  the  structure  fits  into 
the  shuttle  cargo  bay  so  it  needs  no  deployment  assistance  nor  astronaut 
extra  vehicular  activity;  the  CDHS  system  baselined  has  an  automated 
experiment  sequencer  to  prepare,  execute,  monitor,  and  step  through  the 
experiments;  the  structure's  caging  mechanism  is  electromechanical  and  is 
controlled  by  the  sequencer;  and  the  subsystems  failure  detection 
and  shutoff,  data  validation,  and  transmission  of  data  are  all  accounted 
for  in  the  automated  CDHS  architecture.  The  only  roles,  therefore,  that 
one  can  currently  baseline  for  the  astronauts  are  to  monitor  the  execution 
of  the  experiments  as  they  are  cycled;  to  provide  real-time  data  entry/ 
modification  if  failure  is  detected;  and  to  point  the  shuttle  cargo  bay 
away  and  towards  the  earth  during  the  imaging  experiments. 

A. 3  Concept  for  Correlating  Analytical  Models  with  Test  Data 

One  of  the  reasons  for  performing  the  experiments  is  to  learn  how 
to  model  structures  more  accurately.  There  are  two  immediate  goals  in 
mind.  One  is  to  determine  the  relationship  between  ground  and  on-orbit 
structural  behavior.  The  second  is  to  learn  how  to  combine  the  dynamics  of 
components  to  create  the  dynamics  of  the  total  system.  The  questions 
that  must  be  answered  in  order  to  do  this  have  both  a  functional  and  a 
parametric  nature: 


Does  the  open-loop  composed  system  behave  as  predicted  by  the 
modal  synthesis  procedure? 

Is  the  linear  feedback  model  (almost  always  assumed)  a  reasonably 
good  predictor  of  closed-loop  system  behavior?  Or,  does  the 
feedback  enter  in  a  more  complicated  manner? 

For  a  given  system  model  functional,  does  the  parameter  which 
creates  the  best  data  fit  equal  the  physical  values  in  the 
sub-components? 

It  is  necessary  to  determine  that  the  answer  to  the  posed  questions 
is  affirmative,  even  for  micron  and  nanometer  motion. 

Figure  A-2  conceptualizes  a  possible  approach  for  finding  the  answers. 
The  idea  is  to  validate  the  modal  synthesis  procedure  by  comparing  its 
form  with  the  global  FEM  model.  Both  models  are  then  validated  by  showing 
that  the  parameters  providing  the  best  data  fit  are  realistic  physical 
values.  The  analytical  models  will  be  constructed  from  the  material 
properties  of  the  components,  the  actuator  and  sensor  dynamics,  and  the 
control  laws  (things  one  can  model  well).  The  model  construction/ 
functional  is  initially  assumed  to  be  linear  then  adjusted  along  with  the 
parameterization  until  both  the  model  predictions  and  the  parameters  match 
the  true  values. 
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Figure  A-2.  Concept  for  Correlating  Analytical  Models  with  Test  Data 


APPENDIX  B.  EXPERIMENT  MILESTONES,  SCHEDULE  AND  COST 


The  final  objective  of  the  PREMOD  contract  was  to  outline  the  mile¬ 
stones,  schedule  and  rough-order-of-magni tude  (ROM)  cost  estimates 
required  for  the  government  to  accomplish  experiment  objectives. 

B. 1  Milestones  and  Schedule 

Figure  B-l  summarizes  the  experiment  milestones  and  schedule.  The 
total  program  is  baselined  to  last  four  years.  The  first  year  will  be 
spent  primarily  on  subsystem  and  component  design  work,  planning  subsystem 
tests  and  integration,  and  on  specifying  the  experiment  tests  and  test 
conditions.  The  second  year  will  be  used  to  build  the  hardware,  program 
and  debug  software,  test  and  integrate  the  subsystem  components  and  on 
planning  the  system  and  Shuttle  integration  and  test.  The  first  half  of 
the  third  year  will  be  used  to  complete  the  integration  of  the  ground 
experiment  and  the  second  half  will  be  used  to  execute  and  process  the 
ground  experiments.  The  fourth  year  will  be  spent  integrating  the 
experiment  into  the  Shuttle,  with  a  planned  launch  and  on-orbit  experiment 
execution  following  in  the  third  quarter.  The  on-orbit  experiment  data 
processing  is  planned  for  the  fourth  quarter  of  the  fourth  year. 

B. 2  ROM  Costs 

The  roug,i-order-of-magnitude  cost  estimate  for  the  total  program 
is  shown  in  Table  B-i. 

Table  B-l.  ROM  Cost  Estimate 


SIMPLIFIED  GROUND 
EXPERIMENT 
( $ K  IN  1982) 


COMPLETE  GROUND  AND 
ON-ORBIT  EXPERIMENT 
(SK  IN  1982) 


COST  ITEM 


TOTAL  FOR  HARDWARE 
TOTAL  FOR  SHUTTLE  USE 


TOTAL  FOR  LABOR 


PROFIT  (15%) 

MARGIN  FOR  ERROR  (15%) 


BSE 


detailed  telescope  optics  design 

BUILO  4  INTEGRATE  TELESCOPE  OPTICS 
TEST  STRUCTURE  DESIGN 

BUILD  TEST  STRUCTURE  4  ASSOC. SUBSTRUCTURE 


Figure  B-l.  Experiment  Milestones  and  Schedule 


The  cost  breakdown  for  the  total  experiment  is  shown  in  the  Table  B-2 
and  is  approximately  $9M  for  hardware,  $14. 8M  for  Shuttle  costs,  $9M 
for  labor  costs,  S4.9M  profit,  and  S5.7H  (15%)  for  error  margin.  The 
total  cost  is  $43. 5M.  The  Shuttle  cost  estimate  is  detailed  further  in 
Table  B-3.  The  driving  cost  ($10. 3M)  is  primarily  due  to  the  cost  of 
standard  services  for  a  complete  pallet  flight,  (one  user  on  the  pallet, 
all  other  Shuttle  facilities  are  sha'ed).  The  key  assumptions  used  in 
making  the  estimate  are  the  accuracy  of  the  "Civilian  Reimbursement 
Guidelines1  document  dated  May  1980  Qo]  and  the  assumed  10%  inflation 
factor.  All  estimates  quoted  are  in  1982  dollars.  In  reality  this 
experiment  is  a  DOD  flight,  having  a  lower  negotiated  cost  with  NASA 
than  other  government  users.  Therefore,  the  real  cost  of  flight  may  be 
slightly  lower.  The  actual  figures  were  being  negotiated,  however,  by 
DOD  and  NASA  at  the  time  the  estimate  was  made. 


Table  B-2.  ROM  Cost  Estimate  Breakdown 


COST  ITEM 

SIMPLIFIED  GROUNC 
EXPERIMENT 
($K  in  1982) 

COMPLETE  GROUND  & 
ON-ORBIT  EXPERIMENT 
($K  in  1982) 

HARDWARE 

SAM  Sensor 

410 

410 

CCD  Sensor 

-- 

2510 

Experiment  Monitor 

90 

90 

CDHS  System 

155 

155 

Control  System 

291 

291 

Identification  System 

123 

123 

Test  Structure  &  Supporting 
Hardware 

115 

115 

Ground  Test  Equipment  Monitor 

330 

330 

Primary  Mirror  &  Figure  Control 

-- 

4500 

Secondary  Mirror  &  Assoc.  Optics 

-- 

500 

TOTAL  FOR  HARDWARE 

$  1 5 14K 

$9024K 

SHUTTLE  COSTS* 

Standard  Fees 

10247 

El ement  Use  Fees 

1789 

Optional  Services 

2798 

TOTAL  FOR  SHUTTLE  USE 

— 

$ 14834K 

LABOR  COSTS 

$3750 

9000 

SUBTOTAL 

S5264K 

S32858K 

PROFIT  (15%) 

790 

4929 

MARGIN  FOR  ERROR  (15%) 

908 

5668 

TOTAL 

S6962K 

S43455K 

*  Assumes  complete  pallet  flight,  civilian  reimbursement  guidelines 
dated  May  1980,  10%  inflation  factor. 
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Table  B-3.  Breakdown  of  Shuttle  Cost  Estimate* 


STANDARD  SERVICES  -  LOAD  FACTOR  =  0,23 

Shuttle  Operations 
Shuttle  Use  Fee 

Facilities  and  Equipment  Use  Fee 
Payload  Specialist 

TOTAL  COST  OF  STANDARD  SERVICES 

ELEMENT  SERVICES 

Operations  -  Complete  Pallet 
Element  Use  Fee  -  Complete  Pallet 

TOTAL  COST  OF  ELEMENT  SERVICES 

OPTIONAL  SERVICES 

Vandenberg  FI ight 
Unique  Integration  Requirements 
Additional  Resources 
Special  Access  to  Payload  Operation 
Special  Software  in  Command  &  Data 
Management  Subsystem 
Data  Transmission  Outside  of  JSC 
Payload  Control  Center 
Transmission  of  Space! ab  Data  Not 
Contained  in  the  STS  Operational 
Instrumentation  Telemetry  Link 
Payload  Specialist  and  Training 
Thermal  Loads  Analysis 
Structural  Loads  Analysis 
Additional  Days  of  STS  Support  (1  Day) 
Payload  Operations  Control  Center  in 
Most  Mode  -  Est 

Optional  Transportation  Services  -  Est 
Optional  Security  Services  -  Est 
Optional  Photographic  Services  -  Est 
Guaranteed  Launch  Date 

_ TOTAL  COST  OF  OPTIONAL  SERVICES 

TOTAL  SHUTTLE  LAUNCH  COST 


COST  $M 
1975 


$5 . 26M 


S0.92M 


S1.44M 


COST  $M,  1982 
10%  INFLATION/YR 


$10. 3M 


$  1 . 8M 


S2.8M 
$14. 8M 


"Space  Transportation  System  Reimbursement  Guide,"  May  1980  [10]. 


APPENDIX  C.  A  MODEL  OF  THE  SHUTTLE  DIGITAL  AUTOPILOT  (DAP) 


This  appendix  briefly  describes  the  DAP  model  implemented  in  the 
dynamic  interaction  study.  The  model  was  based  on  the  description  given 
in  [9]  with  up-to-date  parameters  obtained  from  Johnson  Space  Center.  In 
the  following,  explanation  of  the  various  model  subroutines  is  given  and 
some  of  the  logic  is  explained.  For  more  detailed  description  see  [9]. 

C. 1  Phase  Plane  Logic 

The  phase  plane  logic  implemented  in  the  Shuttle  Orbiter  DAP  is  based 
on  a  time  optimal  control  law.  It  assumes  low  coupling  between  the  axes 
and,  hence,  the  dynamics  of  the  system  around  each  axis  can  be  approximated 
as  double  integrator.  The  time  optimal  feedback  solution  for  this 
problem  is  well  known.  However,  the  optimal  control  law  cannot  be  im¬ 
plemented  precisely  because  of  the  following  imperfections:  attitude  and 
rate  measurement  errors,  measurement  disturbances  and  noise,  and  quantization 
errors  -  both  quantization  of  the  jet  firing  time  and  the  quantization  of  the 
attitude  and  rate  measurement  in  the  A/D  converter  from  IMU.  Hence, 
implementation  of  the  optimal  control  law  would  lead  to  continuous  limit 
cycle  and  to  fuel  waste.  To  avoid  these  problems,  the  phase  plane  logic 
was  designed  to  hit  tne  smallest  limit  cycle  in  one  pass  and  minimize 
fuel  waste.  Moreover,  the  logic  was  further  developed  to  take  into  account 
the  coupling  between  the  axes  when  the  system  is  in  the  hysteresis  region. 

Figure  C-l  shows  the  various  regions  in  the  phase  plane  of  the 
Orbiter’s  axis.  This  logic  is  implemented  for  each  control  axis,  with 
the  three  commands  forming  the  command  control  vector.  The  control 
vector  is  then  fed  into  the  jet  selection  logic  subroutine  to  choose  the 
jets  which  generate  the  control  torque. 


C.2  Vernier  Jet  Selection  Logic 

Since  the  Shuttle  was  assumed  to  be  in  its  normal  operation  mode, 
only  the  vernier  jet  selection  logic  was  simulated.  The  vernier  jet 
selection  logic  obtains  the  commanded  command  vector,  -AR,  from  the  phase 
plane  logic  module  and  computes  a  dot  product,  D.,  of  AR  with  the  angular 
increment  vector  a.  of  thruster  j  for  j=l, _ ,6.  See  Figure  C-2. 

J 

The  thruster  with  maximum  Dj  is  designated  and  is  fired.  If  the  next 
largest  D ^ ,  call  it  D2,  is  greater  than  T^  •  this  jet  is  fired. 

The  next  largest  D. ,  say  D3,  is  fired  if  D3  >  T2  .  D^.  Hence,  a  maximum 
of  three  vernier  thrusters  fire  at  one  time.  It  should  also  be  noted 
that  when  the  Shuttle  is  in  vernier  operation  mode,  the  element  of 
commanded  vector  AR  can  be  different  from  +  1. 


AR  -  COMMANDED  VERNIER  VECTOR 

a.  -  ATTITUDE  CHANGE  VECTOR  DUE 
TO  THRUSTER  j 


Figure  C-2.  A  Command  Vector  and  Angular  Increment  Directions 
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Table  C-l  summarizes  the  features  of  the  vernier  thrusters. 

The  firing  directions  of  various  vernier  thrusters  are  shown  in  Figure  C-3 
The  Orbiter  possesses  six  vernier  thrusters  which  produce  a  thrust  of 
~24  lbf  each.  Two  thrusters  are  located  in  the  front  and  generate  thrust 
in  the  xy  plane,  forming  an  angle  of  45°  with  the  x-axis.  The  rear 
thrusters  consist  of  two  pairs.  One  pair  fires  in  the  +z  direction  and 
the  other  in  the  +y  direction.  The  vernier  thrusters  are  not  designed 
to  produce  forces  in  -z  or  x  directions.  Thus,  any  altitude  change  or 
drag  makeup  has  to  be  performed  using  the  primary  thrusters. 


Figure  C-3.  Vernier  Thrusters  Firing  Directions 


Table  C-l.  Vernier  Jet  Locations  and  Thrust  Directions 
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C.3  State  Estimator 


The  state  estimator  consists  of  two  parallel  Kalman-tyne  filters 
one  for  acceleration  estimation  and  one  for  rate  estimation.  The 
S-plane  formulation  for  the  acceleration  estimate  is 


where 


Each  of  two  filters  is  divided  into  a  measurement  incorporation  section 
and  an  extrapolation  section.  Measurement  incorporation  for  both  filters 
is  performed  at  the  attitude  measurement  update  rate  (6.25  Hz)  by 
FILTER  1 . 

Extrapolation  for  both  filters  is  performed  at  12.5  Hz  by  FILTER2. 

In  addition,  FILTER1  provides  a  25-Hz  extrapolation  of  the  measured- 
attitude  variable  ATTITUDE.  This  variable  is  reset  to  a  pure  measurement 
value  by  FILTER1  at  6.25  Hz. 
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Abstract 

Because  of  the  cost  and  flexibility  advantages 
offered  by  microprocessors,  any  control  system  to 
be  considered  for  Urge  Space  Structures  (LSS)  will 
likely  be  digital.  This  short  paper  proposes  a 
combined  positivity  and  W— plane,  frequency  domain 
design  procedure  for  synthesizing  the  control  of 
these  structures.  Positivity  theory  is  used  to 
ensure  the  stability  and  stability  robustness  of 
digital  designs.  Inclusive  of  modal  truncation 
and  aliasing.  Shaping  the  characteristic  gains 
in  an  a-shifted  W-plane  synthesizes  the  required 
performance.  The  combined  approach  has  the 
advantage  that  it  leads  directly  to  robust  and 
implementable  low-order  controllers.  A  design 
example  using  an  LSS  like  structure  concludes 
the  paper. 


I.  Introduction 

Of  the  existing  control  problems,  the  active 
control  of  large  space  structures  (LSS)  perhaps 
most  clearly  demonstrates  the  difficulties  in¬ 
volved  in  designing  for  stability  and  performance 
robustness  while  restricting  oneself  to  using  only 
control  laws  of  implementable  complexity.  The 
LSS  "plant"  is  unique  in  the  sense  that  it  exhibits 
high  (theoretically  infinite)  state  dimensionality, 
very  low  damping,  only  poor  "plant"  models  are 
available,  and  the  disturbance  sec  is  widely  varied. 

>i  Briefly,  the  LSS's  dynamics  centered  about  the 
rigid"  body  motion  can  be  approximated  by  the  modal 
state  representation: 


where 


*i  "  Ai*i  +  Biu 


y  "  1  Cixi 

i»l  1  1 


(1) 


«  Ap 

i  L 


position  of  ith  mod 
rate  of  ith  mode 


Ai" 


-2Vi 


■] 


t  IR 


2x2 


input  influence  coeff.  for  ith  mode 


c  R 


2xm 


C  ■  |0  '  ouCPuC  rate  measurement  "1  rx2 
i  I  |  coefficients  for  ith  mode  111 

A 

m  ■  number  of  actuators 

A 

r  *  number  of  sensors 

and  n  is  theoretically  infinite  or,  more  practi¬ 
cally,  very  large  order  (>100  states).  The 
parameters  in  this  representation  are  generally 
computed  by  solving  the  eigenvalue  problem  re¬ 
sulting  from  finite  element  modelling  of  the 
structure  [1].  This  approach  is  currently  the 
best  practical  method  for  modelling  the  structure, 
yet  the  resulting  fidelity  of  the  parameters 
Bi  and  Ci  degrades  drastically  with  increasing  i 
(> 100Z  error  in  some  cases).  Independent  of  the 
finite  element  approach,  the  damping,  5.,  of  the 
structure  is  poorly  known  and  is  very  small 
(.001  -  .01).  Small  uncertainties  in  the  finite 
element  model  may  therefore  induce  very  large, 
non-monotonic  changes  in  the  input/output  map. 

The  phenomenon  occurs  even  when  one  is  careful 
and  uses  a  frequency  dependent  norm  which  restricts 
Che  bandwidth  of  the  input  signal  set.e.g.,  think  of 
Che  response  of  (1)  to  a  sinusoid  just  on  and  just 
off  an  assumed  resonant  frequency.  For  this  reason 
model  order  reduction  within  the  control  band¬ 
width,  based  only  on  open— loop  measures,  is  very 
difficult  and  a  large  order  plant  model  often  has 
to  be  retained.  The  impact  of  this  difficultv  be¬ 
comes  most  evident  whenever  one  attempts  to  base 
and  implement  robust  control  laws  whose  order  is 
fundamentally  tied  to  the  dimension  of  the  plant's 
state  space. 

The  disturbance  set  about  which  the  LSS  is  to 
be  regulated  is  non— trivial.  The  more  common  dis¬ 
turbances  constanc  (e.g.,  gravity  gradients)  and 
random  (e.g.,  equivalent  noise  of  reaction  wheel 
bearings)  are  encountered  together  with  impulsive 
(e.g.,  rapid  slew  maneuvers)  and  periodic*  (e.g., 
reaction  wheel  imbalance) .  This  wide  variety  of 
disturbances  makes  designs  employing  the  placement 
of  disturbance-to-error  transmission  zeros  poten¬ 
tially  too  sensitive  and  complex. 

In  order  to  obtain  robustness  with  an  imple¬ 
mentable  controller,  in  spite  of  the  difficulties 
represented  by  an  LSS,  this  paper  proposes  a  com¬ 
bined  posi tivity/frequency  domain  approach  to  syn¬ 
thesizing  digital  LSS  concrol  systems.  The 

*  not  necessarily  sinusoidal 
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characteristic  gains  of  the  LSS  are  shaped  in  an 
a-shifted  W-plane  In  order  to  synthesize  performance. 
Positivity  theory  is  used  to  ensure  the  stability 
robustness  cf  tbe  design.  The  paper  will  first 
summarize  the  main  results  from  positivity  and 
multi-variable  frequency-domain  design*  theory. 

These  results  are  then  interpreted  in  context  with 
the  LSS  plant.  This  is  followed  by  the  design 
procedure  and  concluded  with  a  design  example  using 
a  tetrahedral  truss. 

II.  Synopsis  of  Positivity  Theory 

First,  some  definitions  are  required. 

Let  f(t)  be  a  real,  possibly  vector  valued 
function  defined  on  [0,»).  The  truncation  of  f(t) 
at  t«T  is  defined  as 

I  f (t),  0  <_  t  <_  T 

fT(t)  |  0  ,  t  >  T 

Define  the  extended  Hilbert  Space  He  such  that 
f(t)  e  Hg  if  and  only  if  fT(t)  £  Hilbert  Space,  H. 

An  operator  with  a  domain  and  range  in  He  is 
then  defined  to  be  (strictly)  positive  if  for 
f(t)  £  He, 

<fT(t),  HfT(t)>  >_  0  (>0) 

Positive  oparators  are  useful  because  they 
have  the  following  properties : 

1)  The  inverse  of  a  positive  operator  is 
positive. 

2)  The  sum  of  positive  operators  is  positive. 

3)  Positivity  Theorem  [2,3,4]:  The  negative 
feedback  systems  shown  in  Figures  1  and  2 
are  bounded- input-bounded-output  BIBO 
stable  for  all  inputs  in  He  if  both  G  and 
H  are  positive  operators  and  at  least  one 
of  them  is  also  strictly  positive.  The  main 
advantage  of  the  theorem  is  that  conditions 
are  Imposed  Individually  on  G  and  H,  en¬ 
suring  stability  when  the  loop  la  closed  by 
negative  feedback. 


Figure  2.  Cascade/Unity  Feedback  Interconnection 

A  positive  operator  may  be  interpreted  physi¬ 
cally  as  one  that  does  not  create  energy  (a 
strictly  positive  operator  always  dissipates 
energy).  For  this  reason,  a  positive  operator  is 
often  referred  to  as  a  passive  operator. 

The  term  "Positive  Real"  is  also  often  used 
to  refer  to  a  positive  operator  chat  is  also  con¬ 
volutional  and  linear  time-invariant  (LTI) .  The 
term  "Hyperstability"  refers  to  the  Positivity 
Theorem  when  at  least  one  of  the  subsystems  G  or  H 
is  positive  real,  while  the  other  may  be  nonlinear 
and/or  time-varying. 

A  frequency  domain  test  exists  for  testing  the 
positivity  of  square  operators  which  have  a  Laplace 
or  Z-Transform.  representation  [4,5],  Positivity 
is  determined,  using  this  test,  by  first  computing 
the  positivity  index: 

[G(s“Jui)+G*(s“Ju) ]|  ,  ut[0,-) 
or  for  discrete  systems 

*  Xmin  {ffC(2-eJ“T)4G*(z-eJuT)]J  ,  UE[0,  *) 

where  {  •  }  denotes  "minimum-eigenvalue-of”  and 

denotes  complex-con jugate-transpose . 

If 

d(u)  >  0  for  all  u>  defined  above,  then  G  is 
strictly  positive  real 

i(w)  >_  0  for  all  u  defined  above,  then  G  is  posi¬ 
tive  real 

5(u)  <_  0  for  some  u  in  the  defined  range,  then  G 
is  not  positive. 

Finally,  there  is  a  strong  relationship  between 
the  Small-Gain  Theorem*,  the  Positivity  Theorem, 
and  the  concept  of  conicity  [2],  The  link  is  thn 
set  of  stability  invariant  transformations  called 
embedding.  As  Appendix  A  shows,  if  G  and  H  are 
LTI,  then  the  stability  of  che  svstem  in  Figure  1 
implies  the  stability  of  the  embedded  system  in 
Figure  3  and  vice-versa.  This  stability  invariant 
property  however,  does  not  in  any  way  imply  input- 
output  equivalence.  The  relationship  between 
small-gain,  positivity,  and  conicity  can  now  be 


*  TP j  system  in  Figures  1  and  2  are  BIBO  if  for 

I H  f  ■' 


*  We  restrict  ourselves  only  to  generalized 
Nyquist/Bode  methods. 


A  I'**  ‘•t'I 

g(H)  -  sup  -p- Ij  ,  g(H)  •  g(G)  <  1. 


Figure  3.  Example  of  Embedded 

Feedback  IncercoimecCion 


demonstrated.  Using  the  small  gain  theorem  on  the 
system  in  Figure  3  (with  D  2  0  and  Sal)  implies  if 

IN!  *  ||  (G-C)  R_1||  <  1  (2) 

and 

llHll  •  II  R  (I+HC)'1  H 1 1  <  1  (3) 

Then  the  system  in  Figure  1  is  stable*.  But  Eq .  (2) 
implies  M  '  ' 

IKg-c>  fTl|  <  Hr  fT||  v  fT 


This  aquation  may  be  geometrically 
cone  with  center  C  and  radius  R. 
constraints 


interpreted  as  a 
If  che  additional 


C  *  I  •  constant 
R  ■  I  •  constant 
C  -  R  -  0 
C  +  R  -»  I  •  - 


are  placed  on  Equation  (2),  then  che  identity  (real- 
inner  product  space  only;  y  (real- 

|(G-C)  fT  II  -  !|RfT;|  <Ge<G-C+R)fT,  (G-C-R)f p  <  Q 

is  equivalent  to  che  positivity  condition.  It  is 

-raanhCir\C  lv  houever.  Chat  had  we  demanded  chat 

j  nd  H  in  Eq.  (2)  and  (1)  satisfy  the  positivity 

"d  ‘  m  (i'*'  ’  eiCher  G  >  0  and  H  >  0,  or  H  >  0 
and  G  >  0),  che  stability  of  the  system  in  Figure  1 


cef  inarFiCgur:frati0n  “  3PPlied  C°  -X- 


could  still  be  inferred  even  chough  che 
can  no  longer  be  interpreted  as  a  cone. 


geometry 


HintT  infCerPretdti°"s  can  he  given  to  each  embed- 
g  transformation.  The  interpretation  depending 
on  .he  stability  theorem  with  which  it  is  used  * 
Some  interpretations  of  che  mathematics  include 
modeling  of  plant  uncertainty,  distance  from  pas¬ 
sive  systems ,  bounding  of  unmodeled  dynamics,  and 
che  model-reference  interpretation  given  in  [?]. 

HI  •  Summary  of  Multivariable 
Frequency-Domain  Design 

^  objective  in  frequency  domain  design  is  to 
meet  and/or  optimize  performance  goals  under  the 
constraint  chat  the  closed-loop  system  be  stable. 

The  Key  to  extending  chis  classical  Nvquist/Bode 
approach  to  che  multivariable  case  [8]  are  '•he 
characteristic  gains  of  che  square  cransfer ’macrix 
G(s),  describing  the  open-loop  plant.  The  charac¬ 
ter  scic  gains  X(s)  are  che  eigenvalues  of  G(s)  as 
function  of  frequency,  i.e.,  che  solution  to 


G(s)  v(s)  -  \(s)  v (s )  (4) 

The  importance  of  X  ( s )  stems  from  the  fact  chat  the 
Xc(s)  of  a  closed-loop  unity  feedback  svscem  are 
related  to  the  open-loop  \0(s)  in  a  manner  chat  is 
analogous  to  single- inpuc-single-oucpuc  systems. 


closed-loop  gain  »  open- loop  gain 

1  +  open-loop  gain 

Using  *nd  shaping  the  open- loop  X (s )  and  v(s): 

1)  Multivariable  Nyquisc  Theorem  [3].  A 
closed-loop  system  is  stable  if  and  only 
if  the  net  sum  of  counter-clockwise  encir¬ 
clements  of  che  (-1. , 0.  )  point  by  the 
characteristic  gains  are  equal  to  the  num¬ 
ber  of  open-loop  unstable  poles. 

2)  With  prooer  interpretation,  any  closed  con¬ 
tour  can  be  taken  while  applying  the  multi¬ 
variable  Nyquisc  Theorem.  The  basic  con¬ 
cept  one  is  exploiting  in  the  Nvquist 
Theorem  is  che  Principal  of  the  Argument  5], 
This  means  that  one  can  eiche-  conformally 
map  any  iesired  closed  contour  into  toe 
Nvquist  Path  and  then  evaluate' Eq.  (4),  or 
evaluate  (-)  along  the  desired  closed  con¬ 
tour  directly.  The  Nvquist  Theorem  then 
gives  conditions  under  which  no  closed-iooo 
poles  will  be  in  the  region  represented  by 
the  closed  contour.  Figure  u  snows  two  of 
Che  possible  mappings  which  can  be  used  Co 
ensure  that  satisfaction  of  che  Theorem  will 
imply  chat  ciosed-loop  transients  exponential! 
decay  faster  than  some  desired  value  a. 

3)  To  achieve  low  control  interaction  in  a 
frequency  band  it  is  sufticienc  to  ensure 
either  that  ail  the  open-loop  X0(s)  have  a 
large  magnitude  over  that  band,  or,  that 
all  the  open-looo  v(s)  are  nearly  ortho¬ 
gonal  over  that  band. 


4)  In  frequency  bands  of  low  interaction,  the 
notion  of  gain  margin  and  phase  margin  may 
be  applied  as  a  qualitative  assessment  of 
performance.  Furthermore,  in  these  bands, 
the  plot  of  the  open-loop  V^s)  in  the 
Nichols  chart  can  provide  a  direct  measure 
of  sensitivity,  tracking  accuracy,  and  dis¬ 
turbance  rejection  along  directions  com¬ 
putable  from  v(s)*. 

IV.  LSS  Plant  from  a  Positivity  and  Frequency- 
Domain  Perspective 

There  are  several  reasons  for  approaching  the 
LSS  controller  design  problem  from  a  Positivity/ 
Frequency-Domain  perspective.  Foremost  is  the  fact 
that  an  LSS  with  collocated,  ideal  (infinite  band¬ 
width)  rate  sensors  and  actuate  :s  is  a  positive 
operator  [4],  This  result  is  Independent  of  the 
parameters  characterizing  the  LSS.  From  equation 
(1),  the  transfer  function  of  the  idealized  LSS 
also  has  reciprocity  properties  (i.e.  in  modal 
coordinates  B~«C,  so  the  transfer  function  from 
input  i  to  output  j  is  the  same  as  from  input  J 
to  output  i).  As  a  consequence  of  positivity  and 
reciprocity  under  ideal  conditions,  tha  LSS 
transfer  function  with  collocated  finite  bandwidth 
actuators  and  rate  sensors  exhibits  nearly  ortho¬ 
gonal  eigenvectors  (v(s)  in  equation  4)  and 
positivity  over  most  of  the  actuator/sensor 
bandwidth.  For  small  negative  excursions  of  4  (u ) 
or  when  the  actuators  and  sensors  begin  to  add 
significant  phase  shift,  the  Positivity  Theorem  is 
applied  with  embedding  [5]. 

It  is  also  reasonable  to  use  conlcity  con¬ 
cepts  in  treating  the  LSS  problem  but  there  are 
some  difficulties.  The  LSS  plant  has  closely 
spaced  modes  with  almost  no  damping.  This  makes 
the  singular  values  of  G(s)  behave  violently  at 
the  resonant  frequencies.  Also,  while  conlcity  has 


in  general  the  potential  for  providing  a  tighter 
bound  (assuming  embedding  is  not  permitted  in  the 
Poaitivity  Theorem),  in  practice,  the  uncertainty 
in  LSS  modal  frequencies  is  so  high  that  the  cone 
parameters  nearly  degenerate  to  the  positivity 
conditions.  In  other  words,  the  LSS  undergoes 
180*  phase  shifts  (-40  dB  to  +40  dB  with  £-.005) 
a  large  number  of  times  and  at  closely  spaced  fre¬ 
quencies  whose  center  location  is  very  uncertain. 
The  uncertainty  of  the  center  and  behavior  bound 
are  therefore  modeled  by  conlcity  parameters  which 
are  practically  equal  to  the  positivity  conditions. 

It  therefore  makes  3ense  to  analyze  LSS 
behavior  in  terms  of  its  proximity  to  a  positive 
system.  This  is  measured  by  the  positivity  index 
4(w)  in  the  sense  that 


G  -  4(w)  •  I 

Is  a  positive  operator.  Hence  5(u) 


1  Is  an 


appropriate  "C"  operator  in  a  positivity  inter¬ 
preted  embedding  procedure.  The^implied  Theorem 
conditions  which  the  controller  H  >  0  satisfies 
(R-I;  S-I;  D-0)  is  that  it  can  phase  stabilize  all 
frequencies  where  4 (w)  >  0,  but,  must  gain  sta¬ 
bilize  wherever  4(w)  <  0.  , 


In  general,  any  skew 

Hermit ian “matrix  can  also  be  added  to  H  with¬ 

out  affecting  stability. 
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Figure  4,  Linear  MdO  Frequency  Domain  Design:  Specified  Degree  of  Stability  Via  a-Shift 


*  It  is  acknowledged  that  sensitivitv  and  disturbance 
rejection  are  generally  best  measured  using  singu¬ 
lar-value  analysis,  but  the  LSS  plant,  as  will  be 
shown  later,  has  some  very  interesting  properties. 


V.  Recipe  for  Sample-Data  LSS  Control 
Design  In  che  W-plane 

1.  Compute  the  characteristic  gains  X (w)  and  the 
positivity  Index  6  (u)  for  the  open-loop,  sampled 
plant,  i.e.,  solve 


and 

5  (w)  -  X 


DET  [G<«)]>  -  X(w)  •  i}  -  0  (5) 

|{'a/|^[G(s)]|+?i/{^[G(s)]f}J  (6) 


where 


^[•]  denotes  "take  the  z-transform  of  [•]" 
Ic/i'}  denotea  "substitute  for  z  In  {•}" 

*  denotes  "take  the  complex  conjugate 
transpose" 

X_.  {•}  denotes  "take  the  minimum  eigen- 
value  of 

and  X (w)  and  6 (w)  is  solved  for  with 
w  -  jv,  v  E [0,») 

2.  Shape  X(w)  to  meet  sensitivity  and  disturbance 
rejection  requirements  in  the  bandwidth  of 
interest*.  Non-mudal  sensitivity  and  disturbance 
rejection  are  approximated  by 


1  +X3(w) 


X 


3(w)  -  "shaped"  X (w) 


6 (w)  can  be  advantageously  used  here  to  indicate 
the  class  of  shaping  compensators  which  will 
preserve  stability  Inclusive  of  modal  errors 
and  any  model  order  reduction  which  may  have 
been  used.  I.e.,  5 (w)  can  be  used  to  measure 
the  permissible  low  frequency  multiplicative 
perturbstions  and  high  frequency  additive  per¬ 
turbations.  The  minimum  of  a (w)  and  5- (w) 

[the  positivity  index  for  the  worst  case  per¬ 
turbed  system)  as  function  of  w  bounds  the 
gains  of  the  permissible  positive  filters/con¬ 
trollers. 

Note:  5  is  never  really  computed.  It  is  esti¬ 
mated,  like  the  conicity  parameters,  from 
physical  information  like  range  in  damping 
ratio,  modal  data  uncertainty,  and  actuator/ 
sensor  bsndwidth  variation. 


If  the  desired  gains  for  Xg(w)  cannot  be 
attained  together  with  stability  and  a  $(w) 
indicative  that  the  worst  case  plant  can  be 
encompassed,  then  a  solution  to  the  problem 
posed  is  not  possible.  One  must  then  trade 
reductions  of  performance  and/or  system  uncer¬ 
tainty  with  increases  in  actuator/sensor  band¬ 
width.  With  ideal,  collocated  actuator/sensor 
pairs,  the  worst  cane  always  has  a  5p(w)>0. 

An  indication  thst  any  positive  controller  will 
preserve  closed-loop  stability,  independent  of 
modal  data  and  order  reduction. 

3.  Determine  the  "a-line"  that  will  reflect  de¬ 
sirable  pole  locations  (e.g.  based  on  desired 
transient  response)  and  compute  the  character¬ 
istic  gains  X(w')  for  the  a-shifted  system, 
i.e.,  evsluate  (6)  using  the-)!/*  operator  instead 
of  UJ. 

lA/’W  denotes  "substitute  — 7-^ 7-  for  z  in  {•}" 

r  (1-w') 


where 


r  -  e~aTSAMP 
0 


*  This  is  done  using  a  Nichols  chart  and  a  com- 
puter-sided  filter  design  program. 


4.  If  necessary,  shape  X (w’)  to  satisfy  the  Nyquist 
Theorem  over  the  bsndwidth  of  interest.  This 
ensures  that  modes  within  the  control  bandwidth 
have  a  desirable  rate  of  decay. 

Map  the  w1  controller  back  to  the  Z-domain. 

5.  Compute  Xc(w)  and  6c(w),  the  open-loop  X  and  5 
for  the  plant  plus  controller.  If  everything 
was  done  correctly,  Xc(w)  will  indicate  sta¬ 
bility  and  5c(w)  will  indicate  s  perturbation 
margin  that  encompasses  che  worst  case  plant. 

Note  chat  because  positivity  retains  the  con¬ 
vexity  properties  of  conicity,  satisfaction  of 
the  nominal  and  worst  case  indies  satisfaction 
of  all  convex  perturbations  in-between. 

6.  Simulate  to  verify  that  disturbance  rejection 
goals  have  been  satisfied.  It  should  be  very 
close.  Otherwise  the  gains  along  the  diagonal 
of  che  compensator  can  be  increased,  without 
loss  of  stability  (including  perturbations) 

by  an  amount  computable  from  6c(w)  and  6p(w). 

VI.  Tetrahedron  Example 

The  Tetrahedral  Truss,  shown  in  Figure  6  will 
now  be  used  to  illustrate  the  ideas  presented. 

This  structure  was  devised  by  Draper  Labs  [10]  as 
one  of  the  simplest  non-planar  geometries  capable 
of  representing  a  large  space  structure.  Despite 
its  apparent  simplicity,  it  models  the  feed-tower 
in  a  generic  class  of  large  antenna  applications. 
The  tetrahedral  apex  represents  the  antenna  feed, 
members  1-6  are  part  of  the  feed  support  structure, 
and  bi-peds  7-0,  9-10  and  11-12  are  supports/con¬ 
trols  attached  to  an  inertially  stabilized  (assumed) 
antenna  dish.  Because  of  the  nature  of  the  resl 
problem,  che  motion  of  the  feed  (apex)  must  be  con¬ 
trolled  through  actuators  and  sensors  placed  at  the 
bi-peds,  none  can  be  placed  at  the  feed.  The  lat¬ 
ter  constitutes  the  control  problem:  i.e.,  control 
the  feed  precisely  despite  large  modeling  uncer¬ 
tainty  and  without  resorting  to  direct  measurement 
or  control  of  its  motion. 


D-6 


(7)  LEO  1  STATION  1 
(S)  LEO  1  STATION  2 
®  LEO  2  STATION  1 
®  LEO  2  STATION  2 
(fi)  LEO  1  STATION  1 
©  LEO  2  STATION  2 


Six  colocated  sensors/accuacors  were  assumed 
to  act  as  member  dampers  in  parallel  with  members  7 
through  12.  The  member  dampers  were  assumed  capable 
of  sensing  relative  position  and  velocity  and  exert¬ 
ing  force  in  the  axi„l  direction  only. 

The  objective  was  to  design  a  robust  controller 
and  evaluate  the  resulting  controller  on  the  full 
order  nominal  model.  A  second  evaluation  was  then 
conducted  using  the  true  tetrahedral  system  param¬ 
eters  shown  in  Table  2. 

Table  2.  Tetrahedron  Parameter  Variation 


•nt  Property  Changes 


Figure  6.  Draper  Tetrahedral  Truss 

The  finite  element  method  was  used  to  model  the 
tetrahedron.  Referring  to  Figure  6,  the  model  con¬ 
tains  ten  nodes,  each  with  three  degrees  of  freedom, 
an4  twelve  truss  members.  Node  coordinates  and 
element  cennecclvenass  for  the  "nominal"  model  are 
listed  in  Table  1.  These  elements  are  assumed  to 
be  capable  of  resisting  axial  forces  only.  The 
masses  were  assumed  to  be  lumped  at  the  nodes  one 
through  four. 

Table  1.  Tetrahedron  Structural  Model 
Parameters  (Nominal) 

Node  Paint  Ucitioni 


X 

V 

I 

0.0 

0.0 

10.155 

-5.0 

-z.  sty 

2.0 

5.0 

-2.887 

2.0 

0.0 

5.7735 

2.0 

-4.0 

-1.1547 

0.0 

-4.0 

-4.6138 

0.0 

4.0 

-4.5188 

0.0 

5.0 

-1.1547 

0.0 

Z.O 

5.7735 

0.0 

-Z.0 

5.7735 

0.0 

Pfltnlnt*  X.r,  4  Z  It  nodi!  !-10 
Flfinint  Conntet1*1t1t!  md  Prooirtlii 


Both  evaluations  were  based  on  the  initial  condi¬ 
tions  given  in  Table  3.  The  goal  was  to  damp  the 
line-of-sight  motion  of  Node  1  in  the  x  and  y  direc¬ 
tions  to  less  than  0.0004  and  0.00025  units, 
respectively,  in  20  seconds. 

Table  3.  Initial  Conditions  for  Control  Design 
Evaluations 

»od»  Displacement  (n)  Velocity!") 


8 

M 

O 

O 

004 

9 

.0 

.0 

10 

.0 

.0 

11 

.0 

.0 

12 

.0 

.0 

Peterlel  E  •  1 .0 
Lunged  Mints 


The  sampling  rate  was  chosen  to  be  0.2  seconds. 
This  race  is  approximately  three  times  faster  than 
the  rate  required  to  reconstruct  the  highest  con¬ 
trolled  mode,  0.909  Hz.  [Note  0.909  Hz  is  not  Che 
highest  mode  modeled.  There  are  several  modes  at 
higher  frequencies  which  even  get  aliased.]  Using 
0.2  sec,  one  can  compute  the  characteristic  gains 
and  positivity  index  for  the  nominal  plant.  These 
are  shown  in  Figures  7  and  8. 


Figure  7.  Characteristic  Gains  of  the 
Tetrahedron  in  the  W-Plane 


rtT*«»*o»w 
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Figure  8.  Positivity  Index  of  Uncompensated 
System  (W-Plane) 

The  disturbances  for  this  system  were  not 
specified  in  [10]  so  one  can  go  directly  into  the 
a-shifted  plane  in  order  to  design  for  the  required 
rate  of  decay.  It  can  be  computed  that  a  «  0.230/sec 
will  decay  the  specified  initial  conditions  to 
levels  desired  after  20  seconds.  So  after  adding  a 
little  feedback,  7.8  •  I  to  a-priori  shift  the 
majority  of  poles  inside  the  a  region,  the  a-shifted 
characteristic  gains  are  computed.  Figure  9  shows 
that  in  order  to  obtain  the  required  rate  of  decay, 
modes  '/I  and  72  must  encircle  the  point  (180°,  0  db). 
The  reason  for  requiring  the  encirclement  is  that 
these  modes  are  open-loop  "unstable"  in  the 
a-shifted  plane,  so  the  encirclement  condition  of 
the  Nyquist  Theorem  applies.  Note  that  mode  #12  is 
also  unstable  in  the  a-plane  (stable  in  the  regular 
plane),  but  does  not  have  to  encircle  the  (-180 
0  db;  point  because  it  lies  outside  the  control 
bandwidth.  The  compensator 

I  •  11  492  -  3 ■  9211622  -  3.9863812 


Figure  9.  Characteristic  Cains  of  the  Tetrahedron 

in  the  a-Shifted  W-Plane,  R  -  0.955 

o 

Figures  10  and  11  show  the  characteristic  gains 
and  the  positivity  index  for  the  open-loop  plant 
plus  the  compensator.  The  characteristic  gains 
show  that  the  system  is  stable.  The  positivity 
index  shows  that  large  perturbation  are  permitted 
at  low  and  high  frequencies. 


2 

Z  -  1. 709212  +  .846578 


can  be  shown  to  shape  the  characteristic  gains  in 
the  desired  manner. 


Figure  10.  Frequency  Response  of  Compensated 
Open-Loop  System  in  the  W-Plane 
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Figure  11.  Positivity  Index  of  Compensated 
System  (W-Plane) 
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At  intermediate  frequencies,  the  system  remains 
relatively  insensitive  to  the  modal  data  (additive 
unmodeled  modal  perturbation),  but,  multiplicative 
phase  shifts  could  destabilize  the  system.  The  per¬ 
turbed  tetrahedron,  however,  primarily  consists  of 
changes  in  the  modal  data  and  it  can  be  computed 
that  the  perturbed  model  plus  all  convex  combina- 
tlons  in  between  are  well  within  the  permissible 
perturbation  region. 

Figures  12  through  17  show  the  simulated  trans¬ 
ient  responses  of  the  nominal  and  perturbed  systems 
to  the  Initial  conditions  given  in  Table  3. 

Figures  12  and  13  show  that  the  unmodelled  modes 
are  stable  despite  the  use  of  the  low  order  de¬ 
centralized  controller.  Figures  14  and  15  show 
that  the  nominal  system  meets  design  specifica¬ 
tions  but  that  there  is  a  small  loss  of  perfor¬ 
mance  using  the  perturbed  model.  Figures  16  and 
17  show  that  actuator  caiman d  forces  are  modest. 
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Figure  15.  Transient  Response  of  Perturbed 
System  to  Initial  Conditions 


Figure  13.  Response  of  Residual  Modes 
(Perturbed  System) 
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Figure  17.  Actuator  Force  at  Leg  1, 

Station  1  (Perturbed  System) 


VII.  Summary /Future  Work 

This  paper  has  summarized  those  aspects  of 
Positivity  and  Frequency  Domain  Theory  which  are 
useful  in  Che  design  of  LSS  control  systems.  A 
design  procedure  combining  the  two  was  described. 

An  example  illustrated  the  approach  and  good 
results  were  obtained. 

A  future  paper  will  further  detail  how  the 
basic  Ideas  have  been  used  to  design  the  control 
system  for  a  much  more  realistic  structure.  Future 
efforts  will  be  directed  at  applying  the  concepts 
to  "non-square"  LSS  plants  where  colocation  of 
actuators  and  sensor  is  not  permitted. 


Figure  14.  Transient  Response  of  Nominal 
System  to  Initial  Conditions 
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Appendix 

Ic  this  Appendix  the  stability  equivalence  of 
the  system  in  Figures  1  and  3  are  established  for 
linear  systems.  The  proof  relies  on  the  following 
two  identities:  s 

1)  If  I  +  A  is  invertible,  then 

I  ;  A  (I  +A)"1  -  (I+A)'1 

Proof:  (I+A)  +  A  2  I 

so  pre-operating  on  both  sides  of  the  identity  by 

(I+A)’1 

yields  the  desired  result. 


2)  If  G  and  rt  are  linear*,  and,  I+HC  and 
I+GH  are  invertible,  then 

G(I+HG)-1  -  (I+GH) _1  G 

& 

Proof:  G+CH/5  G+GHG 

therefore, 

G(I+HG)  -  (I+GH)  G 

so  pre-operating  on  both  sides  by  (I+HG)'1  yields 

G  -  (I+GH)  G  (I+HG)”1 
and  post-operating  on  both  sides  yields 

(I+GH)-1  G  -  G  (I+HG)’1. 

M»*rt ion:  The  stability  of  either  of  the  systems 
in  Figure  1  and  3  imply  the  stability  of  the  other 

a)  (I+CH)  and  [l+(G-C)  D]  have  inverses, 

b)  G,H,C,D,R,  and  S  are  stable, 
and 

c)  The  system  in  both  figures  are  controllable 
and  observable. 

Proof:  From  Figure  3: 

S  -  S  [I+(G-C)D]_1  (G-C)  R"1 
H  -  R  [(I+HC)-1  H-D]  S-1 

Therefore  the  transfer  function  from  U'  to  Y'  is: 
[nci]*1  c  ■ 

{i+s[i+(c-c)prl  <e-c>  [(i+ec)'1  H-Djs-1}*1  i  .  ■ 

s{l+(l+(C-C)D]*‘  (C-C)  [(I*HC)-‘  H-D]}'1  S_1  C  - 
s{l-[l+<G-C)Dr‘  (O-C)D  +{l+<e-C>Drl  (C-C)  (I+HC)”1  h}'1  S'1  a  - 
sftl+lS-OD]-1  +[l*(C-C)Drl  (C-C)(1+«C)*1h}  1  S’1  C  - 

StPMlM)'1!]  1  [l+(5-C)D]  S'1  a  - 
StHCSd+CH)'1  -CHO+OI)*1]  1  [l+<C-C)D]  s'1  a  - 

s[(i+a)  (i+cH)*1]  1  [ i+(c-od]  s'1  a  - 
S(I+CS)(l+CH)_l  [1*(C-C)  .J  s’1  a  - 
s<i-mm>U+ch>*1<o-c>  h'1  .  [140a]*1  c 

*  Linearity  is  not  really  necessary.  All  that  is 
needed  is  that  the  distributive  law  used  be  per¬ 
mitted.  This  fact  can  and  has  been  used  to  derive 
the  nonlinear  error  equation  in  a  reduced-order 
model-reference  adaptive  control  [6]. 
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